”%>-A138  108 

UNCLASSIFIE 

INVESTIGATION  OF  SHEATH  Pin 
GLOW  DISCHARGESCU)  AIR  FOP 
WRIGHT-PATTERSON  AFB  OH  SC 
D  DEC  83  AFIT/DS/PH/83-6 

ENOMENA  IN  ELECTRONEGATIVE 

CE  INST  OF  TECH 

HOOL  OF  ENGINEERING  G  L  DUK 
F/G  20/3 

1/3 

NL 

s 

■ 

■ 

■ 

■ 

1 

k. 

« 


HQ. 


CC 


s' .  t  j?** 

toll! 


■  r:t  0 

M* 

'l  l.  00 

-i' 


icjg!  g 


■•  ,  J#  .  ,•  $ 

■  '‘'rltfiw? 


toil  4  .  differ 


►*  .. 

f-T- 
.  J •,; 

■  T£  -*7\  7r  ■■  ••’ 

3&fc; 

#  * 


t* 


■  %  *  ' 
■  i  * 


•  •  ***£:'' 

'  ■  >*H>  f 
*  *  ‘*4- 
Ni*;  #** 

&  -v  a" ; 

f  i ~2* :  S 

ife«a 


^  *  ' 


AIR  UNIVERSITY 
UNITED  STATES  AIR  FORCE 


INVESTIGATION  OF  SHEATH  PHENOMENA  IN 
ELECTRONEGATIVE  GLOW  DISCHARGES 


DISSERTATION 


AF IT/DS/PH/83-6 


Gary  L.  Duke 
Captain  USAF 


DISTRIBUTION  STATE- MONT  A  - 

Approved  toi  public  ii'looaisj  * 
Di5ttibuHor<  Ur  huited 


WRIGHT-PATTERSON  AIR  FORCE  RASI,  OHIO 


nfir' 

SCHOOL  OF  ENGINEERING  ^electei 

Tk  FE6  2  2  S&4 


r 


i 


AFIT/DS/PH/83-6 


INVESTIGATION  OF  SHEATH  PHENOMENA  IN 
ELECTRONEGATIVE  GLOW  DISCHARGES 

DISSERTATION 

AFIT/DS/PH/83-6  Gary  L.  Duke 

Captain  USAF 


DTIC 

ELECTEB* 
FEB  2  2  1984  jl 

B 


Approved  for  public  release;  distribution  unlimited 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  OECLASSIFICATION/DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

6a.  NAME  OF  PERFORMING  ORGANIZATION 

5b.  OFFICE  SYMBOL 

Air  Force  Institute  of 

(If  applicable) 

Technology 

6c.  ADDRESS  (City.  State  and  ZIP  Code) 

AFIT/EN 

Wright-Patterson  AFB,  OH  45433 

8a.  NAME  OF  FUNOING/SPONSORING 

8b.  OFFICE  SYMBOL 

ORGANIZATION 

(If  applicable) 

AF  Aero  Propulsion  Laboratory 

8c.  ADORESS  (City,  State  and  ZIP  Code ) 

AFWAL/P00C-3 

Wright-Patterson  AFB,  OH  45433 

11.  Ti'TLE  (Include  Security  Classification) 

Investigation  of  Sheath  Phenomena  in  Electronefe 

1b.  RESTRICTIVE  MARKINGS 

NONE 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 


7b.  ADDRESS  (City,  State  and  ZIP  Code) 


B.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


12.  PERSONAL  AUTHOR(S) 

Duke,  Gary  L.,  Capt,  USAF 


13a.  TYPE  OF  REPORT  13b.  TIME  COVEREO 

Dissertation  - 


16.  SUPPLEMENTARY  NOTATION 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 
ELEMENT  NO. 


61102F 


1*.  DATE  OF  REPORT  (Yr„  Mo..  Day) 

1983  December 


TASK 

WORK  UNIT 

NO. 

NO. 

S2 

12 

16.  PAGE  COUNT 
218 


Dean  let  l-.-.f. .  .  ‘  1;  cr.-l  Prcfe-'lonal  C,»al©P*aal 


17. 


FIELD 


COSATI  COOES 


GROUP  SUB.  GR. 


9 


IB  SUBJECT  TERMS  (Continue  on  r everee  if  ntceetaryVHeiidemtHyeba  MBBMMmBWPS 

Anode  Fall,  Argon,  Attachment,  Boltzmann  Transport  Equation, 
Cathode  Fall,  Cathode  Fall  Voltage,  Cathode  Fall  Length, 
Drift  Velocity.  Electron  Kinetics,  Electronegative  Gas, 


IB.  ABSTRACT  (Continue  on  reverie  if  neceuary  and  identify  by  block  number) 

Two  different  methods  of  analyzing  the  cathode  fall  region  of  low  pressure  glow 
discharges  were  developed  and  applied  to  three  different  electronegative  gas  mixtures.  One 
method  was  based  on  a  self-consistent  numerical  solution  to  Poisson's  equation,  the  current 
continuity  equations  for  electrons  and  negative  ions,  and  the  current  conservation  equa¬ 
tion.  This  method  assumes  the  electrons  are  always  in  equilibrium  with  the  electric  field. 
The  other  method  was  based  on  a  self-consistent  numerical  solution  of  the  Boltzmann  trans¬ 
port  equation  for  electrons,  Poisson's  equation,  and  the  current  conservation  equation. 

This  method  allows  the  electrons  not  to  be  in  equilibrium  with  the  field.  Comparing  these 
two  methods  revealed  that  nonequilibrium  prevails  throughout  the  cathode  fall  region. 

The  electronegative  gas  mixtures  investigated  were  small  concentrations  (less  than 
10%)  of  hydrogen  chloride  in  helium,  argon,  or  xenon.  The  electric  field,  Townsend  ioni¬ 
zation  and  attachment  coefficients,  electron  and  negative  ion  current  densities,  and 
electron,  positive  ion,  and  negative  ion  number  densities  are  plotted  as  functions  of 


20.  OIBTRIBUTION/AVAILABILITY  of  abstract 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


UNCLABSIFIEO/UNLIMITED  □  SAME  A»  APT.  B  OTIC  USERS  □  I  UNCLASSIFIED 


22*.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Gary  L.  Duke,  Capt,  USAF 


22b  TELEPHONE  NUMBER 
(Include  Area  Code) 

(513)  255-2168 


22c.  OFFICE  SYMBOL 


AFWAL/POCE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


Item  18.  Equilibrium,  Glow  Discharge,  Helium,  Hydrogen  Chloride,  Ionization,  Negative  Ion, 
Nonequilibrium,  Poisson’s  Equation,  Self-consistent,  Similarity  Rules,  Sheath, 
Xenon 

Item  19.  distance  through  the  cathode  fall  region.  Discharge  current  densities,  cathode 

fall  lengths,  and  voltages  are  compared  to  other  theoretical  and  experimental  data 
through  a  scaling  relationship.  The  experimentally  observed  contraction  of  tie 
cathode  fall  length  in  electronegative  gases  with  helium  as  a  buffer  is  described 
and  not  predicted  for  argon  or  xenon  mixtures.  Other  regions  such  as  the  negative 
glow  and  anode  fall  are  also  briefly  discussed. 


UNCLASS  IF 


SECURITY  CLASSIFICATION  OF  THIS  PA' 


AFIT/DS/PH/83-6 


INVESTIGATION  OF  SHEATH  PHENOMENA  IN 
ELECTRONEGATIVE  GLOW  DISCHARGES 


DISSERTATION 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Doctor  of  Philosophy 


by 

Gary  L.  Duke,  B.A. ,  M.S. 
Capt  USAF 


Approved  for  public  release;  distribution  unlimited 


jnr»  L’tfIMV  vvivi\  vvvv 


AFIT/DS/PH/83-6 


INVESTIGATION  OF  SHEATH  PHENOMENA  IN 
ELECTRONEGATIVE  GLOW  DISCHARGES 


by 

Gary  L.  Duke,  B.A. ,  M.S. 
Captain  USAF 


Approved: 


Accepted: 


mgineering 


^ .  V\«B, 


Preface 


Glow  discharges  have  a  wide  variety  of  applications  such  as  in 
exclmer  lasers,  thyratrons,  or  In  plasma  reactors  for  the  deposition  of 
thin  films  or  for  plasma  etching.  Many  of  these  discharges  of  current 
interest  now  contain  reactive  electronegative  gases.  It  is  hoped  that 
this  study  has  provided  a  better  understanding  of  the  cathode  fall 
region  and  the  influence  of  attachment  within  this  region. 

A  number  of  people  have  provided  guidance  and  encouragement  to  me 
during  the  course  of  this  study.  In  particular,  thanks  are  given  to: 

Dr  Alan  Garscadden  for  his  patience  and  many  critical  discussions;  the 
other  members  of  my  committee,  Lt  Col  William  Bailey,  Maj  Michael 
Stamm,  and  Dr  John  Jones  for  Insuring  a  quality  and  administratively 
correct  document;  Colonels  George  Strand  (ret)  and  James  Johnson  for 
their  encouragement  and  the  time  they  gave  me  to  pursue  my  research  and 
to  put  this  document  together;  Marge  Phillips  and  Rita  Jackson,  two 
excellent  typists;  Bill  Bateson  for  assistance  In  plotting  the  graphs; 
and  especially  my  wife  Mary  and  sons  Brian  and  Eric,  who  helped  me  keep 
my  perspective  throughout  the  ordeal. 
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Abstract 

\  l  - 

-JTwo  different  methods  of  analyzing  the  cathode  fall  region  of  low 
pressure  glow  discharges  were  developed  and  applied  to  three  different 
electronegative  gas  mixtures.  One  method  was  based  on  a  self-consis¬ 
tent  numerical  solution  to  Poisson's  equation,  the  current  continuity 
equations  for  electrons  and  negative  ions,  and  the  current  conservation 
equation.  This  method  assumes  the  electrons  are  always  in  equilibrium 
with  the  electric  field.  The  other  method  was  based  on  a  self-consis¬ 
tent  numerical  solution  of  the  Boltzmann  transport  equation  for  elec¬ 
trons,  Poisson's  equation,  and  the  current  conservation  equation.  This 
method  allows  the  electrons  not  to  be  in  equilibrium  with  the  field. 
Comparing  these  two  methods  revealed  that  nonequilibrium  prevails 
throughout  the  cathode  fall  region. 

The  electronegative  gas  mixtures  investigated  were  small  concen¬ 
trations  (less  than  10%)  of  hydrogen  chloride  in  helium,  argon,  or 
xenon*.  The  electric  field,  Townsend  ionization  and  attachment  coeffi- 

A 

cients,  electron  and  negative  ion  current  densities,  and  electron, 
positive  ion,  and  negative  ion  number  densities  are  plotted  as  func¬ 
tions  of  distance  through  the  cathode  fall  region.  Discharge  current 
densities,  cathode  fall  lengths,  and  voltages  are  compared  to  other 
theoretical  ancf  experimental  data  through  a  scaling  relationship.  The 
experimentally  observed  contraction  of  the  cathode  fall  length  in 
electronegative  gases  with  helium  as  a  buffer  is  described  and  not 
predicted  for  argon  or  xenon  mixtures.  Other  regions  such  as  the 
negative  glow  and  anode  fall  are  also  briefly  discussed. 


Chapter  l_.  Introduction 

Introduction  and  Definition  of  Terms 

An  understanding  of  the  role  of  negative  ions  in  low  pressure  gas 
discharges  is  important  in  several  fields  such  as  electric  discharge 
lasers,  high  power  switches,  and  plasma  deposition  of  thin  films. 
Discharges  such  as  these  which  contain  negative  ion  densities  in 
addition  to  the  usual  electron  and  positive  ion  densities  are  termed 
electronegative  discharges.  They  are  so  named  since  the  attaching  gas 
yielding  negative  ions  is  referred  to  as  an  electronegative  gas. 
Discharges  containing  only  electrons  and  positive  ions  are  termed 
electropositive  discharges. 

Both  types  of  discharges  are  characterized  by  the  distinctive 
luminous  and  dark  regions  shown  in  Fig  1-1.  The  actual  position  of  the 
various  regions  depends  upon  the  gas  or  gas  mixture,  pressure,  and 
current  or  voltage  across  the  discharge.  Starting  at  the  cathode  there 
is  sometimes  visible  a  very  narrow  dark  space  termed  Aston's  dark  space 
followed  by  a  feeble  layer  of  light  called  the  cathode  glow.  However, 
the  term  cathode  glow  often  includes  the  total  region  from  the  cathode 
to  the  beginning  of  the  negative  glow.  The  region  between  the  cathode 
glow  and  negative  glow  is  not  totally  dark,  but  only  appears  so  to  the 
eye  in  comparison  to  its  luminous  neighboring  regions  and  is  termed  the 
cathode  dark  space.  The  negative  glow  is  probably  the  most  luminous 
region  of  the  discharge.  A  sharp  luminous  boundary  exists  with  the 
cathode  nonluminous  space,  but  the  boundary  on  the  anode  side  is 
diffuse  and  gradually  grows  dimmer  as  it  merges  into  the  relatively 
dark  Faraday  dark  space.  This  is  normally  followed  by  a  long  luminous 


Electric 


Voltage 


positive  column  region  extending  up  to  the  anode  glow.  The  anode  glow 
appears  adjacent  to  the  anode  surface.  Depending  on  the  distance 
between  the  electrodes,  only  Aston's  dark  space,  the  cathode  glow,  and 
the  cathode  dark  space  are  necessary  for  a  discharge  to  exist. 

It  has  been  observed  that  the  presence  of  negative  ions  can 
spatially  alter  the  above  regions  and  in  some  gases  such  as  iodine, 
negative  ions  can  become  the  dominant  negative  charge  carrier  through 
most  of  the  discharge  regions.  Formation  of  negative  ions  in  the 
positive  column  of  the  discharge  often  results  in  striations  which  can 
cause  fluctuations  in  the  light  output  of  the  positive  column.  The 
most  undesirable  effect  of  negative  ions  is  in  causing  severe  radial 
construction  of  the  positive  column  limiting  the  efficient  usage  of  the 
available  container  volume.  Up  to  now,  however,  the  exact  role  nega¬ 
tive  ions  play  In  the  cathode  fall  region  of  a  discharge  was  unknown. 
The  literature  frequently  cites  experimental  evidence  indicating  that 
the  cathode  fall  region  contracts  longitudinally  with  the  addition  of 
even  a  trace  of  an  attaching  gas.  This  was  interpreted  to  mean  that 
the  electronegative  gas  can  have  a  significant  impact  in  the  region. 

The  original  question  which  initiated  this  study  was  if  the  formation 
of  negative  ions  in  the  cathode  fall  led  to  the  longitudinal  contrac¬ 
tion  of  the  cathode  fall  region. 

Before  proceeding  further,  definitions  of  equilibrium  and  non¬ 
equilibrium  are  necessary  to  alleviate  any  confusion.  'Equilibrium'  in 
high  pressure  discharges  is  termed  local  thermodynamic  equilibrium 
(LTE)  and  means  that  all  states  of  the  charged  and  neutral  species  can 
be  described  by  one  temperature.  In  contrast  at  low  pressures,  the 
term  'equilibrium',  more  properly  called  'field  equilibrium',  means 


that  the  electron  energy  distribution  function  is  characterized  by  the 
local  value  of  E/N  (E  is  the  electric  field,  N  is  the  gas  number 
density)  and  the  species  in  the  volume  element  of  interest.  The 
electron  distribution  function  need  not  be  Maxwellian,  as  long  as 
it  can  be  defined  in  terms  of  E/N.  Nonequilibrium  at  low  pressures 
occurs  when  the  above  criterion  for  field  equilibrium  does  not  apply. 
For  example,  electrons  are  not  in  equilibrium  in  large  spatial  gradi¬ 
ents  where  the  mean  velocity  of  the  electrons  would  not  be  equal  to  the 
mean  velocity  of  electrons  in  an  equivalent,  but  spatially  uniform 
field. 

Basic  analytical  analysis  (ref  20,  21,  60,  79,  86)  now  indicates 
that  the  electrons  are  not  in  equilibrium  with  the  electric  field  in 
the  cathode  fall  region.  Since  the  electrons  are  being  accelerated 
across  the  cathode  fall  region,  they  are  gaining  energy  from  the  field 
faster  than  they  lose  energy  to  elastic  scattering  or  inelastic 
processes  in  the  gas.  The  electrons  leave  the  cathode  at  low  energies 
of  few  electron  volts  where  the  field  is  actually  the  highest.  The 
electrons  approach  the  negative  glow  with  a  multi-peaked  distribution 
function  as  the  field  is  approaching  a  minimum.  Since  the  electrons 
are  not  in  equilibrium  with  the  field,  their  multiplication  cannot  be 
described  by  the  usual  Townsend  ionization  coefficient.  This  coeffi¬ 
cient  is  measured  in  uniform  fields  where  the  energy  gained  by  the 
electrons  from  the  field  is  balanced  by  their  energy  losses  to  various 
elastic  and  Inelastic  processes  with  the  gas  molecules.  By  analogy, 
the  equilibrium  attachment  coefficient  also  is  not  expected  to  cor¬ 
rectly  describe  the  formation  of  negative  ions  across  the  cathode  fall 
region.  The  attachment  coefficient  is  expected  to  be  large  adjacent  to 


cathode  and  then  to  decrease  as  the  electrons  are  accelerated  to  higher 
energies.  The  formation  of  negative  Ions  In  this  nonequilibrium  sheath 
region  has  not  been  studied  either  analytically  or  experimentally. 
Concentrations  of  negative  ions  could  alter  the  space  charge  in  the 
sheath  region  and  thus  affect  the  voltage  drop  across  the  cathode  fall 
region.  This  phenomena  of  nonequilibrium  in  the  cathode  fall  region 
requires  a  'nonequilibrium  theory',  such  as  calculating  the  electron 
energy  distribution  function  as  a  function  of  distance  in  order  to 
properly  model  the  electron  number  density  and  ionization  and  attach¬ 
ment  coefficients.  However,  'equilibrium  theory'  has  often  been 
applied  to  describe  the  electric  field,  voltage  drop,  and  electron 
current  because  its  results  have  agreed  fortuitously  with  the  macro¬ 
scopic  parameters  that  have  been  able  to  be  measured  across  the  cathode 
fall  region. 

Objectives  of  This  Study 

This  analytical  study  investigates  the  role  small  concentrations 
of  an  attaching  gas  play  in  the  cathode  fall  region  of  a  glow  dis¬ 
charge.  To  do  this,  a  new  equilibrium  model  based  on  current  continu¬ 
ity  equations  and  Poisson's  equations  was  developed.  Additionally,  a 
nonequilibrium  model  based  on  the  Boltzmann  transport  equation  and 
Poisson's  equation  was  modified  to  be  self-consistent,  include  negative 
ions  and  be  more  flexible  The  electron  energy  distribution  function 
was  calculated  as  a  function  of  distance  through  the  cathode  fall 
region  and  the  effects  of  negative  ion  formation  were  examined.  The 
first  and  most  important  objective  was  to  examine  the  effects  in  the 
cathode  fall  region  resulting  from  the  addition  of  small  amounts  of  an 


attaching  gas  such  as  HC1  to  a  rare  gas  discharge.  The  second  objec¬ 
tive  was  to  model  the  nonequilibrium  behavior  of  the  three  charged 
particle  densities  through  the  cathode  fall  region.  The  third  objec¬ 
tive  was  to  compare  and  contrast  the  equilibrium  and  nonequilibrium 
approaches  and  determine  their  regions  of  applicability. 

Overview  of  Contents  and  Results 

Chapter  I  describes  the  purpose  and  objective  of  this  study, 
defines  the  regions  of  a  discharge,  and  outlines  the  contents  of  the 
remaining  chapters. 

In  Chapter  II,  the  sufficient  condition  for  the  existence  of  a 
discharge  is  described  in  terms  of  electron  multiplication  in  the 
cathode  fall  region.  A  historical  review  is  presented  on  how  trace 
amounts  of  negative  ions  affect  an  electropositive  discharge.  A 
qualitative  theory  is  described  which  simplistically  describes  the 
result  of  negative  ion  formation  in  glow  discharges. 

Chapter  III  critically  reviews  the  previous  equilibrium  theories 
of  the  cathode  fall  region,  primarily  for  electropositive  gases.  These 
equilibrium  theories  are  based  on  continuity  type  equations  coupled 
with  Poisson's  equation.  Boundary  conditions  are  also  discussed.  A 
new  equilibrium  model  of  the  cathode  fall  region  allowing  for  negative 
ion  formation  is  described  and  numerical  results  are  presented. 

In  comparison  to  the  equilibrium  approaches  described  in  Chapter 
III,  Chapter  IV  critically  reviews  the  previous  nonequilibrium  theories 
of  the  cathode  fall  region  based  on  calculating  the  electron  dis¬ 
tribution  as  a  function  of  position  through  the  cathode  fall  region.  A 
numerical  technique  for  calculating  the  electron  distribution  as  a 


function  of  position  in  the  cathode  fall  is  described.  This  study 
extends  this  technique  to  include  attachment,  the  presence  of  negative 
ions  in  the  relevant  equations,  gas  mixtures  with  a  variable  number  of 
collision  processes,  and  the  calculation  of  the  Townsend  ionization  and 
attachment  rates  from  the  distribution  function.  A  convergent  tech¬ 
nique  was  developed  which  leads  to  a  self-consistent  solution  between 
Poisson's  equation  and  the  electron  distribution  function.  The  boun¬ 
dary  conditions,  electron  and  ion  kinetics,  and  the  self-consistency  of 
the  solution  are  discussed.  Characteristic  cathode  fall  parameters  are 
calculated  as  a  function  of  distance  and  percentage  of  attaching 
species. 

Chapter  V  compares  the  equilibrium  theory  with  the  nonequilibrium 
theory  and  discusses  the  discrepancies  between  the  two  theories  in  the 
cathode  fall  region.  The  anode  fall  region  is  also  discussed. 

Finally,  Chapter  VI  summarizes  the  new  results  of  this  study  and 
presents  some  considerations  for  future  study. 


Chapter  II.  Review  of 


Qua! itative  Theory 

Description  of  the  Cathode  Fall  Region 

The  regions  of  a  discharge  were  identified  in  Fig  1-1.  Of  these, 
the  cathode  fall  region  is  the  most  important  segment  of  the  discharge 
in  that  the  largest  potential  drop  per  unit  distance  usually  occurs 
there.  However,  it  is  probably  the  region  least  understood  theoret¬ 
ically.  It  is  the  region  of  the  discharge  closest  to  the  cathode  and 
consists  of  the  three  subregions  identified  in  the  previous  chapter  as 
the  primary  dark  space  or  Aston's  dark  space,  the  cathode  glow,  and  the 
cathode  dark  space.  The  properties  of  the  cathode  fall  region  are 
almost  independent  of  the  rest  of  the  discharge,  since  a  discharge  can 
exist  without  a  positive  column,  Faraday  dark  space,  or  negative  glow. 

A  discharge  cannot  exist,  however,  without  a  cathode  dark  space  (ref 
6).  Specifically,  the  length  of  the  cathode  dark  space,  the  voltage 
drop  across  this  distance,  and  the  discharge  current  do  not  depend  on 
any  other  external  physical  parameters  such  as  the  length  of  the 
discharge  (providing  the  anode  is  not  so  small  as  to  inhibit  the  axial 
current).  The  cathode  dark  space  is  the  longest  of  the  three  cathode 
fall  regions. 

The  essential  processes  in  the  cathode  fall  region  are  as  follows. 
The  purpose  of  the  cathode  is  to  provide  electrons  to  maintain  the 
discharge.  Several  processes  exist  by  which  the  cathode  can  emit 
electrons.  Direct  processes  include  thermionic  emission  and  field 
emission.  For  applied  voltages  up  to  several  thousand  volts  these 
processes  can  be  ignored.  Indirect  emission  of  electrons,  usually 


called  secondary  emission,  results  from  bombardment  of  positive  ions, 
photons,  excited  and  neutral  molecules  on  the  cathode.  Due  to  the  low 
collision  rate  of  excited  and  neutral  molecules  with  the  cathode  in 
most  discharges,  the  physical  processes  contributing  most  to  the 
emission  of  electrons  from  the  cathode  are  ion  and  photon  bombardment. 
The  electrons  produced  by  secondary  emission  are  accelerated  across  the 
cathode  dark  space,  whereupon  they  begin  to  excite  and  ionize  the  gas 
creating  new  electrons  and  positive  ions.  The  positive  ions  in  turn 
are  accelerated  back  towards  the  cathode  where  they  produce  new  elec¬ 
trons.  For  every  (M-l)  electrons  produced  through  ionization  in  the 
gas,  y  electrons  are  released  at  the  cathode.  The  maintenance  condi¬ 
tion  for  a  self  sustained  discharge  becomes 

m-D*i  ot  m  ■  (ii-D 

The  multiplication  M  of  electrons  depends  upon  the  nature  and  pressure 
of  the  gas.  It  is  commonly  expressed  (in  an  electropositive  gas)  as 

/A*  e.xp  ( ii-2) 

where  d(x)  is  Townsend's  ionization  coefficient.  The  number  of  elec¬ 
trons  released  at  the  cathode  per  incident  particle  (y)  is  influenced 
by  the  size  and  energy  of  the  ion  or  photon  as  well  as  the  cathode 
material  and  the  state  of  its  surface.  Thus,  the  cathode  fall  region 
is  where  processes  essential  to  the  maintenance  of  the  rest  of  the 
discharge  occur.  The  next  section  describes  the  scaling  relationships 
for  glow  discharges. 
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Although  the  various  parameters  describing  a  glow  discharge  are 
all  interrelated  and  often  depend  on  each  other  through  complex  rela¬ 
tionships,  the  parameters  describing  two  different  discharges  can  be 
related  to  each  other  using  similarity  or  scaling  rules.  Two  dis¬ 
charges  are  said  to  be  similar  if  they  are  identical  in  one  or  two 
parameters  and  different  in  others.  For  example,  two  discharges  in  the 
same  gas,  at  the  same  gas  temperature,  with  the  same  electrode  materi¬ 
al,  and  with  the  same  potential  difference  between  the  electrodes  are 
similar  even  if  all  the  corresponding  linear  dimensions  differ  by  a 
constant  factor.  This  includes  the  vessel  dimensions,  the  electrode 
dimensions,  and  the  properties  of  the  gas  such  as  the  mean  free  path. 
Francis  (ref  6)  has  comprehensively  derived  a  set  of  similarity  rules 
which  result  in  several  invariant  quantities.  If  the  discharge 
voltage,  current,  and  gas  temperature  are  the  same  for  both  discharges, 
he  has  shown  that  the  following  quantities  are  also  invariant: 

E/p  -  electric  field/pressure 

pd  -  pressure  *  distance 

J/pz  -  current  density/pressure2 

These  parameters  remain  invariant  for  respective  points  in  the  dis¬ 
charges  as  long  as  the  following  processes  dominate  the  kinetics. 

a)  Single  stage  processes  such  as  electron  impact  ionization, 
attachment,  or  detachment. 

b)  Penning  ionization 

c)  Charge  transfer 

d)  Drift  and  diffusion  from  a  volume 

If  two  stage  processes,  photoionization,  or  recombination  are  the 


dominant  processes  in  the  discharge,  then  the  similarity  rules  will  not 
be  followed  and  the  scaling  parameters  above  may  no  longer  be  appro¬ 
priate.  The  above  invariant  quantities  apply  both  to  the  whole 
discharge  as  well  as  to  specific  regions  of  the  discharge.  For  exam¬ 
ple,  given  two  discharges  with  the  same  discharge  voltage,  current,  and 
gas  temperature,  the  invariant  quantities  for  the  cathode  fall  region 
would  be: 

Vc  -  cathode  fall  voltage  drop 

E/p  -  ratio  of  electric  field  to  pressure  at  respective 
positions  in  the  cathode  fall 
pdc  -  pressure  *  cathode  fall  length 

J/p2  -  current  density/(pressure)2 

Basically,  given  V,  I,  and  T,  then  Vc>  E/p,  pdc,  and  J/p2  are 
constants  and  relate  how  the  electric  field,  cathode  fall  length  and 
current  density  scale  with  pressure.  The  next  section  describes  the 
experimental  effects  observed  and  attributed  to  negative  ions  in  low 
pressure  glow  discharges  including  the  cathode  fall  region. 

Effects  Attributed  to  Negative  Ions  in  Gas  Discharges 

Several  authors  have  reviewed  the  effects  of  adding  small  amounts 
of  negative  ions  to  a  glow  discharge  (ref  5,  6,  8),  but  all  agree  with 
the  initial  description  advanced  by  Emeleus  and  Sayers  in  1938  (ref 
42).  Emeleus,  et  al.  observed  the  effects  on  the  discharge  when  a 
trace  of  an  attaching  gas  such  as  chlorine  was  added  to  a  neon  dis¬ 
charge  containing  a  little  helium.  With  the  addition  of  a  trace  of 
chlorine,  the  fluorescence  from  the  pure  neon  discharge  illustrated  in 
Fig  1 1 - 1  a  changed  to  that  illustrated  in  Fig  Il-lb.  The  following 


results  were  observed  at  a  pressure  of  2  torr  and  discharge  potential 


of  800v: 
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Fig  II-l  Effect  of  Adding  an  Electronegative  Gas  to  a 

Rare  Gas  Discharge 


a.  All  the  negative  zones  (cathode  dark  space,  4;  cathode  glow, 
3;  negative  glow,  5;  and  Faraday  dark  space,  6)  contracted  towards  the 
cathode  such  that  Aston's  dark  space,  or  the  primary  dark  space,  2, 
became  so  small  as  to  be  unobservable. 

b.  Curved  striations  appeared  in  the  positive  column,  8,  often 
accompanied  by  an  increase  in  the  potential  gradient. 

c.  A  new  dark  sheath,  7,  appeared  between  the  Faraday  dark 
space,  6,  which  was  moderately  luminous  and  the  positive  column,  8. 

d.  A  new  dark  sheath,  9,  also  appeared  between  the  positive 
column,  8,  and  anode  glow,  10. 
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Spencer-Smith  (ref  81)  and  more  recently,  Davis  and  King  (ref  34) 
observed  similar  but  more  pronounced  effects  in  pure  iodine  discharges, 
and  in  helium  containing  a  trace  of  iodine.  The  positive  column  was 
observed  to  constrict  radially  to  either  a  stationary  form  or  a  mobile 
ribbon.  This  type  of  constricted  positive  column  has  often  been  found 
in  other  discharges  containing  electronegative  gases. 

Following  the  discussion  by  Emeleus  and  Sayers  (ref  39),  the 
origin  of  the  characteristic  effects  of  electronegative  gases  was 
attributed  to  the  decrease  in  the  average  mobility  of  the  negative 
charge  carriers  as  electrons  are  converted  to  heavier  negative  ions. 
Three  almost  independent  effects  ensue: 

a.  Creation  or  addition  of  negative  ions  corresponds  to  an 
increase  in  pressure  in  an  electropositive  discharge.  The  reduction  of 
the  drift  velocity  of  the  electrons  by  their  attachment  to  gas  mole¬ 
cules  is  equivalent  to  the  more  continuous  retardation  caused  by  more 
frequent  collisions  from  an  increase  in  gas  pressure.  The  contraction 
of  the  negative  zones,  the  occurrence  of  striations  in  the  positive 
column,  and  the  constriction  of  the  positive  column  are  all  charac¬ 
teristics  of  an  electropositive  discharge  at  higher  pressures,  for 
example  100  torr  compared  to  1  torr. 

b.  Since  negative  ions  are  more  massive  than  electrons,  they 
diffuse  more  slowly  from  regions  of  space  charge.  Thus  regions  of 
space  charge  tend  to  grow  once  negative  ions  have  formed.  This  leads 
to  the  formation  of  double  space  charge  layers  of  positive  ions  next  to 
negative  ions.  The  high  field  gradient  between  these  layers  enhances 
the  formation  of  stationary  striations  and  thus  the  enhanced  boundaries 
between  dark  and  luminous  regions.  The  axial  field  resulting  from 


these  stationary  double  (+-)  space  charge  regions  accelerates  electrons 
to  sufficient  energy  to  produce  maximum  excitation  of  the  gas  atoms  and 
molecules  leading  to  increased  light  emission.  The  dark  sheaths,  7, 
between  the  positive  column,  8,  and  the  Faraday  dark  space,  6,  and  the 
anode  glow,  9,  are  believed  to  be  the  result  of  space  charge  due  to 
negative  ions  created  in  previous  regions  of  the  tube  (ref  5,  6,  42). 

c.  Since  some  of  the  electrons  become  attached  to  gas  molecules, 
the  radial  distribution  of  negative  charge  carriers  will  change.  Due 
to  the  lower  negative  ion  mobility,  the  negative  ions  will  be  less 
likely  to  diffuse  away  from  the  slightly  positive  axis  in  the  positive 
column.  Since  detachment  of  electrons  from  negative  ions  can  often  be 
an  inefficient  process  in  many  gases,  large  concentrations  of  negative 
ions  can  build  up  near  the  axis.  This  causes  positive  ions  also  to 
drift  to  the  axis  to  balance  the  total  charge  resulting  in  a  con¬ 
stricted  discharge. 

The  above  description  is  based  mainly  on  qualitative  observations 
which  have  been  reviewed  in  reference  5.  Several  quantitative  measure¬ 
ments  have  been  made  most  of  which  have  been  accomplished  since  1938, 
confirming  the  above  hypotheses.  These  will  be  reviewed  briefly. 

Spencer-Smith  (ref  81)  and  Woolsey  (ref  96)  used  probe  measure¬ 
ments  and  magnetic  methods  for  measuring  charge  to  mass  ratios  to 
investigate  the  negative  glow  and  positive  column  regions  of  low 
pressure  iodine  discharge.  Probe  measurements  in  the  presence  of  large 
quantities  of  negative  ions  are  very  difficult  to  interpret,  but 
Spencer-Smith  found  the  positive  and  negative  ion  concentration  in  the 
negative  glow  to  be  about  equal,  and  about  100  times  larger  than  the 
electron  concentration.  In  the  positive  column,  the  ratio  of  positive 


ions  to  electrons  remained  about  the  same,  but  the  absolute  concentra¬ 
tions  were  lower.  Spencer-Smith  also  found  the  electron  temperature  in 
the  negative  glow  to  be  normal  (~1  to  1.5  eV)  for  E/N  =  620Td.  The 
electron  temperature  In  the  positive  column  was  much  greater,  ~13  eV. 
This  can  be  accounted  for  by  some  of  the  electrons  being  accelerated 
through  several  potential  drops  as  a  result  of  striations.  However, 
the  conventional  concept  of  electron  temperature  under  such  nonequi¬ 
librium  circumstances  probably  is  not  valid.  Woolsey's  results  were 
similar  to  Spencer-Smith's  in  that  he  found  the  ratio  of  negative  ions 
to  electrons  to  be  200  in  both  the  negative  glow  and  the  positive 
column.  Woolsey  also  determined  the  radial  distribution  of  negative 
charge  carriers  to  be  much  flatter  than  the  normal  zero  order  Bessel 
function  found  in  the  negative  glow.  The  number  of  negative  ions  was 
found  to  be  much  larger  and  more  concentrated  around  the  axis  in  the 
positive  column  than  the  negative  glow,  in  agreement  with  Emeleus  and 
Sayer's  speculation.  In  both  experiments  the  positive  column 
fluorescence  about  the  center  of  the  discharge  tube  was  constricted. 

In  a  similar  experiment  using  probes,  Zimmermann  (ref  98)  found  the 
positive  ion  concentrations  more  than  fifty  times  greater  than  the 
electron  concentration  in  the  constricted  positive  column  of  a  hot 
cathode  discharge  in  hydrogen  chloride  vapor. 

Lunt  and  Gregg  (ref  63),  and  Thompson  (ref  84),  have  investigated 
oxygen  discharges.  Again  using  probes  and  mass  analysis,  Lunt  and 
Gregg  found  the  concentration  of  0~  in  the  positive  column  with  an 
E/N  -  250  Td  to  be  about  equal  to  the  concentration  of  electrons.  In 
the  Faraday  dark  space  and  negative  glow,  however,  hardly  any  negative 
ions  were  detected  (<!%).  Thompson  used  a  radio  frequency  mass  spec- 


trometer  probe  to  measure  the  particle  concentrations.  He  found  the 
concentration  of  0"  was  about  equal  to  the  concentration  of  02+ ,  thus 
confirming  Lunt  and  Greg's  observation  that  0"  was  the  most  abundant 
negative  ion.  Thompson  also  measured  negative  particle  energy  dis¬ 
tributions  in  various  regions  of  the  discharge  and  compared  them  to  the 
electron  energy  distributions  in  the  same  regions  in  a  pure  nitrogen 
discharge.  The  energy  distributions  in  oxygen  contained  an  apparent 
large  low  energy  peak  due  to  negative  ions,  as  expected.  He  determined 
the  radial  space  potential  in  the  positive  column  to  be  much  flatter 
for  oxygen  than  for  nitrogen  and  the  axial  space  potential  changed  much 
more  abruptly  through  a  striation  than  in  nitrogen.  Although  both 
discharges  contained  standing  striations,  the  striations  in  the  oxygen 
discharge  were  much  closer  together.  The  radial  distribution  of  the 
electron  concentration  in  the  oxygen  discharge  was  almost  constant 
whereas  the  radial  variation  of  positive  and  negative  ions  both 
decreased  almost  to  zero  supporting  the  idea  that  the  positive  and 
negative  ions  would  concentrate  near  the  axis. 

More  recently  Davis  and  King  (ref  4}  observed  the  striations 
formed  in  a  He-I2  glow  discharge  laser.  They  calculated  the  ratio  of 
negative  ions  to  electrons  to  be  on  the  order  of  20-110  for  a  current 
density  of  100  amp/m2  and  on  the  order  of  2-11  for  a  current  density  of 

4 

10  amp/m2.  The  negative  ions  with  their  lower  mobility  readily  led  to 
the  formation  of  both  stable  and  unstable  stationary  and  moving 
striations  as  well  as  the  constriction  of  the  positive  column  when 
small  quantities  (<.1  torr)  of  iodine  was  admitted  to  pure  helium 
discharges  (2-10  torr).  The  forward  wave  (anode  to  cathode)  striations 


occurred  at  a  10kHz  frequency  and  results  in  a  partial  modulation  of 
the  laser  intensity.  The  radius  of  the  constricted  positive  column 
decreased  with  an  increase  in  iodine  pressure.  Even  when  the  column 
remained  stable,  it  was  not  always  axially  symetric,  but  often  adopted 
a  twisted  configuration  in  the  discharge  tube  caused  by  fluctuations  in 
the  local  wall  temperature  and  iodine  density.  Thus  constricted 
discharges  are  undesirable  for  laser  applications,  since  they  prevent 
efficient  usage  of  the  available  resonator  volume.  The  instabilities 
which  result  in  these  striations  and  constriction  of  the  positive 
column  has  been  investigated  extensively  by  Nighan  and  Wiegand  (ref 
70).  This  topic  will  not  be  investigated  in  this  study.  The  occur¬ 
rence  of  negative  ions  in  the  negative  glow  and  positive  column  is, 
however,  relevant  to  the  cathode  fall  region  and  the  boundary  con¬ 
ditions  relating  these  regions. 


Chapter  III.  Electron  Kinetics  in  the  Cathode 
Fall  Region:  Equil ibrium  Analysis 

Review  of  Equilibrium  Theories 

There  are  several  theories  for  the  cathode  fall  region  in  an 
electropositive  gas  which  relate  the  cathode  fall  voltage  (Vc),  the 
discharge  current  (J),  and  the  cathode  fall  distance  (dc).  These 
theories  will  be  referenced  later  as  they  are  discussed.  In  compari¬ 
son,  there  is  a  distinct  lack  of  any  kind  of  analytical  analysis  for 
electronegative  discharges.  These  theories  for  the  cathode  fall  region 
can  be  divided  into  equilibrium  analyses  based  on  Poisson's  equation 
and  hydrodynamic  current  continuity  equations,  and  nonequilibrium 
analyses  based  on  calculating  an  electron  distribution  function.  The 
equilibrium  analyses  will  be  critically  analyzed  and  expanded  upon  in 
this  chapter.  The  nonequilibrium  analyses  will  be  considered  in  the 
next  chapter. 

The  equilibrium  analyses  are  of  two  types;  those  which  regard  the 
cathode  dark  space  as  an  isolated  independent  unit  considered  separate 
from  the  negative  glow  and  those  which  include  the  negative  glow  so 
that  together  their  processes  sustain  the  discharge.  These  theories  of 
glow  discharges  are  thoroughly  summarized,  but  not  necessarily  criti¬ 
cally  reviewed  in  several  references  (ref  6,  7,  89).  The  cathode  fall 
current-voltage  relationship  derived  by  von  Engel  and  Steenbeck  (ref 
90)  for  a  self-sustaining  glow  discharge  has  provided  a  benchmark  for 
all  following  papers  and  will  be  briefly  reviewed  first.  The  emphasis 


will  then  be  on  results  presented  since  1975. 

The  von  Engel  and  Steenbeck  theory  assumed  the  electric  field  in 


the  cathode  fall  to  be  a  linearly  decreasing  function,  of  the  form 


E  on-') 


where  E  is  the  electric  field  at  the  cathode  and  d„  is  the  cathode 
o  c 

fall  distance.  The  ion  flow  to  the  cathode  is  assumed  to  be  mobility 
limited  so  that 

v,*^,  E  UII-2) 


where  v+  is  the  positive  ion  velocity  and  y+  their  mobility.  The 
boundary  conditions  that  von  Engel  and  Steenbeck  used  set  the  positive 
ion  current  to  zero  at  the  cathode  fall -negative  glow  boundary  and 
assumed  each  positive  ion  incident  on  the  cathode  produced  y  electrons, 
yielding 
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(II 1-3 ) 

jt(o)=a(o) 

( III -4 ) 

equation  ( 1 1 1-3 )  is  not  entirely  correct  in  that  the  positive  ion 
current  does  decrease  at  the  cathode  fall -negative  glow  boundary  but 
does  not  necessarily  to  to  zero.  The  question  of  whether  the  number  of 
positive  ions  entering  the  cathode  fall  from  the  negative  glow  is 
significant  or  not  has  developed  into  an  area  of  controversy  and  is 
described  later. 


Using  the  empirical  form  of  Townsend's  ionization  coefficient, 


(II 1—5) 


they  derived  equations  relating  the  total  current  J,  the  cathode  fall 
voltage  Vc,  and  the  cathode  fall  width  dc.  The  cathode  fall  thickness 
is  derived  by  applying  the  maintenance  condition,  equation  1 1 -1 ,  to  the 
cathode  fall  region  resulting  in 

J 

j€  -  Jx  (t*  y)  ( II 1-6) 

o 

Definition  of  the  width  of  the  cathode  fall  is  not  exact  but  it  is 
usually  defined  as  the  distance  from  the  cathode  to  the  point  in  the 
discharge  where  the  electric  field  is  extrapolated  to  zero.  The 
voltage  is  determined  from  the  electric  field. 

J 

V  **  f  E(k)  Jg  ft  ( 1 1 1-7) 

if  E(x)  is  given  by  equation  III-l.  Since  the  ion  current  to  the 
cathode  is  generally  mobility  limited,  they  derived  a  current-voltage 
relationship  for  the  cathode  fall  from  the  high  pressure 
space-charge-limited  current  voltage  relationship  yielding 

%xs  (l*t)  8  ^*P  (fa)*  (III -8 ) 

Equations  1 1 1-6 ,  III-7,  and  1 1 1-8  uniquely  describe  the  cathode  fall 
region  and  obey  the  similarity  or  scaling  rules  described  in  Chapter  I. 
Solving  these  equations  numerically  for  a  linearly  decreasing  electric 
field,  provides  good  agreement  with  experiment  over  a  wide  range  of 


conditions  in  electropositive  gases.  Helium  is  an  exception  and  one 
reason  may  be  that  the  Townsend  ionization  coefficient  cannot  be  fitted 
with  an  exponential  form  like  equation  1 1 1-5 .  However,  the  above 
assumptions  have  recently  become  controversial  areas  (ref  7). 

Druyvesteyn  and  Penning  (ref  38)  pointed  out  that  if  a  significant 
number  of  positive  ions  entered  the  cathode  fall  region  from  the 
negative  glow,  the  maintenance  condition  for  the  discharge  was  altered 
to  the  following: 


-  /  y  \*  H 

7T  (  II 1-9) 

where  6  is  the  ratio  of  the  positive  ion  current  to  the  electron 
current  at  the  cathode  fall-negative  glow  interface.  Thus  the  value  of  a 
d£  could  be  reduced  if  the  positive  ion  current  was  comparable  to  the 
electron  current  at  the  cathode  fall-negative  glow  boundary. 

Little  and  Von  Engel  (ref  59)  refute  the  idea  that  a  significant 
number  of  ions  enter  the  cathode  fall  region  from  the  negative  glow.  A 
linear  electric  field  requires  a  constant  net  space  charge  density  in 
the  cathode  fall,  whereas  if  a  significant  number  of  positive  ions 
originated  in  the  negative  glow,  then  the  positive  ion  density  would  be 
expected  to  decrease  towards  the  cathode.  As  Ingold  (ref  7)  pointed 
out,  these  are  not  conclusive  arguments,  since  such  field  measurements 
are  very  difficult  and  any  slowly  varying  field  would  appear  approxi¬ 
mately  linear  to  the  experimentalist.  Reviews  by  Francis  (ref  6)  and 
Weston  (ref  93)  survey  many  attempts  to  measure  the  electric  field  in 
the  cathode  fall  region.  The  weight  of  the  evidence  suggests  an 
electric  field  decreasing  linearly  through  the  cathode  dark  space 
becoming  almost  zero  in  the  negative  glow.  Several  recent  theoretical 


papers  have  tried  to  show  the  electric  field  to  be  linear  in  the 
cathode  fall.  These  papers  will  be  reviewed  later,  in  chronological 
order. 

The  other  area  of  controversy  is  the  concept  of  Townsend's  ioni¬ 
zation  coefficient  (°)  and  the  assumption  that  the  electrons  are  in 
equilibrium  with  the  field.  a  is  the  number  of  ionizing  collisions 
per  centimeter  of  path  length  and  is  usually  determined  from  drift  tube 
experiments.  These  experiments  assume  the  electrons  are  in  equilibrium 
with  a  uniform  electric  field.  This  condition  probably  is  not  met  in 
the  cathode  fall  where  the  electric  field  changes  rapidly.  Evidence  of 
this  appears  in  the  nonequilibrium  analyses  of  the  cathode  fall  by  Long 
(ref  60)  and  Tran  Ngoc  An  (ref  86)  which  are  discussed  in  Chapter  IV. 
Thus  the  theoretical  models  based  on  the  concept  that  the  electron 
distribution  function  is  in  equilibrium  with  the  field  may  not  be 
valid.  As  von  Engel  noted  about  his  own  theory,  the  equilibrium 
theories  give  "physically  the  correct  picture,  but  the  numerical 
agreement  with  observations  is  largely  fortuitous"  (ref  89:224).  This 
observation  is  also  discussed  later  in  Chapter  V  where  the  equilibrium 
and  nonequilibrium  results  are  compared. 

Returning  to  the  review  of  the  literature,  since  1975,  U1 'yanov 
(ref  87)  analytically  investigated  the  effect  of  the  positive  ion 
contribution  from  the  negative  glow.  He  based  his  approach  on  experi¬ 
mental  evidence  of  Gunthershultze  (ref  46)  indicating  that  as  the  anode 
approaches  the  plasma  side  of  the  cathode  sheath,  the  potential  drop 
across  the  cathode  fall  increases.  This  indicates  that  some  of  the 
positive  Ions  required  to  carry  the  current  in  the  cathode  fall  region 
originate  in  the  negative  glow.  U1 'yanov  solved  the  following  set  of 


equations  for  both  high  and  low  ion  mobility  cases,  defined  below 
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JtsC.n^E  (low  field  case) 

( 1 1 1-1 3a ) 

L 

J+scn.k,EJ  (high  field  case) 

( 111-1 3b) 

The  boundary  conditions  which  neglect  any  potential  drop  in  the  nega¬ 
tive  glow  are  given  by 


£(</)  =  0  (hi-'4) 

(III-15) 

where  6  was  a  variable  parameter  representing  in  his  theory  the  ratio 
of  positive  ion  current  density  to  electron  current  density  at  the 
cathode  fall -negative  glow  boundary.  For  both  mobility  cases,  his 
final  transcendental  equations  contained  integrals  which  had  to  be 
solved  numerically.  His  results  indicated  that  if  the  fraction  of  ion 
current  coming  from  the  negative  glow  (  6)  decreases,  then  the  electric 
field  at  the  cathode  increases.  For  a  given  secondary  emission  coeffi¬ 
cient  (  a),  he  determined  that  the  normal  cathode  fall  potential  (VN), 


normal  current  density  (J^),  and  normal  cathode  fall  length  (x^) 
changed  quite  slowly  over  a  wide  range  of  6  values.  He  found  the 
greatest  discrepancy  between  his  theory  and  von  Engle's  theory  at  small  6 
where  a  singularity  existed  in  his  equations  for  and  J^.  His 
equations  for  V^,  J^,  and  x^  depended  only  on  gas,  cathode  parameters, 
and  6,  once  the  transcendental  equation  for  the  electric  field  was 
solved.  From  the  fact  that  V^,  J^,  and  x^  were  a  function  of  6,  he 
concluded  that  a  complete  theory  for  the  cathode  fall  normal  glow  dis¬ 
charge  must  take  into  account  ion  production  in  the  negative  glow. 

Neuringer  (ref  69)  made  an  analytical  investigation  of  the  cathode 
fall  region  in  high-voltage  low-current  discharges.  He  determined  that 
to  first  order  the  electric  field  decreased  linearly  in  the  cathode 
fall  region,  the  voltage  drop  depended  on  the  two-thirds  power  of  J/p2 
and  the  cathode  fall  thickness  was  independent  of  the  operating  current 
and  varied  inversely  with  pressure.  His  model  combined  the  electron 
current  continuity  equation, 

< (E(*))Jt(x)  ( II 1-16) 


Poisson's  equation, 
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with  the  requirement  for  total  current  conservation. 


V'JsO  or  T*T++Te  *  constant 


( 1 1 1-18) 


J  and  J,  are  the  electron  and  positive  ion  current  densities  and  v 
e  +  e 

and  v+  are  the  electron  and  positive  ion  drift  velocities.  He  used 
Townsend's  ionization  coefficient  (equation  III-5),  Ward's  (ref  92) 
modified  formula  for  the  positive  ion  velocity  for  ions  in  a  high 
field. 


( III-19) 


and  the  fact  that  v+/vg  <<  1  throughout  the  cathode  fall  to  obtain  a 
differential  equation  separable  in  and  E.  This  equation 
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(III- 20 ) 


can  be  integrated  assuming  the  electric  field  goes  to  zero  and  the 
electron  current  density  is  equal  to  the  total  current  density  at  the 
cathode  fal 1 -negative  glow  boundary. 


EU)=0 


( 1 11-21 ) 
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(III-22) 


He  was  able  to  derive  an  expression  for  the  electric  field  as  a  func¬ 
tion  of  distance  by  assuming  Jg  «  0  close  to  the  cathode.  Keeping 
only  first  order  terms,  he  derived 


Assuming  no  externally  initiated  electron  current  density,  such  that 


( III -24 ) 


he  also  obtained 
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The  derivation  of  equation  1 1 1-23  is  important  because  it  was  the  first 
time  the  linear  functional  dependence  of  the  electric  field  on  distance 
had  been  derived,  and  not  assumed.  He  claims  good  agreement  with 
experimental  results  quoted  in  Francis'  article  (ref  6)  and  with  a 
computer  model  using  a  differential  equation  solver  technique  which  was 
previously  developed  by  Ward  (ref  92)  to  solve  equations  III-16  and 
1 1 1 - 1 7 .  Even  though  he  uses  valid  boundary  conditions  (equations 
1 1 1-21  and  1 1 1-22),  his  model  is  correct  only  in  region  very  near  the 
cathode  since  he  has  assumed  v+  ~  k+(E/p)*  and  «  J  in  order  to 
drop  terms.  Additional  analysis  in  the  next  section  includes  negative 
ions  and  reveals  that  the  first  order  approximations  do  not  have  to  be 
made.  Nevertheless,  his  results  add  credence  to  von  Engel's  assumption 
that  the  electric  field  decreases  linearly  in  the  cathode  fall  region. 

Davies  and  Evans  (ref  32)  extended  Neuringer's  analysis  and  showed 
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( 1 1 1 - 25 ) 


it  could  be  applied  to  the  whole  discharge  by  changing  the  boundary 
conditions  and  assuming  the  ratio  of  positive  ion  to  electron  mobility 
to  be  a  constant. 


This  assumption  was  found  to  breakdown  only  at  electric  fields  for 


which  1  -  Jg  «  rJg,  which  normally  was  not  the  case.  The  boundary 
condition  at  the  cathode  was  the  same  as  Neuringer's  i.e.  equation 
1 1 1-24 ,  but  the  other  boundary  condition  was  changed  to  represent  the 
electric  field  at  the  anode  or  to  specify  a  particular  applied  voltage. 


E6/)s£*  (II 1-27 ) 

Je(x)Jx  (II 1-28) 

Based  upon  combining  Poisson's  equation  1 1 1 -1 7  and  the  electron  current 
continuity  equation  III-16  into  a  transcendental  type  of  equation, 
including  the  above  changes,  they  developed  a  new,  fast  numerical 
technique  for  predicting  the  variation  of  E  and  j  as  a  function  of 
distance  from  cathode  to  anode.  Although  equation  III-16  is  true  in 
the  cathode  fall  region  where  ionization  is  the  dominant  process  and 
the  electrons  are  scattered  predominantly  in  the  forward  direction,  it 
is  not  correct  in  the  positive  column  region  where  the  electrons  are 
scattered  isotropically  and  the  gain  of  electrons  is  balanced  by  losses 
to  radial  diffusion,  recombination,  or  attachment.  In  essence,  when  r 
is  held  constant,  this  is  tantamount  to  including  diffusion  in  the 
positive  column.  However,  diffusion  was  not  included  as  a  loss  process 
in  the  electron  current  continuity  equation.  In  addition  to  this 
discrepancy,  the  anode  fall  voltage  appears  to  be  much  too  small. 
Estimating  from  their  Fig  1,  the  anode  fall  in  helium  appears  to  be 


less  than  5  volts  which  is  much  less  than  the  nominal  value  of  26  volts 
(ref  6:148).  Later  in  this  chapter  the  same  problem  was  encountered 
when  this  method  was  used.  Further  discussion  of  the  anode  fall  is 
reserved  for  Chapter  V. 

Recently,  Mitchell,  Kline,  and  Denes  (ref  67)  calculated  the 
electric  field,  particle  fluxes,  and  particle  densities  as  a  function 
of  position  in  the  cathode  and  anode  regions  of  a  high  pressure  glow 
discharge  containing  an  attaching  gas  mixture.  They  used  a  one  dimen¬ 
sional  steady  state  continuity  equation  model  including  Poisson's 
equation.  The  following  boundary  conditions  were  used: 

Cathode  Anode  Positive  Column 

7.  =  yX  n,  •  n.  •  n.  •  0  (III-29) 

J >0 

where  Jg/J_  and  E  were  adjustable  initial  parameters.  They  reported 
that  the  cathode  sheath  characteristics  are  insensitive  to  variations 
in  the  concentrations  of  the  attaching  species,  whether  diffusion  of 
electrons  were  included  or  not.  However,  the  anode  sheath  was  affected 
by  the  inclusion  of  electron  diffusion  as  well  as  the  attaching  concen¬ 
tration.  However,  as  will  be  shown  in  the  next  section,  the  boundary 
condition  for  the  positive  column  in  equation  1 1 1-29  may  not  lead  to  a 
unique  solution.  Thus  there  remains  some  question  on  the  convergence 
and  hence  self-consistency  of  their  solutions. 


kv 


New  Equi  1  ibriuir,  Analysis  of  the  Cathode  Region  Including  Negative  Ions 
This  section  expands  Neuringer's  analysis  of  the  cathode  fall 
region  to  include  negative  ions.  Since  the  ratio  of  positive  ions  to 
electron  mobility  is  held  constant,  the  transverse  diffusion  of  elec¬ 
trons  was  added  to  the  electron  continuity  equation.  A  one-dimensional 
steady  state  distribution  is  assumed  so  all  quantities  are  functions  of 
the  distance  x  from  the  cathode  surface  and  not  of  the  transverse 
geometry.  When  attachment  and  transverse  diffusion  processes  are 
included,  the  equations  governing  the  steady-state  discharge  are: 

Poisson's  equation 

JE. mzL  fJk  _  JL_ 

Jt  U  l  Y,  * 

the  electron  current  continuity  equation 

the  negative  ion  current  continuity  equation 

T-  n  T  kr  J-  J* 
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( 1 1 1-31 ) 


and  the  total  current  conservation  equation 


V*J~0  or  *L  ( II 1-33) 
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where  E  is  the  electric  field  strength,  J+,  Jg,  J_  and  J  are  the 
positive  ion,  electron,  negative  ion  and  total  current  densities 
respectively,  v+,  vg,  and  v_  are  the  positive  ion,  electron,  and 
negative  ion  drift  velocities,  is  the  permittivity  of  free  space,  a 
and  n  are  Townsend's  first  ionization  and  attachment  coefficients,  De 
is  the  electron  free  diffusion  coefficient  (cm2/sec),  *  is  the  electron 
mean  free  path,  is  the  ion-ion  recombination  rate  (cm3/sec),  and  x 
is  the  distance  from  the  cathode  to  the  anode.  These  equations  differ 
from  Davies  and  Evan's  equations  in  that  negative  ions  have  been  added 
to  Poisson's  equation,  attachment  and  electron  diffusion  have  been 
added  to  the  electron  continuity  equation,  the  negative  ion  current 
continuity  equation  is  now  required,  and  the  negative  ion  current 
density  has  been  added  to  the  total  current  conservation  equation. 

Normalizing  the  current  density  to  unity  implies: 
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Then  dividing  equation  I I 1-30  by  I I 1-31  and  using  I I 1-34  to  eliminate 
the  positive  ion  current  density  yields 
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( (III-35) 
Equations  1 1 1 -31 ,  III -32 ,  and  1 1 1 -35  are  the  set  which  will  be  analyzed 


instead  of  the  original  system  of  equations  1 1 1-30,  1 1 1-31 ,  and  1 1 1-32 . 
Because  a,  n>  v+,  vg  are  strong  functions  of  the  electric  field, 
this  system  of  coupled  equations  is  highly  nonlinear.  The  system 
becomes  solvable  if  the  following  assumptions  are  made: 

and  ^  a  J  ( 1 1 1-36 ) 

The  variation  in  r  with  the  electric  field  is  very  small  and  becomes 
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significant  only  when  J+  »rJe.  Thus  these  are  good  assumptions  over 
the  complete  range  of  E.  Making  use  of  these  assumptions,  separating 
variables  and  using  the  boundary  conditions  E/p  =  EQ/p  and  jg  =  jgo  at 
the  cathode  results  in  the  following  integral  equation 

H  A. 


Jh+Hh'M  -z  j.  if 


( II 1-37) 


The  Townsend  first  ionization  coefficient  that  was  used,  was  the 
more  general  form  given  by 


< . « -f«) 


( II 1-38) 


where  A  and  B  are  empirically  determined  constants  characteristic  of 
the  gas,  and  s*l,  for  molecular  gases  and  s=4  for  monatomic  gases. 
Similarly,  the  attachment  coefficient  can  be  expressed  as  an  exponen- 


%  -  c 


(II 1-39 ) 


or  as  a  series 


yp  -  «• 


(III -40 ) 


where  C,  D,  and  t  or  a^  are  empirically  determined  constants  depending 
on  the  attaching  gas.  These  are  good  assumptions  provided  the  elec¬ 
trons  are  in  equilibrium  with  the  field.  The  ion  drift  velocity  in  a 
high  field  region  can  be  described  by  Ward's  modified  formula. 


Substituting  for  a  ,  n,  and  v+,  equation  1 1 1-37  becomes: 

( Cpt**1'-  m -42) 


Changing  variables  by  letting  4>  =  B(p/E),  the  first  term  on  the  left 
side  of  the  equation  becomes  the  difference  of  the  two  incomplete  gamma 


functions. 


J Ap^yk.(%){pJ(%)  =  A p'L  bhJff 


*  t  Jf 


v(-h,bmi] 


(111-43) 


Where  Vf*-;x)a  (ref  9:260).  Similarly  changing 

variables  so  that  <J>  =  D(  /p)  ,  the  second  term  becomes: 

i^d(E/r>)  -  Tj0*J  0^ 


=7^[rao^)')-v(^Dff/,))] 


(111-44) 


The  third  term  integrates  easily  to  yield 

a  (III -45 ) 

4 

Integrating  the  right  hand  side  of  1 1 1-42 ,  yields 

j" 

'Trjtt’h-J-z  tJJJ- ~  t &(£)+(br)(i-JJ+zfj;  A  <  1 1 1-46) 


II 1-15 


Therefore  III -42  can  now  be  rewritten  as 


TZbtra,  mfj  4^‘fra  mb 

_  4,  (III -47 ) 

r^tra,^  -  jfe  m  -  ^uaui^L-i.)  *[$  jq 

The  gamma  functions  can  be  simplified  for  numerical  calculations  with 
the  help  of  the  following  recurrence  relation  (ref  15:457). 


*  *  v(*>*)  *  x 


q.  -Jr 
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( III -48 ) 


For  s  =  i,v(3/2s,  x)  reduces  to 

where  £,(*}*  ff*  ,  the  exponential  integral  (ref  9:228).  For 

X 

numerical  calculations  E^x)  can  easily  be  represented  by  polynomial 
approximation  (ref  9:231).  For  s  =  1,  r(3/2s,  x)  reduces  to 

+$-£*(■&  -  %*)  (III-50) 

00 
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where  ttfc(x)n^p)C  is  the  complementary  error  function  for 

which  there  are  tables. 

Similarly  for  t  =  1,  v (+3/2t,  x)  reduces  to 


( 1 1 1-51) 


These  simplifications  for  the  gamma  functions  above  were  used  in  the 
numerical  calculations  to  ease  computation.  Error  functions,  exponen¬ 
tials,  and  square  roots  were  more  accessible  in  computer  libraries  than 
incomplete  gamma  functions. 


These  terms  still  form  a  trancendental  equation  involving  j  and 


E.  Given  j  at  the  cathode,  various  mathematical  methods  can  be  used 


to  find  E  which  is  a  root  to  the  equation.  This  field  can  then  be  used 
in  solving  the  continuity  equations,  1 1 1 -31  and  1 1 1-32 ,  very  near  the 
cathode.  Iteratively,  one  can  spatially  step  from  the  cathode  to  the 


anode  by  alternatively  solving  for  E(jg)  and  j  .  This  method  elimi¬ 


nates  the  integration  problems  conventional  differential  equation 
solvers  have  in  solving  the  original  set  of  equations,  which  included 


Poisson's  equation  explicitly.  E(jg)  is  a  smooth  function  and  can  be 


determined  using  successive  numerical  approximations  to  solve  the 
transcendental  equation  1 1 1-47 .  See  Fig  III-l. 

An  expression  for  the  variation  of  the  electric  field  can  be 
obtained  by  first  taking  the  derivative  of  equation  III -47 .  The 
derivative  of  each  of  the  terms  in  I 11-47  is  given  below. 
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( II 1-52) 
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The  derivatives  involving  j  can  be  eliminated  by  substituting  in 


ELECTRIC  FIELD  AS  A  FUNCTION  OF  ELECTRON  CURRENT  DENSITY 


equation  1 1 1 -31  in  normalized  form  yielding 


( III -56 ) 


In  an  electropositive  gas,  j_  =  0  and  j  «1  very  near  the  cathode  so 
the  last  two  terms  in  1 1 1-56  are  negligible.  Integrating  equation 
1 1 1-56  under  these  conditions  yields 


( II 1-57) 


or  taking  the  first  two  terms  in  a  binomial  expansion  yields 


( II 1-58) 


This  is  the  same  result  Neuringer  obtained,  1 1 1-23 ,  using  many  more 
approximations  in  his  derivation.  He  retained  only  the  first  term  in 
equation  1 1 1-46  as  he  assumed  j  «  J  near  the  cathode.  This  resulted 
in  an  equation  similar  to  I I 1-57  which  contained  terms  proportional  to 
(E/p)*  which  he  conveniently  dropped.  The  binomial  expression  leading 
to  equation  1 1 1-58  converges  as  long  as  >  Tkfpc, . 

This  inequality  holds  in  the  cathode  fall  region  but  may  not  hold  in 
the  negative  glow  or  positive  column  region  where  EQ/p  becomes  very 
small  and  no  longer  a  linearly  decreasing  function. 

Returning  to  equation  1 1 1 -56 ,  the  dominant  terms  for  an  electro¬ 
positive  gas  are 

J_  T  T  U_  _  Jrut  (III-59) 
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which  is  consistent  with  Poisson's  equation  and  the  fact  that  the  field 
in  the  cathode  fall  region  is  determined  primarily  by  the  positive  ion 
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number  density. 

In  an  electronegative  gas,  the  situation  is  more  complex.  How¬ 
ever  very  near  the  cathode  even  in  an  electronegative  gas,  j_«j  and 
j  «1.  Thus  the  field  there  is  determined  primarily  by  the  positive 
ion  current  in  both  electronegative  and  positive  gases.  As  one  pro¬ 
ceeds  across  the  cathode  fall,  the  electron  and  negative  ion  current 
densities  increase  in  a  related  manner  given  by  1 1 1-32 .  This  has  the 
affect  of  increasing  the  slope  of  the  electric  field  as  can  be  seen 
from  the  following  analysis. 

It  was  found  from  the  results  of  the  next  section  that  the  nega¬ 
tive  ion  current  density  through  the  cathode  fall  region  could  be 
related  to  the  electron  current  density  by  an  exponential  function  of 
the  following  form, 

J.  =  «■  -  c  Ull-60) 

where  c  is  a  constant  determined  by  the  boundary  condition  that  j_  =  0 
at  the  cathode  and  m  is  an  empirically  determined  constant.  Results  of 
the  nonequilibrium  analysis  in  the  next  Chapter  also  confirm  that 
equation  1 1 1-60  is  a  good  approximation  throughout  the  cathode  fall 
region  except  for  a  very  small  region  a  few  electron  mean  free  paths  in 
length  in  front  of  the  cathode.  Using  this  approximation,  the  deriva¬ 
tive  of  the  integral  on  the  right  side  of  111-56  becomes 

z  (*  ‘  p  A  tmit)  (ii  I-6D 
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After  integrating  1 1 1 -56  and  using  this  expression,  one  can  solve  for 


I I 1-20 


the  field  as  a  function  of  distance 

*  t 

■4)^--  (%)  /ju«  -l  /e”4*  J  (111-62) 
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Taking  the  first  two  terms  in  a  binomial  expansion  as  before  yields 

T  /  * 

%  ~CfC(E7Jt  h(l*ZcHi*r) JUi-i  fsik  (III-63> 

In  the  electronegative  case,  the  field  is  shifted  higher  at  each 
corresponding  x  near  the  cathode  by  the  factor  (1  +  2c).  This  shift  is 
small  but  can  be  observed  in  the  results  presented  in  the  next  section 
as  well  as  the  results  of  the  nonequilibrium  analysis  in  Chapter  IV  for 
mixtures  of  Xe  and  HC1 .  Thus  this  analysis  predicts  that  the  attaching 
gas  slightly  lengthens  the  cathode  fall  and  slightly  increases  the 
cathode  fall  voltage.  As  the  electron  current  density  grows  with 
respect  to  distance  the  integral  terms  become  larger  and  have  the 
effect  of  increasing  the  field  in  the  positive  column  region  of  the 
discharge.  Thus  increasing  the  amount  of  an  attacher  also  increases 
the  field  in  the  positive  column  region  in  an  attempt  to  balance  the 
gain  and  loss  of  electrons. 

Returning  to  equation  1 1 1 -56  and  solving  for  the  minimum  slope  of 
the  electric  field  yields 

0=  011-64) 


The  minimum  field  in  an  electropositive  gas  occurs  when  j  =  1/1+r.  In 


the  electronegative  case,  the  minimum  field  occurs  when 


/  =  ii-u.) 
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( 1 11-65) 


Both  of  these  expressions  are  equivalent  to  charge  neutrality. 
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( 1 11-66) 


Description  of  Program  GLOW 

The  GLOW  program  was  developed  during  this  study  to  provide  a 
fast,  yet  complete  description  of  the  electric  field  in  a  gas  discharge 
from  cathode  to  anode,  similar  to  the  Davies  and  Evan's  model  (ref  32). 
This  type  of  model  was  chosen  because  it  can  give  a  complete  descrip¬ 
tion  of  the  electric  field  as  a  function  of  distance  from  cathode  to 
anode.  The  advantage  of  this  type  of  technique  is  that  it  minimizes 
the  number  of  boundary  conditions.  This  code  differs  from  most  other 
approaches  including  Davies  and  Evans  by  starting  at  the  cathode  and 
proceeding  towards  the  anode.  Most  other  approaches  start  at  the  anode 
or  positive  column  and  work  back  towards  the  cathode  due  to  stability 
problems.  This  method  is  sufficiently  stable  to  start  at  the  cathode 
where  the  boundary  conditions  on  the  negative  ion  current  are  well 
known.  The  theory  upon  which  the  model  is  based  has  already  been 
described  in  the  preceding  section.  This  section  will  describe  the 
boundary  conditions  and  the  iteration  techniques  in  detail. 

The  boundary  conditions  used  at  the  cathode  and  anode  are  ones  on 
which  there  is  fair  agreement.  The  boundary  conditions  used  at  the 


cathode  are 


J.  (O)  =>  jjy  J  (II 1-67 ) 

l(0}-0  (HI-68) 

where  yis  the  secondary  emission  coefficient  at  the  cathode.  The 
other  boundary  condition  requires  current  continuity  at  the  anode  such 
that 


It  +  T.mT  »»•  J*sO  (  hi -69 ) 

Other  techniques  such  as  starting  in  the  positive  column  using  a 
differential  equation  solver  and  starting  toward  the  cathode  usually 
require  one  or  more  boundary  conditions  in  addition  to  equation  1 1 1-67 
and  1 1 1-68  to  match  solutions  for  each  region  being  described. 

The  code  is  based  on  a  shooting  technique  which  varies  the  elec¬ 
tric  field  at  the  cathode  to  converge  upon  a  calculated  value  of  the 
electric  field  at  the  anode.  The  value  of  the  electric  field  at  which 
the  electron  current  density  plus  the  negative  ion  current  density 
equals  the  total  discharge  current  density  was  found  to  converge  very 
quickly  to  a  constant  value.  This  value  of  the  electric  field,  at 
which  the  anode  boundary  conditions  hold  was  the  field  that  was  the 
target  in  the  shooting  technique.  After  each  iteration  from  cathode  to 
anode  the  target  value  of  the  electric  field  at  the  anode  was  updated 
with  a  refined  value  for  which  equation  III -69  holds.  If  E(d)  was  too 


high  such  that 

I-J.  >J  (m-70) 

then  E (0 )  was  decreased.  If  E(d)  was  too  low,  or  the  calculation  did 
not  reach  the  anode,  then  E ( 0 )  was  increased.  Fig  1 1 1-2  shows  the 
electric  field  as  the  program  converged.  For  the  second  and  twelfth 
iterations  in  Fig  1 1 1-2 ,  no  solution  existed  for  the  transcendental 
equation  past  the  points  plotted.  This  implies  that  this  technique 
converges  not  only  to  a  solution  which  satisfies  the  boundary  equa¬ 
tions,  but  is  also  the  lowest  voltage  for  the  discharge  for  the  initial 
conditions. 

Any  numerical  technique  for  finding  the  roots  of  an  equation  could 
be  used  in  solving  the  transcendental  equation,  such  as  the  method  of 
halving  the  interval.  However,  a  modified  linear  interpolation  method 
(ref  45:10-11)  was  used  since  it  is  supposed  to  be  sligthly  faster  than 
the  halving  the  interval  method.  The  only  requirement  was  that  the 
first  two  estimates  for  E(0)  should  bound  the  solution,  i.e.  one 
estimate  be  too  large  and  the  other  too  small.  Rapid  convergence  to 
4-6  digit  accuracy  was  achieved  within  20  iterations  from  electric 
fields  differing  as  much  as  10-20v/cm  at  the  cathode.  Numerical 
calculations  were  terminated  when  the  target,  E(d),  was  within  10%  of 
the  calculated  value. 

Within  each  of  these  iterations  from  cathode  to  anode,  a  differen¬ 
tial  equation  solver,  and  this  modified  linear  interpolation  method  was 
used  to  solve  the  system  of  current  continuity  equations.  III -31  and 
III-37,  and  find  the  root  to  the  transcendental  equation,  III-47, 
respectively.  The  variable  order  Adams  predictor-corrector  method 


I TI-24 


which  was  part  of  a  differential  equation  solver  package  called  DGEAR 
(ref  36)  for  solving  stiff  differential  equations  was  used  to  solve  for 
j  given  an  electric  field.  The  overall  iteration  scheme  from  cathode 
to  anode  went  as  follows: 

a.  Given  E(0),  DGEAR  was  called  to  solve  the  set  of  equations 
1 1 1 -31  and  1 1 1-32  for  jg(x)  a  short  distance  from  the  cathode. 

b.  Given  this  value  of  je(x),  the  modified  linear  interpolation 
method  was  used  to  find  the  root  or  equilibrium  value  for  E(x)  from  the 
transcendental  equation,  III  -47. 

c.  Given  this  value  of  E(x),  DGEAR  was  called  again  to  solve  for 
a  new  j g ( x )  at  the  next  small  increment  in  space.  The  spatial  step 
size  h  was  usually  picked  such  that  (E(0)  -  E(h))/E(0)<  1%.  This 
cycle  of  calculating  jg(x)  and  E(x)  was  repeated  as  the  program  stepped 
from  cathode  to  anode. 

The  program  can  be  run  without  the  help  of  the  differential 
equation  solver,  using  the  field  at  the  last  step  to  calculate  the 
field  from  equation  1 1 1-47 .  However,  this  method  was  very  sensitive  to 
the  spatial  step  size  used.  The  addition  of  a  differential  equation 
solver  to  predict  the  electron  current  density  at  the  next  step  where 
the  transcendental  equation  III -47  was  to  be  solved  next,  enabled  much 
larger  steps  to  be  taken  with  a  savings  in  computer  time. 

The  GLOW  program,  incorporating  both  a  differential  equation 
solver  and  the  modified  linear  interpolation  method  was  usually  run  on 
a  Cyber  750  computer  system.  Each  iteration  took  on  an  average  less 
than  4  sec  of  execution  time  for  250-300  spatial  steps.  It  required 
less  than  52k  of  memory  to  load  and  run.  Parameters  used  in  the 
numerical  calculations  are  described  in  Appendix  C. 


To  establish  a  baseline  and  verify  the  accuracy  of  the  technique* 
Davies  and  Evans'  boundary  conditions  and  initial  values  were  used  and 
the  resultant  field  plotted  in  Fig  III-3.  They  set  the  electron 
current  density  equal  to  the  total  discharge  current  and  then  varied 
the  electric  field  at  the  anode  until  a  fixed  value  of  the  field  was 
obtained  at  the  cathode.  By  starting  at  the  anode,  r  was  not  required 
as  a  boundary  condition.  The  GLOW  code  was  first  run  from  anode  to 
cathode  to  determine  Y  at  the  cathode.  The  code  was  then  run  from 
cathode  to  anode  using  this  value  of  Y=  .293.  In  this  instance,  Y  was 
fixed  and  the  field  at  the  cathode  was  varied  until  the  electron 
current  density  equaled  the  total  discharge  current  density  at  the 
anode.  As  can  be  seen  form  Fig  III-3,  very  good  agreement  of  the  GLOW 
code  (cathode  to  anode)  is  obtained  with  Davies  and  Evan's  result 
(anode  to  cathode). 

Results  of  Numerical  Calculations 

This  section  presents  the  results  of  calculations  using  the  GLOW 

code  in  various  mixtures  of  rare  gases  with  HC1.  All  data  presented 

for  He  and  Ar  was  obtained  for  a  gas  pressure  of  1  torr.  The  data  for 

Xe  was  calculated  for  a  gas  pressure  of  100  torr.  However,  it  was 

found  that  for  all  three  rare-gas  mixtures,  the  results  at  different 

pressures  obeyed  the  following  pressure  scaling  relationships. 

V  =  constant 
c 

J/p2  =  constant 
pd  =  constant 


This  finding  is  in  accordance  with  the  similarity  rules  or  scaling  laws 


FIG  III-3  COMPARISON  OF  THE  ELECTRIC  FIELD  WITH  OTHER  DATA 


for  gas  discharges  discussed  in  Chapter  I.  Data  was  taken  at  1,  10, 
and  100  torr  for  two  gas  mixtures  and  was  found  to  deviate  less  than 
.01%  from  these  scaling  relationships. 

A  current-voltage  relationship,  similar  to  equation  1 1 1-8  can  be 
derived  for  the  high  field  case  where  the  positive  ion  velocity  is 
proportional  to  the  square  root  of  the  field.  Following  von  Engel  and 
Steenbeck's  derivation,  but  using  equation  III -41  for  the  ion  drift 
velocity  results  in  the  following  current-voltage  relationship. 


Thus,  this  new  scaling  function 
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is  a  constant  for  a  given  y  and  k+  in  the  high  field  case.  This 
scaling  relationship  is  useful  in  comparing  data  at  different  voltages 
and  currents. 

Results  will  be  presented  in  the  sequence  He  mixtures,  Ar  mix¬ 
tures,  and  then  Xe  mixtures.  Results  will  be  discussed  in  this  section 
for  the  cathode  fall  and  positive  column  regions,  but  not  for  the  anode 
fall  region.  At  present  the  GLOW  code  does  not  adequately  describe  the 
anode  region.  The  anode  fall  region  is  discussed  in  Chapter  V. 


Helium  Mixtures 


This  section  presents  the  results  obtained  for  He/HCl  mixtures. 


given  p  =  1  torr,  d  =  3cm,  J  =  1.6E-5  amp/cm2,  and  y  =  .2.  A  summary 
of  numerical  values  for  the  following  figures  is  given  in  Table  1 1 1 - 1 


for  EQ/p,  Emi- p/p ,  pdc>  J  /p2,  and  Vc.  Eq  is  the  electric  field  at  the 
cathode,  Emin  is  the  minimum  electric  field  which  in  these  calculations 
occurs  in  the  positive  column,  p  is  the  total  gas  pressure,  dc  is  the 
length  of  the  cathode  fall  region,  J  is  the  total  current  through  the 
discharge,  and  Vc  is  the  cathode  fall  voltage.  Consistent  with  the 
somewhat  arbitrary  but  common  scheme  for  calculating  the  cathode  fall 
distance,  the  linear  cathode  field  region  was  extrapolated  to  zero  to 
determine  the  cathode  fall  length  dc<  The  cathode  fall  voltage  was 
obtained  from  the  voltage  between  the  cathode  and  this  point.  (Other 
definitions  can  be  used,  but  for  the  comparison  with  other  published 
models  this  definition  was  applied.  Experimental  comparisons  would 
have  slightly  lower  values).  The  last  column  in  Table  III-l  is  the 
scaling  function  (SF)  given  in  equation  1 1 1-72.  This  scaling  relation¬ 
ship  helps  to  put  in  perspective  the  various  data  on  cathode  fall 
parameters.  The  value  of  SF(y.  k+)  for  100%  He  agrees  only  within 
62.5%  of  Ward's  data  (ref  92).  This  discrepancy  is  largely  due  to  the 
different  approximation  for  the  positive  ion  drift  velocity  which  he 
used: 
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Data  was  taken  from  a  summary  presented  in  ref  6. 
Data  was  taken  from  ref  92. 

Data  was  taken  from  ref  32. 

Data  was  calculated  in  this  study. 


(where  for  He 


=  8.E3  cm2-torr/volt-sec. 


V  + 

k+  =  4.1E4  cm^-torr*/vo1t*-sec,  C  =  8.E-3  torr-cm/volt, 

D  =  27.44  volt*/cm^-torr*  and  Wj  =  25v/cm-torr). 

At  200v/cm-torr,  Ward's  approximation  yields  a  positive  ion  velocity 
that  is  only  1%  larger,  whereas  at  30  v/cm-torr  it  is  17%  larger. 

Since  the  slope  of  the  electric  field  is  inversely  proportional  to  the 
positive  ion  drift  velocity,  the  cathode  fall  length  and  voltage  are 
also  sensitive  functions  of  the  positive  ion  drift  velocity  and  largely 
accounts  for  the  above  discrepancy.  Very  good  agreement  is  reached 
with  Davies  and  Evans  data  (ref  32)  for  y  =  .23.  Care  should  be  taken 
in  comparing  the  values  of  SF(Y,k+)  since  this  'constant'  is  propor¬ 
tional  to  (1  +  Y). 

Fig  III -4  shows  a  3%  decrease  in  the  electric  field  through  most 
of  the  cathode  fall  region  (0.-1.26  cm)  as  the  amount  of  HC1  is  in¬ 
creased  from  0  to  5%.  This  is  contrary  to  the  results  that  might  be 
anticipated  from  equation  1 1 1-63.  However,  HC1  has  an  ionization 
potential  approximately  half  of  that  for  He  (12.74eV  as  compared  to 
24.59  eV),  the  HC1  contributes  more  electrons  through  ionization  than 
the  discharge  loses  to  attachment  to  HC1  resulting  in  a  slightly 
reduced  cathode  fall  voltage  of  up  to  5%  for  the  95/5  He/HCl  mixture. 
There  is  also  a  trend  for  the  cathode  fall  width  and  voltage  to  de¬ 
crease:  1%  and  1.4%  respectively  for  99/1  mixture  and  5%  and  6.3%  for 
the  95/5  mixture.  In  the  positive  column  region  (around  2.5  cm)  the 
electric  field  is  5%  higher  for  99/1  mixture,  and  25%  higher  for  a  95/5 
mixture  of  He/HCl  as  a  result  of  the  increase  in  the  number  of  attach¬ 
ing  molecules.  The  fields  in  the  positive  column  region  are  consistent 
with  the  value  of  the  field  required  to  have  the  electron  gain  and  loss 
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processes  equal.  In  the  nearly  uniform  field  region  of  the  positive 
column  any  additional  loss  of  electrons  must  be  balanced  by  an  increase 
in  ionization  in  order  to  maintain  constant  current. 

Fig  1 1 1 -5  shows  the  same  calculation  in  95/5  He/HCl  with  and 
without  the  contribution  of  ionization  from  HC1 .  As  expected,  the 
electric  field  is  increased  in  the  cathode  fall  region  when  ionization 
of  HC1  is  not  included.  This  additional  ionization  in  the  cathode  fall 
region  as  the  percentage  of  HC1  is  increased  results  in  a  small  rela¬ 
tive  increase  in  the  electron  current  density  in  Fig  1 1 1 -6  in  the 
cathode  fall  region.  The  current  densities  plotted  in  Fig  1 1 1 -6  and 
1 1 1-7  have  been  normalized  as  in  equation  1 1 1-34 .  The  electron  current 
density  decreases  in  the  positive  column  region  as  the  percentage  of 
HC1  is  increased  indicating  electrons  are  being  lost  to  the  formation 
of  negative  ions.  As  expected  the  negative  ion  current  density  shown 
in  Fig  1 1 1-7  becomes  a  larger  fraction  of  the  discharge  current  density 
through  both  the  cathode  fall  region  and  the  positive  column  region. 
However  it  still  contributes  less  than  1%  to  the  total  discharge 
current. 

Argon 

The  results  for  Ar/HCl  mixtures  will  be  presented  next  where 
p=ltorr,  d  =  1.5cm,  J  =  l.E-5  amp/cm2  and  y  =  .04167.  A  summary  of 
numerical  values  for  the  following  figures  is  given  in  Table  1 1 1 -2  as 
was  done  for  He/HCl  mixtures. 

In  100%  Ar  the  value  of  SF(y,  k+)  agrees  very  well  with  Ward's 
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FIG  III-6  NORMALIZED  ELECTRON  CURRENT  DENSITY  IN  He/HCl  MIXTURES 


value,  within  1%.  Better  agreement  is  achieved  in  Ar  than  in  He 
because  the  slope  of  the  electric  field  is  steeper  in  Ar  resulting  in 
less  error  in  the  cathode  fall  length  and  voltage.  Approximately  the 
same  discrepancy  exists  in  the  positive  ion  drift  velocity  in  Ar  as  in 
He.  At  400v/cm-torr,  Ward's  approximation  in  equation  1 1 1 -73 
(where  for  Ar  p  +  =  1.  E3cm2-torr/vol  t-sec , 

k+  =  8.25E3  cm^-torrVvolt^-sec,  C  =  2.22E-3  torr-cm/vol t, 

D  =  86.52  vol t^/cm^-torr* ,  and  =  60v/cm  torr) 

yields  a  positive  ion  velocity  that  again  is  about  1%  larger,  whereas 
at  60v/cm-torr  it  is  18%  larger. 

Contrary  to  He  mixtures  there  is  a  negligible  effect  on  the 
electric  field  in  the  cathode  fall  region  (0.-.66cm)  in  Fig  III -8  when 
up  to  5%  HC1  is  added  to  an  Ar  discharge.  In  this  case,  since  HC1  has 
a  slightly  lower  ionization  potential  than  that  for  Ar  (12.7eV  as 
compared  to  15.7eV),  the  small  additional  ionization  due  to  the  small 
percentage  of  HC1  is  insufficient  to  be  seen  as  any  effect  on  the 
electric  field  in  the  cathode  fall  region.  Fig  I I I -9  shows  the  same 
calculation  in  95/5  Ar/HCl  with  and  without  the  contribution  of  ioni¬ 
zation  from  HC1.  As  expected,  the  electric  field  is  increased,  but  not 
as  much  as  for  the  He/HCl  case  since  the  thresholds  for  ionization  are 
much  closer.  Similarly  no  change  is  seen  in  Fig  1 1 1 - 1 0  in  the  electron 
current  density  in  the  cathode  fall  region.  The  negative  ion  current 
density  in  Fig  1 1 1 -11  never  becomes  a  significant  fraction  of  the  total 
discharge  current  density.  The  negative  ion  density  increases  much 
more  slowly  than  for  the  He  case,  since  the  fields  in  the  cathode  fall 
and  positive  column  are  higher  in  Ar  than  in  He. 
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Data  was  taken  from  a  summary  presented  in  ref  6. 
Data  was  taken  from  ref  85. 

Data  was  calculated  in  this  study. 
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fig  III-9  EFFECT  OF  IONIZATION  OF  HC1  IN  Ar/HCl  MIXTURES 
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NORMALIZED  ELECTRON  CURRENT  DENSITY  IN  Ar/HCl  MIXTURES 
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In  the  positive  column  (around  1.2cm)  the  electric  field  is 
increased  approximately  3%  for  a  99/1  Ar/HCl  mixture  (not  plotted)  and 
almost  12%  increase  for  a  95/5  mixture.  Again  this  is  the  result  of  an 
increase  in  the  number  of  attaching  molecules  in  the  positive  column 
region. 


Xenon 

The  results  for  Xe/HCl  mixtures  with  p  =  100  torr,  d  =  ,01cm, 

J  =  .1  amp/cm2,  and  y  =  .004  will  be  presented  next.  This  value  of  y 
was  chosen  so  results  could  be  compared  with  Ward's  results  (ref  85). 

A  summary  of  numerical  values  for  the  following  figures  is  given  in 
Table  1 1 1-3  as  was  done  for  the  previous  mixtures. 

In  100%  Xe,  there  is  apparent  agreement  (6%)  between  the  scaling 
factors  calculated  from  Wards  data  and  the  data  from  the  GLOW  code. 
Approximately  the  same  discrepancy  exists  in  the  positive  ion  drift 
velocity  as  in  the  previous  rare  gases.  At  690  v/cm-torr.  Ward's 
approximation  using  equation  I I 1-73 

(where  in  Xe  y  +  -  4.E2  cm  torr/volt-sec, 

k+  =  4.E3  cm^-TorrVvolt^-sec,  C  =  2.25E-3  torr-cm/volt, 

D  =  225  volt^/cm^-torr^,  and  W^  =  100  v/cm-torr) 

yields  a  positive  ion  velocity  that  is  about  1%  larger,  whereas  at  100 
v/cm-torr,  it  is  23%  larger. 

In  Xe  mixtures  there  is  a  very  slight  trend  for  the  electric  field 
to  increase  in  the  cathode  fall  region  (0.-.6cm)  as  the  percentage  of 
HC1  is  increased.  This  is  observable  in  Fig  1 1 1 - 12 .  Xe/HCl  mixtures 
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Data  was  taken  from  ref  85. 

Data  was  calculated  in  this  study. 
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FIG  III-12  ELECTRIC  FIELD  IN  Xe/HCl  MIXTURES 


are  unique  in  that  both  constituents  have  almost  the  same  ionization 
potential  ( 12 . 74eV  for  HC1  and  12. 13eV  for  Xe).  Thus  only  effects  of 
attachment  are  observed  and  the  results  agree  with  the  analytical 
description  given  earlier.  The  effect  of  attachment  is  to  raise  the 
field  slightly  throughout  the  cathode  fall  region. 

Even  though  the  field  increases  as  HC1  is  added  to  Xe,  the  elec¬ 
tron  current  density  decreases  slightly  throughout  most  of  the  cathode 
fall  region  in  Fig  III-13  due  to  the  pure  loss  of  electrons  to  attach¬ 
ment.  As  in  Ar  mixtures,  the  negative  ion  current  density  in  Fig 
1 1 1 -1 4  grows  very  slowly  near  the  cathode  but  begins  to  grow  rapidly  as 
the  electric  field  levels  off  into  the  positive  column.  It  is  still 
insignificant  in  comparison  to  the  contribution  from  the  electron  and 
positive  ion  current  densities  to  the  total  discharge  current  density. 

In  the  positive  column  (around  .009  cm)  the  electric  field  in¬ 
creases  approximately  2%  for  99/1  mixture,  4%  for  95/5  mixture,  and  9% 
for  90/10  mixture.  Again  this  is  the  result  of  an  increase  in  the  num¬ 
ber  of  attaching  species  in  the  positive  column  and  the  imposed  boun¬ 
dary  condition  of  constant  current. 

In  summary,  the  above  calculations,  which  implicitly  assumed  that 
the  electrons  were  in  equilibrium  with  the  local  field,  have  shown 
that: 

a.  The  formation  of  small  amounts  of  negative  ions  in  the 
cathode  fall  region  is  not  the  reason  that  the  cathode  a  length  con¬ 
tracts  in  the  axial  direction. 

b.  The  observed  contraction  of  the  cathode  full  length  is  a  re¬ 
sult  of  the  increased  ionization  from  the  more  efficient  ionization  of 
the  attaching  gas  in  some  gas  mixtures.  The  contraction  will  only 
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appear  in  gas  mixtures  in  which  the  ionization  rate  of  the  attaching 
gas  is  larger  than  that  of  the  background  gas. 

c.  When  the  ionization  rates  and  thresholds  of  both  gases  in  the 
mixture  are  similar,  then  the  formation  of  negative  ions  tends  to  in¬ 
crease  the  voltage  drop  across  the  cathode  fall  region  and  to  increase 
the  field  at  the  cathode. 

d.  The  field  in  the  positive  column  region  increases  as  the 
percentage  of  attacher  is  increased.  This  is  required  in  order  to 
balance  the  electron  gain  and  loss  processes  and  to  maintain  a  constant 
current. 
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Chapter  IV.  Electron  Kinetics  in  the 
Cathode  Fall  Region:  Nonequi librium  Analysis 


Review  of  Nonequilibrium  Analyses 

The  analyses  in  Chapter  III  were  all  based  on  Poisson's  equation 
and  current  continuity  or  flux  equations.  These  analyses  presuppose 
that  the  electron  energy  distribution  function  is  in  equilibrium  with 
the  local  electric  field.  As  mentioned  in  Chapter  I,  there  is  analyti¬ 
cal  evidence  that  the  electrons  are  not  in  equilibrium  with  the  electric 
field  in  the  cathode  fall  region.  Since  the  electron  current  or  flux 
is  an  average  over  the  product  of  the  electron  number  density  times  the 
electron  velocity,  <nv>,  it  is  not  possible  to  derive  n  or  v  from  this 
average  unless  the  nonequilibrium  electron  energy  distribution  function 
is  known.  There  are  two  methods  of  calculating  the  electron  energy  dis¬ 
tribution  function: 

1)  by  using  a  Monte  Carlo  simulation  of  the  electrons  as  they  drift 
from  cathode  to  anode,  or 

2)  by  calculating  the  distribution  function  from  the  Boltzmann 
equation. 

This  chapter  will  critically  review  and,  where  appropriate,  point 
out  defects  in  the  few  papers  that  have  investigated  the  nonequilibrium 
electron  kinetics  using  either  of  these  two  methods.  The  new  results  of 
this  study,  which  used  a  self-consistent  Boltzmann  method  and  included 
dissociative  attachment,  will  be  presented  at  the  end  of  the  chapter. 

First,  a  comparison  will  be  made  between  the  Monte  Carlo  method  and 
the  Boltzmann  method.  Although  both  methods  should  yield  similar 
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results,  there  are  important  computer  and  scientific  differences.  The 
Monte  Carlo  method  is  much  more  costly  in  terms  of  computer  time.  In 
order  to  obtain  good  statistics,  a  large  number  of  collisions  are 
required,  typically  10  or  greater.  The  advantage  of  the  Monte  Carlo 
method  lies  in  its  ease  of  implementation  as  a  computer  code.  However, 
the  Boltzmann  method  normally  enables  a  greater  insight  into  the  physics 
of  the  situation.  The  sensitivity  to  changes  in  inputs  is  more  easily 
investigated  primarily  due  to  the  shorter  run  time.  This  enables 
approximate  analytical  solutions  to  be  found  which  illustrate  the  depen¬ 
dence  on  the  relevant  parameters.  In  uniform  electric  fields,  the 
agreement  between  Monte  Carlo  results  and  Boltzmann  analyses  have  been 
excellent  (ref  62,  72,  77). 

The  theory  and  first  calculation  of  the  electron  distribution 
function  as  a  function  of  distance  through  the  cathode  fall  region  was 
reported  by  Allis  and  co-workers  (ref  13,  14)  in  1975  and  1977.  Their 
approach  involved  transforming  Boltzmann's  equation  into  energy  space 
where  the  new  variables  are  the  electron  potential  energy,  <j>=e/  E(x)dx 
and  the  electron  kinetic  energy  £  =  imv2.  In  addition,  they  treated 
the  cathode  fall  region  as  a  separate  entity  allowing  for  no  flux  of 
positive  ions  from  the  negative  glow  and  forcing  the  electric  field  to 
zero  at  the  anode.  Their  calculations  were  based  on  an  'ideal'  gas 
defined  as  having  a  constant  momentum  transfer  cross  section  and  con¬ 
stant  excitation  and  ionization  cross  sections  with  thresholds  of  8eV 
and  16eV  respectively.  In  addition  their  calculations  are  not  con¬ 
sidered  self-consistent,  meaning  the  calculations  of  the  electron  num¬ 
ber  density  and  electric  field  were  not  iterated  on  successively  to 
on verge  to  a  sel f-consi stent  solution.  The  nonequilibrium  phenomena 
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which  their  calculations  demonstrated  included: 

a.  A  dip  in  the  electron  number  density  adjacent  to  the  cathode. 
This  res.  Its  from  the  electrons  being  accelerated  away  from  the  cathode 
and  not  yet  having  sufficient  energy  for  ionization. 

b.  A  peaked  ionization  coefficient  which  is  small  near  the  cath¬ 
ode  where  the  field  is  highest  and  peaking  between  one  half  and  two 
thirds  the  total  cathode  fall  length. 

This  technique  of  transforming  Boltzmann's  equation  into  energy  space 
has  become  the  basis  of  several  investigations  discussed  later. 

Almost  simultaneously  with  Allis,  Tran  Ngoc,  et  al .  (ref  86)  used 
a  one  dimensional  Monte  Carlo  method  to  investigate  the  nonequilibrium 
electron  kinetics  of  the  cathode  fall.  They  began  by  assuming  a  lin¬ 
early  decreasing  electric  field.  Although  no  attempt  was  made  to  cor¬ 
rect  the  field  through  Poisson's  equation,  their  results  did  give  a  good 
physical  picture  of  the  nonequilibrium  process  occurring  in  a  high  field 
gradient.  They  found  the  electron  distribution  function  evolved  from  a 
sharply  peaked  distribution  near  the  cathode  to  a  multiply  peaked 
distribution  entering  the  negative  glow.  This  was  the  first  theoretical 
prediction  of  the  three  groups  of  electrons  observed  experimentally  by 
Pringle  and  Farvis  (ref  73)  in  the  negative  glow  (high  energy  electrons 
which  have  undergone  few  collisions,  intermediate  energy  electrons  which 
form  the  bulk  of  the  distribution,  and  low  energy  electrons  which  were 
primarily  a  result  of  ionization  processes).  Tran  Ngoc,  et  al.  conclu¬ 
ded  that  the  use  of  Townsend's  ionization  coefficient  given  by  equation 
1 1 1-38  was  inappropriate  in  the  cathode  fall  region  since  the  electrons 
never  reached  equilibrium  with  the  field.  In  addition  they  determined 
that  the  major  features  of  the  cathode  fall  are  actually  determined  by  a 


small  number  of  inelastic  collisions  and  that  elastic  collisions  play 
only  a  subsidiary  role. 

The  most  definitive  work  using  a  Boltzmann  analysis  was  accom¬ 
plished  by  Long  (ref  60).  It  was  based  upon  an  approach  outlined  by 
Allis  (ref  13)  of  a  direct  solution  of  the  Boltzmann  equation  without 
the  conventional  expansion  methods.  Long  simplified  Allis'  result  by 
simplifying  the  angular  dependence  of  the  distribution  function  in  the 
cathode  fall  region.  His  method  of  solution  was  first  to  estimate  the 
distribution  function  and,  to  use  this  in  the  collision  terms  to  inte¬ 
grate  the  distribution  function  along  the  various  energy  characteris¬ 
tics.  The  distribution  function  was  iterated  until  there  was  less  than 
.21  change  in  the  distribution  function  at  each  energy  bin  and  then  the 
electron  number  densities  and  electron  current  densities  were  calculated 
from  the  distribution  function  as  a  function  of  distance.  He  attempted 
to  use  these  values  to  solve  Poisson's  equation  for  the  electric  field 
before  repeating  the  process  to  converge  on  a  self-consistent  solution. 
Long's  results  and  conclusions  for  mercury  and  argon  agreed  very  well 
with  those  of  Tran  Ngoc  for  helium.  However,  in  his  integration  of 
Poisson's  equation  for  the  electric  field,  he  failed  to  include  any 
change  in  the  ion  drift  velocity  due  to  changes  in  the  electric  field 
that  occurred  from  assuming  current  continuity.  This  limited  the  con¬ 
vergence  of  his  code  to  unique  cases  where  the  anode  was  placed  inside 
the  cathode  sheath.  The  resolution  of  this  problem  and  further  improve¬ 
ments  of  Long's  computer  code  are  described  in  the  following  section. 

His  program  required  202K  of  memory  to  load  and  run.  Run  times  were 
typically  5-10sec  per  field  iteration  on  a  CDC  6600  computer. 

More  recently  Segur,  et  al .  (ref  79)  coupled  a  Boltzmann  solution 


for  the  electrons  with  an  analytical  approximation  to  the  Boltzmann 
equation  for  positive  ions  in  the  cathode  fall  region  for  various 
hel ium-mercury  mixtures.  They  assumed  the  electron  collisions  were 
isotropic  and  only  considered  charge  transfer  collisions  with  a  constant 
cross  section  for  the  positive  ions.  They  did  not  have  to  use  conserva¬ 
tion  of  current  to  calculate  the  positive  ion  density  as  a  function  of 
position  as  Long  had  done.  Segur,  et  al.  initially  assumed  a  linearly 
decreasing  electric  field  vanishing  at  the  anode.  Within  3-4  iterations 
of  calculating  the  electron  and  ion  densities,  followed  by  solving 
Poisson's  equation  for  the  field,  their  numerical  technique  achieved  a 
stable  convergence.  Thus  their  technique  was  self-consistent.  As  the 
percentage  of  mercury  was  increased  from  0%  to  9.1%,  the  peak  of  the 
ionization  coefficient  increased  and  the  threshold  for  ionization 
shifted  towards  the  cathode.  This  is  due  to  the  fact  that  mercury  atoms 
have  a  larger  ionization  cross  section  as  well  as  a  lower  ionization 
threshold  than  helium  atoms.  They  found  the  ionization  in  the  cathode 
fall  region  took  place  in  two  stages:  mercury  is  ionized  first  near  the 
cathode  and  helium  is  ionized  slightly  further  away  from  the  cathode. 
Eventually  both  processes  exist  together  during  the  remaining  part  of 
the  discharge.  Since  the  mobility  of  the  mercury  atoms  is  much  lower 
than  helium  atoms,  they  also  found  a  smaller  ionic  current  near  the 
cathode  and  a  larger  ionic  current  near  the  anode  when  mercury  was 
added. 

At  the  same  time,  Boeuf,  et  al .  (ref  21)  investigated  the  cathode 
fall  in  a  helium  discharge  using  a  self-consistent  Monte  Carlo  tech¬ 
nique.  They  also  investigated  some  features  in  an  oxygen  discharge. 

They  performed  one  Monte  Carlo  calculation  for  positive  ions  and  found 


that  the  positive  ion  drift  velocity  reached  equilibrium  with  the  field 

_2 

in  less  than  10  cm.  Thereafter  they  assumed  that  the  ions  drifted 
towards  the  cathode  with  their  equilibrium  drift  velocity.  In  their 
self-consistent  iterations  between  the  number  density  and  the  field, 
their  calculations  normally  reached  convergence  after  5  iterations. 

They  concluded  from  their  ionization  coefficient  and  electron  drift 
velocity  data  that  the  nonequilibrium  region  in  a  helium  discharge 
extends  just  over  half  (54%)  of  the  cathode  fall  length  which  agreed 
with  the  previous  investigations. 

A  three  dimensional  Monte  Carlo  calculation  was  accomplished  very 
recently  by  Boeuf,  et  al.  (ref  20).  They  made  extensive  use  of  the  null 
collision  technique  using  a  fictitious  total  cross  section  to  maintain 
a  pseudo-constant  collision  frequency  so  that  numerical  integrations 
could  be  avoided  in  determining  the  time  interval  between  two  collisions. 
This  allowed  them  to  calculate  a  three-dimensional  electron  distribution 
function  more  efficiently  taking  angular  scattering  into  account. 

Radial  fields  were  not  included  and  no  attempt  was  made  to  iterate 
successively  on  the  electron  number  density  and  electric  field.  Thus 
their  calculations  were  not  self-consistent.  They  examined  two  cases: 

1)  a  normal  glow  discharge  in  helium  at  one  torr  with  an  absorbing 
anode  placed  in  the  negative  glow,  and  2)  an  abnormal  discharge  in 
helium  at  one  torr.  They  found  that  there  is  little  error  if  isotropic 
scattering  is  assumed  and  the  momentum  transfer  cross  section  is  used  as 
the  elastic  cross  section  instead  of  assuming  angular  scattering  and 
using  the  elastic  differential  cross  section.  Boeuf  et  al.  showed  that 
the  original  assumption  of  forward  scattering  made  by  Tran  Ngoc  (ref  86) 
led  to  too  high  of  electron  energies  in  the  latter  two  thirds  of  the 


cathode  fall.  As  a  result  of  all  the  electrons  being  directed  in  the 
same  direction  and  at  too  high  of  energy,  the  Townsend  ionization  coef¬ 
ficient  for  the  forward  scattering  case  was  as  much  as  25%  too  low 
through  most  of  the  latter  two  thirds  of  the  cathode  fall  region.  They 
calculated  the  Townsend  ionization  coefficient  by  integrating  the  ioni¬ 
zation  cross  section  and  electron  distribution  function.  Its  maximum 
occurred  somewhat  beyond  the  maximum  of  the  mean  energy  distribution. 
They  concluded  that  as  a  result  of  the  bimodal  character  of  the  elec¬ 
tron  distribution  function, the  mean  energy  and  ionization  coefficient 
were  not  related  to  the  same  part  of  the  distribution  function. 

Finally,  their  angular  distribution  functions  show  that  the  high  energy 
segment  of  the  distribution  function  is  largely  forward  directed,  while 
that  of  the  low  energy  part  is  much  more  isotropic.  As  a  function  of 
distance  through  the  cathode  fall  region,  this  trend  became  much  more 
evident  close  to  the  negative  glow.  This  concludes  the  critical  review 
of  nonequilibrium  techniques. 


Description  and  Modifications  to  Proqram  SHEATH 


The  SHEATH  program  was  developed  by  W.  H.  Long  to  solve  the 
Boltzmann  equation  for  electrons  in  non-uniform  fields  without  requir 
ing  the  electron  distribution  function  to  be  in  equilibrium  with  the 
local  field.  The  collisional  Boltzmann  equation  describes  the  evolu¬ 
tion  of  a  swarm  of  electrons  in  both  space  and  time,  and  is  given  by 


«r 


where  f  =  f(r,v,t)  is  the  electron  distribution  function  which 
represents  a  complete  statistical  description  of  the  electron  in  a 
weakly  ionized  gas.  The  variable  t  is  time,  v  is  the  electron  velocity 
vector,  r  is  the  electron  spatial  vector,  a  =  -e/m  (E  +  v  +  B)  is  the 
acceleration  of  the  electrons  due  to  electric  and  magnetic  fields. 


wLi  represents  the  change  in  the  distribution  function  due  to 
elastic  and  inelastic  electron  collisions. 

Considering  only  spatial  variations  in  the  direction  of  the  elec¬ 


tric  field,  then  IV-1  reduces  to 


*  if +m‘£Wif=ifl 


(IV-2) 


Dividing  both  sides  by  vx  and  eE(x),  equation  IV-2  can  be  transformed 


into  energy  space. 


ii.  ,  )f-  1  If  ) 

^  Jf  ~  &YtE(x)  H  L 


(IV-3) 


where  <J>  =  e  /E(x)dx  is  the  electron  potential  energy  and  £  =imv  2  is 

A 

the  electron  kinetic  energy  in  the  field  direction.  By  defining  the 
total  kinetic  energy  as  e  =  £  +  h  ,  then  f  can  be  represented  as  the  sum 
of  two  functions,  f_(  5  $  )  which  is  the  distribution  of  electrons 

with  v  <0  moving  against  the  field  and  f  (  e,  &  <p  )  which  is  the 

A  ' 

distribution  of  electrons  with  vx  > 0  moving  with  the  field.  The 
variable  n  represents  the  electron  random  kinetic  energy.  Equations 
IV-2  and  IV-3  can  be  used  to  mathematically  define  equilibrium  and 
nonequilibrium  regions.  Equilibrium  can  be  defined  as  occurring  when 
jj’or  when  85  ev/E6<)  ji" )c,l  •  Thus  the  kinetic  energy  gained 
from  the  field  is  balanced  by  the  loss  of  energy  to  collisions.  In 
nonequilibrium  regions,  the  decrease  in  the  potential  energy  of  the 


IV-8 


electrons  is  converted  predominantly  into  an  increase  in  their  kinetic 
energy.  In  equation  IV-2,  equilibrium  can  similarly  be  defined  as 
occurring  when  or  when  *  e£(J]  •  In  other 

words,  equilibrium  exists  when  the  distribution  function  changes  faster 
in  velocity  space  than  it  does  as  a  function  of  position.  This  was  the 
definition  used  in  Chapter  I. 

In  the  absence  of  collisions,  f  is  constant  along  the  energy 
characteristics,  C  -  4>  =  constant.  These  characteristics  correspond  to 
electron  trajectories  in  the  two-dimensional  energy  space  (£,<f>  ).  The 
right  hand  term  in  equation  IV-3  represents  the  change  in  f  along  the 
energy  characteristics  when  collisions  are  included.  This  collision 
term  is  a  complicated  function  involving  f,  and  the  cross-sections  for 
all  the  various  collision  processes.  The  theory  and  derivation  of  this 
term  is  described  in  Long's  report  (ref  60:12-21).  He  showed  that  IV-3 
can  be  put  in  operator  form 

where  K  is  an  integral  operator  representing  a  sum  over  all  collisional 
processes  and  Q(  e)  is  the  total  cross  section  for  electrons  with  total 
kinetic  energy  e.  Using  this  technique,  the  distribution  function  f 
does  not  need  to  be  stored  between  iterations,  only  Kf  must  be  re¬ 
tained.  The  collision  term  can  be  further  simplified  by  assuming  Kf  is 
nearly  isotropic  in  velocity  space  and  expanding  Kf  in  Legendre  poly¬ 
nomials,  P^(cos  0)  for  example. 

0O 

=  H  Ajlv)  R  («•* ©) 

Jt*  0  *  x 


(IV-5) 


where 


This  expansion  is  quite  different  from  the  usual  two  term  expansion  of 
Boltzmann's  equation  in  spherical  harmonics.  This  expansion  is  an 
expansion  of  the  product  of  the  integral  operator  times  f,  and  not  an 
expansion  of  f.  Assuming  that  the  secondary  electrons  are  emitted 
isotropically  and  neglecting  superelastic  electron  collisions  and  the 
recoil  term  for  electrons,  Long  derived  the  following  form  of  the 
Boltzmann  equation  in  terms  of  the  total  electron  energy: 
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{  Hr*  S  B)J(cos&) 

0/c)s  v*  j <r(Kf  i)  $(coslf)j (fotH) 


( IV— 6) 


=  inelastic  collision  energy  loss 
e  i  =  ionization  energy  loss 

(V  =  collision  cross  section  integrated  over  all  scattering 
angles 

superscript  h  refers  to  a  sum  over  inelastic  processes 
superscript  i  refers  to  a  sum  over  ionization  processes 
superscript  e  refers  to  a  sum  over  elastic  collision  processes 
superscript  a  refers  to  an  attachment 
superscript  r  refers  to  recombination 
n+  =  positive  ion  density 


IV- 10 


Thus,  in  equation  IV-4,  the  electrons  scattered-in  is  represented  by 

K  f-tih  [{Mtffe)  U<  <f)  Q*  U‘) 

*  f~t  %  (t,^l  Q(tft)dc'  j  (/(«•*  ©)  (IV-7) 

and  the  electrons  scattered-out  is  represented  by 

TT* 

q:u)  *  $  o^fj  ( iv-8) 

Equation  IV-4  is  solved  numerically  using  an  iterative  technique 
for  evaluating  Kf.  It  is  assumed  initially  that  there  is  no  contribu¬ 
tion  from  scattered  in  electrons,  ie  Kf  =  0  and  equation  IV-6  is  then 
integrated  both  for  forward  moving  electrons  and  backward  moving 
electrons  for  each  value  of  n=e-C.  This  yields  a  first  approximation 
to  the  distribution  function  from  which  a  new  Kf  can  be  calculated. 

The  procedure  is  repeated  until  convergence  is  reached,  typically  when 
there  is  less  than  .2%  change  in  the  distribution  function.  By  per¬ 
forming  the  integration  in  a  specific  order  the  same  array  can  be  used 
to  store  both  the  old  and  the  new  values  of  Kf,  thus  minimizing  com¬ 
puter  storage  requirements.  The  arrays  for  Kf+  and  Kf_  are  trape¬ 
zoidal  in  shape  since  f+(e,£,4>)  =  0  for  S  >  ,  where  £  is  the 

initial  kinetic  energy  in  the  field  direction  of  the  electrons  leaving 
the  cathode.  Fig  IV-1  shows  how  the  arrays  for  Kf  are  arranged  in  core 
to  further  minimize  storage  requirements.  The  total  area  of  the  array 
in  Fig  IV-1  is  48,240jq  ( 136 , 160g)  words.  For  a  ’00V  cathode  fall, 
these  dimensions  give  a  leV  energy  resolution.  The  storage  of  only  Kf 
prevents  the  distribution  function  from  being  known  as  a  function  of 
distance,  except  at  the  cathode  and  anode  where  it  is  stored  in  sepa- 


rate  arrays. 


Boundary  conditions  are  usually  specified  at  the  cathode  and 
anode.  At  the  cathode,  the  energy  distribution  of  electrons,  f+(e,C,0) 
is  a  variable  input.  The  function  f_(e, £;,<)>)  can  represent  the  dis¬ 
tribution  of  electrons  at  the  head  of  the  positive  column  when 
<P  =  <J>  „  =  V  .  the  voltage  drop  across  the  cathode  fall,  negative  glow 

and  Faraday  dark  space,  or  it  can  represent  the  distribution  of  elec¬ 
trons  at  the  anode  then  f_  =  0,  or  it  can  represent  the  distribution 
along  an  axis  of  symmetry  such  that  f_  =  f  .  Along  the  locus  of 
turning  points  where  the  electron  kinetic  energy  in  the  field  direction 
£=  0,  the  forward  and  backward  fluxes  are  equal,  f+(e ,0,d> )  =  f _ (e , 0 ,4> ) . 

The  algorithm  used  in  Long's  code  for  the  numerical  integration  of 
equation  IV-6  is  second  order  and  stable.  The  distribution  function  is 
calculated  via 


(IV-S) 


where  c  ' 
each  step  is 


and  is  the  energy  step. 


The  relative  error  at 

(IV-10) 


which  implies  that  the  error  is  small  as  long  as 


Near  the  cathode  this  condition  is  easily  satisfied  since  «£.. 
which  is  typically  20eV  or  less. 

In  the  negative  glow  and  Faraday  dark  space  however  the  field  be¬ 
comes  very  small.  In  this  region  of  the  discharge,  the  distribution 
function  is  changing  more  rapidly  with  respect  to  position  (x)  than  it 
is  with  respect  to  potential  (<j>  ).  Thus  $  is  no  longer  an  appropriate 
variable,  and  x  becomes  the  variable  of  choice.  Changing  variables, 
equation  I V-4  becomes 

77  *  eEM  jy  =  (K'QM)f  w-i*) 

To  correctly  determine  f  through  these  regions,  the  code  first  starts 
with  coordinates  (e,  F  4>)  and  iterates  on  equation  I V-5 .  When 
A£>2E/NQ,  then  the  coordinates  are  changed  to  (e,  £,x)  and  equation 
IV- 12  is  then  iterated  on  until  convergence. 

The  code  is  presently  configured  to  calculate  the  electron  energy 
distribution  function  as  a  function  of  distance  from  the  cathode  to  the 
anode.  The  electrons  are  constrained  to  move  either  with  or  against 
the  field,  i.e.  all  integrations  are  in  the  plane  n  =  0.  This  simpli¬ 
fication  significantly  reduces  computer  run  time. 

The  electron  number  density  is  normalized  to  the  electron  number 
density  at  the  cathode  such  that 


where  subscript  0  indicates  the  boundary  condition  at  the  cathode. 


Similarly  the  electron  current  density  is  normalized  to  the  electron 
current  density  at  the  cathode  and  is  given  by  the  following 


Townsend  rates  for  the  various  inelastic  processes  including  ionization 
are  normalized  to  he  electron  current  density  at  the  cathode  as  well 


°c(4)-- N($ 


( I V- 1 5} 


In  order  to  obtain  a  completely  self-consistent  solution  in  which 
Poisson's  equation  is  obeyed  in  addition  to  the  boundary  conditions  at 
the  cathode  and  anode,  successive  calculations  for  the  electron  dis¬ 
tribution  function,  electron  number  density  and  the  resulting  electric 
field  were  accomplished. 

Long  attempted  to  do  this  for  the  case  where  the  anode  is  just 
inside  the  cathode  fall-negative  glow  boundary.  By  integrating  the 
electron  distribution  function  he  obtained  the  electron  number  density 
and  electron  current  density  and  obtained  the  positive  ion  number 
density  using  current  conservation. 
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In  this  calculation,  however  he  had  to  use  his  initial  guess  for  the 
electric  field,  so  the  ion  drift  velocity  (p+E)  was  not  consistent  with 
the  new  calculation  of  the  field. 

He  normalized  his  number  densities,  not  through  the  previous  defi¬ 
nitions,  but  by  the  ratio  of  Eq/E ^  where 


E  =  electric  field  at  the  cathode  prior  to  the  calculation 
o 

of  the  distribution  function 


** 

Ei  =  /("a.  -  n  )dx  which  is  the  electric  field  at  the  cathode 
after  the  unnormalized  densities  are  known.  Xj  is  the 
spatial  point  at  which  the  field  is  zero. 

This  approach  does  not  permit  the  secondary  emission  coefficient  to  be 
used  as  an  input  parameter,  however  the  coefficient  was  calculated  via 


the  following  equation: 


Je„ 


(IV-17) 


This  normalization  gave  good  results  and  converged  as  long  as  the  anode 
was  at  the  boundary  or  within  the  cathode  fall  region. 


It  was  found  by  this  study  that  Long's  technique  diverged  when  the 
anode  was  more  realistically  placed  outside  the  cathode  fall,  so  a  new 


technique  to  reach  a  self-consistent  solution  was  derived.  This 


technique  is  based  on  putting  Poisson's  equation  in  difference  form  so 
the  value  of  the  field  at  each  point  is  based  on  the  value  at  the 
previous  increment  with  the  field  at  the  cathode  given  as  a  boundary 
condition.  Writing  Poisson's  equation  in  difference  form  yields. 


( IV-18) 


where  n  . ,  n  .,  n  .  are  the  positive  ion,  electron,  and  negative  ion 

'  J  ^  J  ”  J 

number  densities  at  position  j,  and  x.  is  the  spatial  point  correspond- 

J 

ing  to  increment  j.  The  positive  ion  number  density  can  be  written  in 
terms  of  the  current  conservation  such  that 
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( IV-19) 


The  negative  ion  number  density  can  be  calculated  from  the  negative  ion 
current  density. 


J-j  ~  fl,  Jtj 


( IV-20) 


( IV-21 ) 


where  n  is  the  attachment  rate.  Equations  IV-21  and  IV-19  can  be 
substituted  into  IV-18.  This  results  in  a  fifth  order  polynomial  whose 
solution  can  be  found  by  numerical  means. 


Steinman's  method  (ref  66)  was  used  in  this  study  to  solve  for  the 
roots.  Thus  E  is  calculated  based  on  ng,  je,  j_  and  the  previous  value 
of  the  electric  field.  This  technique  converged  after  successive 
iterations  of  first  calculating  f(  e),  ng,  j  and  then  E(x).  Fig  IV-2 
gives  an  example  of  the  convergence  as  the  electric  field  and  electron 
distribution  are  successively  iterated. 

The  initial  discharge  parameters,  Eq,  Em.n,  and  dc>  were  varied  to 
achieve  the  best  convergence.  These  parameters  describe  the  initial 
field  distribution  according  to  the  formulas 


X  <  J. 

(IV-23a) 

EM -EL. 

(IV-23b) 

It  was  observed  that  the  optimum  convergence  was  achieved  when  the 
initial  E^n  was  slightly  greater  than  (what  was  later  found  to  be)  the 
converged  Em^n  and  the  initial  slope  EQ/dc  was  equal  to  or  slightly 
shallower  than  the  shallowest  slope  of  the  converged  curve.  This 
generally  was  the  slope  occurring  between  Eq  and  3EQ/4  or  7EQ/8. 

Fig  IV-3  shows  the  effect  of  changing  the  value  of  the  electric 
field  at  the  cathode.  Note  that  if  EQ  is  assumed  too  low,  In  this  case 
for  E  <757v/cm,  then  the  solution  to  equation  IV-22  becomes  imagi- 


FIG  TV-2  CONVERGENCE  OF  THE  SHEATH  CODE  TO  A  SELF-CONSISTENT  ELECTRIC  FIELD 


CONVERGENCE  OF  THE  SHEATH  CODE  TO  THE  LOWEST  VOLTAGE 


nary,  representing  an  unphysical  solution  for  the  field  between  the 
electrodes.  This  corresponds  quite  closely  to  the  phenomena  seen  in 
the  convergence  of  the  previous  GLOW  model  (Fig  1 1 1-2 ) .  The  calcu¬ 
lation  for  EQ  *  758v/cm  was  deemed  to  be  the  solution  most  representa¬ 
tive  of  a  discharge  since  it  represented  the  lowest  voltage  in  which  a 
complete  real  solution  existed  from  cathode  to  anode.  Although  EQ  was 
varied  manually,  a  technique  similar  to  the  modified  linear  interpo¬ 
lation  or  halving  the  interval  technique  could  be  programmed  to  con¬ 
verge  to  the  lowest  value  of  EQ  for  which  a  complete  solution  exists. 

Summarizing  the  changes  made  since  Long  published  his  report,  this 
study  has  added  attachment  as  an  electron  loss  process  in  the  numerical 
calculation  of  the  distribution  function,  modified  the  code  to  Include 
a  variable  number  of  gases,  and  modified  the  code  to  include  a  variable 
number  of  collision  processes  per  gas.  The  original  program  was 
limited  to  analyzing  a  pure  gas  with  one  momentum  transfer,  one  exci¬ 
tation,  and  one  ionization  cross  section.  The  Townsend  ionization  and 
attachment  rates  are  now  calculated  from  the  distribution  function. 
Instead  of  a  net  ionization  rate  from  the  electron  current  continuity 
equation.  In  addition,  in  the  calculation  of  the  positive  ion  drift 
velocity.  Long  had  assumed  It  was  proportional  to  the  field.  This  has 
been  changed  so  that  the  positive  ion  drift  velocity  is  more  appropri¬ 
ately  represented  at  high  fields  by  equation  I 11-41  in  which  it  is 
proportional  to  the  square  root  of  the  field.  Also  a  new 
self-consistent  technique  was  found  for  calculating  the  electric  field 
from  the  electron  density,  the  electron  current  density,  and  the 
negative  Ion  current  density. 

The  SHEATH  program  was  normally  run  on  a  Cyber  74  computer  system. 


Each  Iteration  to  converge  to  a  solution  for  a  given  EQ  took  an  average 
of  15sec.  Normally  4  to  6  iterations  were  required  to  reach  three 
significant  figure  convergence  in  E(d).  The  code  required  216K  of 
memory  to  load  and  run.  The  majority  of  this  space  was  taken  up  by  the 
220  x  201  array  described  earlier. 

Results  of  Numerical  Calculations 

This  section  presents  the  results  of  calculations  using  the  SHEATH 
code  in  the  same  rare  gas-HCl  mixtures  as  in  the  previous  chapter.  All 
the  data  presented  was  obtained  for  a  gas  pressure  of  1  torr.  Again 
the  results  for  He  will  be  be  presented  first,  followed  by  Ar  mixtures, 
and  then  Xe  mixtures. 

Only  the  cathode  fall  region  will  be  discussed  because  the  anode 
was  placed  just  outside  the  cathode  fall  region  where  the  electric 
field  leveled  off  to  a  constant  value  as  in  an  obstructed  discharge. 

The  SHEATH  code,  while  improved,  still  had  difficulty  converging  at 
very  low  electric  fields  where  the  code  transformed  from  integrations 
over  energy  to  integrations  over  position.  This  transformation  was 
discussed  in  the  preceding  section.  All  the  SHEATH  data  presented  in 
this  study  was  obtained  before  this  transformation  of  variables 
occurred,  1e.  while  integrations  are  conducted  in  energy  space. 

Further  comments  on  the  anode  fall  will  be  made  in  Chapter  V. 

Helium 

This  section  presents  the  results  of  the  nonequilibrium  calcula¬ 
tions  for  He/HCl  mixtures.  A  sunmary  is  given  in  Table  IV-1  of  the 
same  scaling  parameters  listed  In  tables  in  the  previous  chapter.  The 


Table  IV-1 


i  Data  was  taken  from  a  summary  presented  in  ref  6.  iv  Data  was  taken  from  ref  20 

ii  Data  was  taken  from  ref  86  (not  self-consistent).  (not  selt-consi stent) 

iii  Data  was  taken  from  ref  79.  v  Data  was  taken  from  ref  21 

vi  Data  was  calculated  in  this  study 


cathode  fall  width  and  voltage  was  calculated  from  the  electric  field 
in  the  same  manner  as  in  Chapter  III  for  the  GLOW  model. 

There  is  a  wide  range  of  values  for  SF(y,  k+)  in  Table  IV— 1 .  Some 
of  these  differences  are  due  to  the  differences  in  input  data,  such  as 
the  Townsend  secondary  emission  coefficient  or  the  cross  section  set 
used.  Also,  there  were  different  approximations  used  for  the  positive 
ion  drift  velocity  dependence  on  the  electric  field.  However,  there 
were  also  important  model  differences  that  influenced  the  results  as 
well.  The  data  from  calculations  by  Boeuf  and  Marode  (ref  20)  and  Tran 
Ngoc  et  al.  (ref  86)  are  not  self-consistent.  They  picked  values  of 
pd,  J/p2,  and  Vc  which  agreed  well  with  experiment  and  assumed  a 
linearly  decreasing  field  before  calculating  the  electron  distribution 
function  as  a  function  of  distance  through  the  cathode  fall  using  a 
Monte  Carlo  technique.  Segur  et  al.  (ref  79)  calculated  the  positive 
ion  density  and  current  from  the  Boltzmann  equation  for  ions,  so  an 
analytical  fit  to  experimental  data  was  not  used  for  the  positive  ion 
drift  velocity.  Boeuf  et  al.  (ref  21)  used  a  Monte  Carlo  technique  for 
both  electrons  and  positive  ions.  The  positive  ion  "nonequilibrium 
drift  velocity"  predicted  by  the  Monte  Carlo  technique  was  within  20% 
of  experimental  data  which  had  been  measured  in  equilibrium.  In  the 
SHEATH  code,  equation  I 11-41  was  used  to  calculate  the  positive  ion 
drift  velocity.  This  analytic  approximation  agrees  more  closely  with 
the  drift  velocity  predicted  by  the  Monte  Carlo  results  of  Boeuf  et  al. 
than  the  experimental  data  they  presented.  As  a  result,  the  scaling 
factor  SF(y,k+)  for  the  data  from  the  Monte  Carlo  technique  and  the 
SHEATH  code  agree  within  1%.  This  agreement  with  the  Monte  Carlo 
technique  was  used  as  a  benchmark  in  comparing  the  accuracy  of  the 
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SHEATH  model  to  other  nonequil ibrium  techniques.  Note  that  SF(y,  k+) 
changes  very  little  for  the  different  He/HCl  mixtures,  indicating  that 
SF(Y,  k+)  is  still  an  invariant  even  when  the  ionization  efficiency  of 
the  gas  changes.  These  changes  in  the  ionization  coefficients  will  be 
seen  in  later  figures.  The  observation  that  SF  (Y,  k+)  is  not  a 
function  of  the  ionization  coefficient  is  consistent  with  equation 
I I 1-72 .  The  variation  in  SF  (Y,  k+)  in  this  author's  results  is  due  to 
round  off  and  extrapolation  errors  in  pdc  and  Vc. 

The  SHEATH  program  predicts  a  15%  decrease  in  the  electric  field 
in  the  cathode  fall  region  as  the  amount  of  HC1  is  increased  form  0  to 
5%.  See  Fig  IV-4.  This  corresponds  to  a  32%  decrease  in  the  cathode 
fall  voltage.  This  trend  agrees  with  the  results  of  the  GLOW  model 
including  the  linearity  of  the  fjeld  near  the  cathode  in  both  cases. 

The  shift  in  the  electric  field,  however,  is  much  greater  in  the  SHEATH 
model  and  can  be  observed  in  the  GLOW  model  in  Fig  III -4.  A  more 
detailed  comparison  of  the  two  techniques  will  be  described  in  the  next 
chapter.  This  decrease  in  the  field  and  voltage  as  a  function  of  HC1 
is  due  to  the  increase  in  the  Townsend  ionization  coefficient  which  can 
be  seen  in  Fig  IV-5.  The  increase  in  the  peak  region  which  corresponds 
to  the  cathode  glow  region  is  due  to  the  fact  that  HC1  is  much  easier 
to  ionize  than  He.  HC1  has  an  ionization  threshold  of  12.74eV  which  is 
about  half  of  that  for  He  whose  ionization  threshold  is  24.59eV.  Note 
that  the  spatial  threshold  for  the  Townsend  ionization  coefficient  is 
halved  when  HC1  is  added  corresponding  to  the  change  in  the  ionization 
cross  section  threshold.  This  indicates  that  very  close  to  the  cath¬ 
ode,  the  HC1  is  being  ionized  before  the  He.  These  results  are  consis¬ 
tent  with  the  results  of  Segur,  et  al.  with  Hg  in  He.  The  decrease  in 
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the  linear  electric  field  with  the  addition  of  HC1  results  in  similar 
decreases  in  the  cathode  fall  length,  decreasing  6.4%  for  the  99/1 
mixture  and  21%  for  the  95/5  mixture.  These  results  are  summarized  in 
Table  IV-1. 

The  curves  for  the  various  electric  fields  in  Fig  IV-4  are  almost 
parallel  because  the  slope  of  the  electric  field  is  determined  through 
Poisson's  equation  by  the  positive  ion  density.  The  positive  ion 
density  at  the  cathode  is  determined  by  y,  0,  and  the  positive  ion 
drift  velocity.  Since  these  variables  and  relations  remained  constant 
for  the  various  HC1  mixtures,  the  slopes  of  the  electric  fields  re¬ 
mained  constant. 

Fig  IV-6  and  IV-7  compare  the  Townsend  ionization  and  attachment 
coefficients  as  a  function  of  distance  for  99/1  and  95/5  concentrations 
of  He/HCl.  The  magnitude  of  the  attachment  rate  is  much  smaller  than 
the  ionization  rate  throughout  the  cathode  fall  region  as  one  would 
expect,  except  immediately  in  front  of  the  cathode.  In  this  region 
which  corresponds  to  the  Primary  dark  space,  the  electrons  have  just 
left  the  cathode  at  a  few  eV  of  energy  and  are  not  yet  in  equilibrium 
with  the  field.  The  second  peak  seen  in  the  attachment  curve  at  l.cm 
is  probably  due  to  the  electron  distribution  function  first  reaching 
the  ionization  energy  and  producing  a  new  crop  of  slow  secondary  elec¬ 
trons. 

Fig  IV-8  and  IV-9  display  the  electron  and  negative  ion  current  dens¬ 
ities  as  a  function  of  distance  and  as  a  function  of  percent  of  HC1 . 

The  increase  in  the  ionization  rate  also  manifests  itself  in  an  increase 
in  the  electron  current  density  throughout  the  cathode  fall  region.  The 
electron  current  density  stays  constant  for  the  first  few  mean  free 
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FIG  IV- 7  COMPARISON  OF  IONIZATION  AND  ATTACHMENT  COEFFICIENT  IN  95/5  He/HCl 
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paths  as  the  electrons  are  being  accelerated  up  to  the  ionization  poten¬ 
tial.  It  then  increases  almost  exponentially  through  the  remainder  of 
the  cathode  fall.  As  expected,  the  negative  ion  current  density  in  Fig 
IV-9  increases  as  the  proportion  of  HC1  increases.  It  is  characterized 
by  a  very  steep  onset  due  to  the  large  attachment  rate  near  the  cathode. 
As  the  total  distribution  of  electrons  is  accelerated  away  from  the 
cathode,  they  reach  the  ionization  threshold  of  the  gas.  Most  electrons 
are  then  past  the  peak  region  of  the  attachment  cross  section  and  the 
negative  ion  current  density  grows  at  a  much  slower  rate.  Note  the 
logarithmic  scale,  indicating  the  growth  of  the  negative  ion  current 
density  is  still  slightly  faster  than  an  exponential  function  of  dis¬ 
tance. 

Fig  IV-10,  IV-11,  and  IV-12  display  the  electron,  positive  ion,  and 
negative  ion  number  densities  as  a  function  of  distance  through  the 
cathode  fall  region.  The  result  of  nonequilibrium  electron  kinetics  is 
easily  seen  in  these  figures.  The  electron  density  initially  decreases 
because  the  electrons  are  rapidly  accelerated  away  from  the  cathode. 

Just  before  they  have  gained  sufficient  energy  for  ionization,  their 
density  reaches  a  minimum  and  begins  to  increase.  The  fact  that  by 
mid-cathode  fall  the  field  has  also  decreased  substantially  allows  the 
electron  density  to  increase.  Similarly  the  positive  ion  density  in¬ 
creases  away  from  the  cathode  since  there  is  very  little  ionization 
close  to  the  cathode  and  the  positive  ions  are  being  accelerated  towards 
it.  Although  negative  ions  are  present  throughout  the  cathode  fall 
region  in  Fig  IV-11  and  IV-12,  their  density  is  generally  much  too  small 
to  affect  the  net  space  charge  to  any  significant  degree.  Note  that  the 
negative  ion  density  is  more  than  four  orders  of  magnitude  smaller  than 
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FIG  IV-10  CHARGED  PARTICLE  DENSITIES  IN  100’.  He 
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FIG  IV-11  CHARGED  PARTICLE  DENSITIES  IN  99/1  He/HCl 
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FIG  IV-12  CHARGED  PARTICLE  DENSITIES  IN  95/5  He/HCl 


the  positive  ion  density.  The  positive  ion  density  actually  determines 
the  slope  of  the  field  throughout  the  cathode  fall.  The  magnitude  of 
the  electron  density  approaches  that  of  the  positive  ion  density  only 
after  the  negative  glow  is  reached  or,  in  these  examples  just  prior  to 
the  anode.  Equilibrium  is  never  reached  in  these  figures  because  a 
slightly  bimodal  electron  distribution  still  exists  at  the  anode. 

Since  the  cathode  fall  voltage  decreases  as  HC1  is  added,  the  ions 
arrive  at  the  cathode  at  a  lower  velocity.  Similarly,  the  electrons  are 
not  accelerated  away  as  quickly,  so  both  electron  and  positive  ion  dens¬ 
ities  increase  near  the  cathode  (20%  and  8%  respectively  for  up  to  5% 
HC1).  The  opposite  trend  can  be  seen  in  Fig  IV-10,  IV-11,  and  IV-12 
near  the  cathode  glow-negative  glow  boundary  (1.5-1.75  cm).  The  elec¬ 
tron  density  increases  as  a  function  of  HC1  in  this  region  due  to  the 
increase  in  the  ionization  coefficient  closer  to  the  cathode.  Since 
the  field  is  lower  in  this  region  and  the  ionization  coefficient  is 
lower,  the  positive  ion  density  decreases  as  the  percentage  of  HC1  is 
increased,  which  in  a  self-consistent  (or  interrelated)  manner  causes 
the  cathode  fall  region  to  contract. 

Argon 

The  results  for  Ar/HCl  mixtures  are  summarized  in  Table  IV-2 
as  was  done  for  previous  mixtures.  As  in  He,  there  is  a  wide  range  of 
values  for  SF(y,  k+)  in  Table  IV-2  for  Ar.  The  only  other  nonequili¬ 
brium  data  available  for  Ar  is  that  of  Long  (ref  60)  which  was  not 
entirely  self-consistent.  The  lack  of  self-consistency  of  his  method 
was  discussed  earlier  in  this  Chapter.  Also,  he  had  assumed  the  posi¬ 
tive  ion  drift  velocity  was  proportional  to  the  field.  This  assumption 
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i  Data  was  taken  from  a  summary  presented  in  ref  6. 

ii  Data  was  taken  from  ref  60  (not  self-consistent). 

iii  Data  was  calculated  in  this  study. 


is  only  appropriate  in  the  low  field  region,  i.e.  E/p<60v/cm-torr. 

Note  that  as  in  He  mixtures  the  scaling  factor  SF(y.  k+)  remains  almost 
constant  as  the  concentration  of  HC1  increases. 

In  contrast  to  He  mixtures,  there  is  a  negligible  effect  on  the 
electric  field  in  the  cathode  fall  region  when  up  to  5%  HC1  is  added. 
This  can  be  seen  in  Fig  IV-13  and  these  results  agree  with  the  GLOW 
model.  This  electric  field  distribution  through  the  cathode  fall  cor¬ 
responds  to  approximately  a  184v  potential  drop  across  the  cathode  fall 
region.  With  1%  HC1 ,  the  difference  in  total  ionization  is  too  small 
to  have  any  observable  effect  on  the  electric  field.  The  ionization 
thresholds  for  Ar  (15. 7eV )  and  HC1  (12. 7eV )  are  much  closer  than  for  He 
and  HC1.  This  difference  is  too  small  to  be  observed  in  Fig  IV-14  in 
the  Townsend  ionization  coefficient.  With  5%  HC1 ,  there  is  sufficient 
HC1  to  increase  the  ionization  rate  about  1.6%  throughout  most  of  the 
cathode  fall  region.  Fig  IV-15  and  IV-16  compare  the  Townsend  ioniza¬ 
tion  and  attachment  coefficients  as  a  function  of  distance  for  99/1  and 
95/5  concentrations  in  Ar/HCl.  As  in  He/HCl  mixtures,  the  attachment 
rate  is  larger  than  the  ionization  rate  only  for  a  few  electron  colli¬ 
sions  close  to  the  cathode  where  the  electrons  have  insufficient  energy 
for  ionization. 

Fig  IV- 17  and  IV- 18  display  the  electron  and  negative  ion  current 
densities  as  a  function  of  distance  and  concentration  of  HC1.  Again 
very  little  change  is  observed  in  the  electron  current  density  since 
the  change  In  HC1  concentration  had  only  a  minor  effect  on  the  ioniza¬ 
tion  coefficients  in  Fig  IV-14.  The  electron  current  density  stays 
constant  as  the  electrons  are  accelerated  up  to  the  ionization  poten¬ 
tial.  Again  the  negative  ion  current  density  in  Fig  IV-18  i.  charac- 
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FIG  IV- 13  ELECTRIC  FIELD  IN  Ar/HCl  MIXTURES 
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FIG  IV- 14  COMPARISON  OF  TOWNSEND  IONIZATION  COEFFICIENTS  IN  Ar/HCl  MIXTURES 
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FIG  IV-17  ELECTRON  CURRENT  DENSITY  IN  Ar/HCl  MIXTURES 


terized  by  a  very  steep  onset  followed  by  a  slightly  faster  than  expo¬ 
nential  growth  through  most  of  the  cathode  fall  region.  Even  though 
the  negative  ion  current  density  increases  with  HC1  concentration,  its 
contribution  to  the  total  discharge  current  through  the  cathode  fall  is 
still  negligible. 

The  same  nonequilibrium  phenomena  that  occurred  in  He  mixtures  in 
the  number  densities  can  be  observed  in  Ar  mixtures  in  Fig  IV-19,  I V-20 
and  IV-21.  The  positive  ion  density  increases  in  the  direction  moving 
away  from  the  cathode,  while  the  electron  number  density  initially 
decreases.  The  number  of  negative  ions  is  more  than  two  orders  of 
magnitude  smaller  than  the  number  of  positive  ions,  and  actually  their 
density  does  not  begin  to  grow  exponentially  until  the  cathode  fall 
begins  to  merge  into  the  negative  glow. 

Since  there  is  very  little  difference  in  the  fields  adjacent  to 
the  cathode,  there  is  very  little  difference  in  the  electron  and  posi¬ 
tive  ion  number  densities  in  this  region.  The  positive  ion  density  is 
still  several  orders  of  magnitude  larger  than  the  electron  and  negative 
ion  densities  and  thus  determines  the  slope  of  the  electric  field. 

Only  in  the  last  third  of  the  cathode  fall  length  (.5-. 75  cm)  is  there 
a  change  when  Fig  IV-19,  IV-20,  and  IV-21  are  compared.  As  in  He/HCl 
mixtures,  the  electron  number  density  increases  and  the  positive  ion 
number  density  decreases  but  to  a  lesser  degree  since  there  is  less 
change  in  the  Townsend  ionization  coefficient. 

Xenon 

The  results  for  Xe/HCl  mixtures  are  summarized  in  Table  IV-3. 


FIG  IV-19  CHARGED  PARTICLE  DENSITIES  IN  100%  Ar 
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FIG  IV-20  CHARGED  PARTICLE  DENSITIES  IN  99/1  Ar/HCl 


There  is  very  little  data  for  comparison  in  Xe.  The  cathode  fall  volt¬ 
age  Vc  and  length  dc  are  sensitive  functions  of  the  secondary  emission 
coefficient  Y.  The  value  of  Y  =  .004  was  chosen  not  to  compare  with 
an  iron  cathode  but  to  be  able  to  compare  with  the  equilbrium  results 
of  the  GLOW  code  which  were  compared  with  Ward's  data.  Even  though 
this  value  of  Y  is  low,  it  can  not  totally  account  for  the  large 
discrepancy  between  theory  and  experiment  in  SF(Y,k+).  Of  the  three 
rare  gases  studied,  the  cross  sections  for  Xe  are  probably  the  least 
well  known.  The  derivation  of  the  set  of  cross  sections  for  Xe  is 
described  in  Appendix  C.  The  cross  sections  could  be  in  error  by  as 
much  as  20-25%.  Since  a  parametric  study  of  the  cross  section  thresh¬ 
olds  and  initial  slopes  was  not  a  part  of  this  study,  the  magnitude  of 
any  change  in  the  cross  sections  is  as  yet  unknown.  Since  SF(Y,k+)  is 
larger  for  the  theoretical  data,  this  implies  that  a  smaller  pdc  is 
required  to  make  the  theoretical  data  agree  more  closely  with  experi¬ 
ment.  A  reduction  in  pdc>  probably  will  also  be  accompanied  by  a  de¬ 
crease  in  Vc<  Changing  the  excitation  and  ionization  cross  sections  in 
order  to  make  the  Townsend  ionization  coefficient  larger  would  tend  to 

make  pd  and  V  smaller.  Another  source  of  error  could  be  in  the  value 
c  c 

of  k+.  The  value  of  k+  was  obtained  from  Ward  (ref  92:2791)  who  sub¬ 
jectively  fit  equations  1 1 1 -71  to  experimental  data.  He  made  no  sys¬ 
tematic  effort  to  obtain  the  best  parameters.  In  addition,  only  the 
first  term  in  equation  I I 1-71 b  was  used  by  the  SHEATH  model.  In 
effect,  a  smaller  positive  ion  drift  velocity  would  result  in  a  steeper 
slope  for  the  electric  field.  This  also  would  tend  to  make  pdc  and  Vc 
smaller.  The  minor  variation  in  S F (  Y, k  + )  in  this  author's  results  is 
due  to  round  off  and  extrapolation  errors  in  pdc  and  Vc. 


In  contrast  to  He  and  Ar  mixtures,  there  is  a  slight  increase  in 
the  electric  field  in  the  cathode  fall  region  when  up  to  5%  HC1  is 
added.  This  trend  of  the  field  shifting  vertically  can  be  seen  in  Fig 
IV-22.  It  results  in  a  slight  increase  in  the  voltage  which  is  sum¬ 
marized  in  Table  IV-3.  This  trend  is  also  the  same  as  in  Fig  III-ll 
from  the  GLOW  model.  Slightly  higher  fields  are  required  to  obtain  the 
same  net  ionization  near  the  cathode.  This  results  in  a  uniform  ver¬ 
tical  shift  in  the  electric  field  throughout  the  cathode  fall  region. 
Note  that  even  though  the  field  is  shifted  to  slightly  higher  values, 
there  is  no  observable  change  in  the  ionization  coefficient  in  the 
cathode  fall  region  as  illustrated  in  Fig  IV-23.  Again  only  near  the 
minimum  field  region  in  the  negative  glow  is  a  difference  observed. 
Xe/HCl  mixtures  are  unique  in  that  both  constituents  have  almost  the 
same  ionization  potential  (12.74  for  HC1  and  12.13  for  Xe)  and  similar 
cross  section  magnitudes.  A  description  of  these  cross  sections  occurs 
in  Appendix  C.  Thus  the  slight  expansion  seen  here  is  probably  due  to 
the  formation  of  negative  ions  close  to  the  cathode  where  there  is  a 
peak  in  the  attachment  rate.  A  loss  of  electrons  close  to  the  cathode 
is  similar  to  the  effect  of  reducing  the  secondary  emission  coefficient 
of  the  cathode.  A  decrease  in  Y  similarly  results  in  a  higher  cathode 
fall  voltage.  However,  the  slope  of  the  electric  field  stays  the  same 
because  the  positive  ions  still  dominate  the  other  number  densities. 
Thus  the  trend  of  the  field  shifting  vertically  with  an  increase  in 
voltage  would  be  expected  to  continue  for  larger  concentrations  of  HC1 
in  Xe.  In  this  study  the  positive  ion  velocity  of  only  the  rare  gas 
was  included.  For  Xe/HCl  mixtures,  the  addition  of  HC1  would  tend  to 
increase  the  average  positive  ion  drift  velocity.  This  would  result  in 
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FIG  IV-22  ELECTRIC  FIELD  IN  Xe/HCl  MIXTURES 
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FIG  IV-23  COMPARISON  OF  TOWNSEND  IONIZATION  COEFFICIENT  TN  He/HCl  MIXTURES 


a  shallower  slope  of  the  electric  field  and  a  further  expansion  in  pdc> 

Fig  IV-24  and  IV-25  compare  the  Townsend  ionization  and  attachment 
coefficients  as  a  function  of  distance  for  99/1  and  95/5  concentrations 
in  Xe/HCl.  Again  the  attachment  coefficient  is  much  smaller  than  the 
ionization  rate  throughout  the  cathode  fall  region  except  immediately 
in  front  of  the  cathode.  Note  that  the  region  where  the  attachment 
coefficient  is  greater  than  the  ionization  rate  is  smaller  in  Xe  mix¬ 
tures  than  in  He  or  Ar  mixtures.  This  is  due  to  Xe  having  a  lower 
ionization  threshold  than  the  other  buffer  gases.  The  first  peak  in 
the  attachment  coefficient  is  a  result  of  the  slow  electrons  leaving 
the  cathode  and  the  second  peak  is  a  result  of  the  formation  of  slow 
electrons  that  have  lost  their  energy  to  excitation  or  ionization.  The 
first  peak  has  the  effect  of  reducing  the  secondary  emission  coeffic¬ 
ient  of  the  cathode. 

Fig  IV-26  and  IV-27  display  the  electron  and  negative  ion  current 
densities  as  a  function  of  distance  and  concentration  of  HC1 .  Very 
little  change  is  observed  in  the  electron  current  density  until  near 
the  cathode  glow-negative  glow  boundary.  Here  the  electron  current 
density  is  observed  to  decrease  instead  of  increasing  as  a  function  of 
HC1  concentration  as  in  previous  rare  gas  mixtures.  This  is  another 
example  of  the  effect  of  attachment  in  the  cathode  •fall  region.  Even 
though  the  negative  ion  current  density  in  Fig  IV-27  increases  signifi¬ 
cantly  with  HC1  concentration,  its  contribution  to  the  total  discharge 
current  is  still  negligible. 

The  same  nonequilibrium  phenomena  that  occurred  in  previous  rare 
gas  mixtures  for  the  number  densities  can  be  observed  in  Xe  mixtures. 
The  results  are  illustrated  in  Fig  IV-28,  through  IV-30.  Again  the 


FIG  IV-24  COMPARISON  OF  IONIZATION  AND  ATTACHMENT  COEFFICIENT  IN  99/1  Xe/HCl 


COMPARISON  OF  IONIZATION  AND  ATTACHMENT  COEFFICIENT  IN  95/5  Xe/HCl 
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FIG  IV- 2 6  ELECTRON  CURRENT  DENSITY  IN  Xe/HCl  MIXTURES 
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NEGATIVE  ION  CURRENT  DENSITY  IN  Xe/HCl  MIXTURES 
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FIG  IV-28  CHARGED  PARTICLE  DENSITIES  IN  100  %  Xe 
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-29  CHARGED  PARTICLE  DENSITIES  IN  99/1  Xe/HCl 
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FIG  IV-30  CHARGED  PARTICLE  DENSITIES  IN  95/5  Xe/HCl 


trends  of  each  of  the  number  densities  is  similar  to  the  trends  seen 
before.  The  positive  ion  density  is  still  several  orders  of  magnitude 
larger  than  the  electron  and  negative  ion  densities  and  thus  determines 
the  slope  of  the  electric  field  in  Xe  as  well.  The  number  of  negative 
ions  is  more  than  three  orders  of  magnitude  smaller  than  the  number  of 
positive  ions,  and  actually  their  density  does  not  begin  to  grow  expo¬ 
nentially  until  the  cathode  fall  begins  to  merge  into  the  negative 
glow.  As  in  Ar,  there  is  almost  no  difference  in  the  number  densities 
of  electrons  and  positive  ions  adjacent  to  the  cathode,  since  the  field 
changes  only  slightly  there  when  HC1  is  added.  As  in  Ar,  only  in  the 
last  third  of  the  cathode  fall  length  (.45-. 65cm)  do  the  number  densi¬ 
ties  differ  as  a  function  of  HC1.  This  is  a  result  of  a  change  in  the 
net  ionization  prior  to  this  region. 

In  summary,  we  have  examined  three  cases  represented  by  He/HCl , 
Ar/HCl  and  Xe/HCl  mixtures.  In  all  three  gas  mixtures,  there  was  a 
peak  in  the  Townsend  attachment  coefficient  as  a  result  of  electron 
nonequilibrium.  However,  throughout  the  cathode  fall  region  the  posi¬ 
tive  ion  density  dominates  the  other  charge  densities  and  so  determines 
the  slope  of  the  electric  field  even  in  electronegative  gas  mixtures. 
He/HCl  mixtures  were  characterized  by  a  significant  contraction  of  the 
cathode  fall  length  and  a  decrease  in  the  cathode  fall  voltage.  This 
was  a  result  of  the  HC1  being  easier  to  ionize  than  the  He.  As  a 
result  of  HC1  having  a  lower  ionization  cross  section  threshold  than 
He,  the  threshold  for  the  total  Townsend  ionization  coefficient  shifted 
closer  to  the  cathode.  There  was  a  negligible  effect  in  the  electric 
field,  cathode  fall  length,  or  cathode  fall  voltage  when  HC1  was  added 
to  Ar  mixtures.  The  Townsend  ionization  coefficient  did  increase 


slightly  with  the  addition  of  5%  HC1.  This  is  because  the  ionization 
rate  is  very  sensitive  to  the  gas  species  and  the  shape  of  the  electron 
energy  distribution  function.  In  contrast,  in  Xe/HCl  mixtures,  the 
electric  field,  cathode  fall  length,  and  cathode  fall  voltage  increased 
slightly,  probably  as  a  result  of  the  nonequilibrium  electrons  being 
attached  to  form  negative  ions  in  the  primary  dark  space. 

In  conclusion,  the  formation  of  negative  ions  in  the  cathode  fall 
region  does  not  lead  to  a  contraction  of  the  sheath.  The  contraction 
that  is  often  observed  and  that  was  seen  by  Emeleus  and  Sayers,  was 
actually  due  to  an  increase  in  the  ionization  efficiency  of  the  gas 
mixture  when  an  electronegative  gas  is  added. 
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This  section  compares  and  contrasts  the  GLOW  and  SHEATH  models. 
First,  their  equations  and  assumptions  are  briefly  reviewed.  Second, 
the  results  from  these  models  are  compared.  In  order  to  determine  the 
influence  of  nonequilibrium  without  further  complications  due  to  nega¬ 
tive  ions,  results  for  the  pure  rare  gases  will  be  emphasized.  These 
results  include  the  electric  field,  Townsend  ionization  and  attachment 
coefficients,  the  electron  current  densities,  and  the  charged  particle 
densities.  The  trends  predicted  by  both  models  are  summarized  along 
with  a  criticism  of  each  technique  on  its  appl icabil ity  in  the  cathode 
fall  region.  Finally,  the  applicability  of  each  model  to  the  anode 
region  is  discussed. 


Review  of  Basic  Assumptions 

The  GLOW  and  SHEATH  models  are  founded  on  two  entirely  different 
methods  for  describing  the  kinetics  and  fields  in  a  glow  discharge.  As 
described  in  Chapter  III,  the  GLOW  model  is  based  on  a  self-consistent 
solution  to  a  set  of  equilibrium  equations:  Poisson's  equation,  the 
hydrodynamic  current  continuity  equations  for  electrons  and  negative 
ions,  and  the  current  conservation  equation.  The  SHEATH  model  des¬ 
cribed  in  Chapter  IV  is  based  on  a  solution  of  the  Boltzmann  transport 
equation  for  electrons  as  a  function  of  distance  coupled  with  Poisson's 
equation  and  the  current  conservation  equation.  The  equations  are 
listed  in  Table  V-l  for  comparison.  The  fundamental  difference  behind 


these  equations  is  that  the  glow  model  assumes  the  electrons  are  in 
equilibrium  with  the  field,  whereas  the  SHEATH  model  does  not.  Since 
both  methods  are  describing  the  same  discharge,  the  boundary  conditions 
in  Table  V-l  are  almost  identical. 

In  addition,  both  methods  predict  results  that  appear  to  scale 
according  to  equation  I 11-71.  In  the  equilibrium  model,  limited  data 
was  taken  at  10  and  100  Torr  to  check  the  pressure  scaling.  There  was 
also  good  agreement  with  other  theoretical  data  when  the  same  value  of 
the  secondary  emission  coefficient  and  the  same  relationship  of  the 
positive  ion  drift  velocity  to  the  electric  field  was  used.  In  the 
nonequilibrium  model,  the  scaling  factor  in  Table  I V-l  for  100%  He 
agreed  very  well  with  the  Monte  Carlo  results  of  Boeuf,  et  al.  (ref 
21).  Scaling  agreement  of  both  models  can  be  expected  because  both 
models  support  the  assumptions  used  in  deriving  equation  I I 1-71 ,  which 
are: 

a.  Positive  ions  dominate  the  other  charged  particle  species 
through  most  of  the  cathode  fall  and  therefore  determine  the  slope  of 
the  electric  field. 

b.  The  positive  ion  drift  velocity  is  proportional  to  the  square 
root  of  the  electric  field, 

c.  The  electric  field  can  be  approximated  as  a  linearly  decreas¬ 
ing  function  through  most  of  the  cathode  fall  region. 

Despite  the  differences  in  assuming  equilibrium  or  not  assuming 
equilibrium,  both  models  scale  appropriately  and  predict  the  same 
trends  in  all  three  gas  mixtures. 


Comparing  the  results  of  the  GLOW  and  SHEATH  models  from  Chapters 
III  and  IV  respectively  indicates  that  both  methods  do  predict  the  same 
trends  for  each  gas  mixture.  In  He  mixtures  both  models  predict  a 
contraction  of  the  cathode  fall  with  the  addition  of  HC1  as  seen  in  Fig 
III -3  and  IV-4.  Data  from  both  models  indicate  this  contraction  is  due 
to  the  more  efficient  ionization  of  the  HC1  in  the  mixture,  rather  than 
the  formation  of  negative  ions  in  the  cathode  fall  as  was  speculated  by 
Emeleus,  et  al.  (ref  42).  Similarly  in  Ar/HCl  mixtures,  both  methods 
agree  that  the  addition  of  up  to  5%  HC1  has  an  insignificant  effect  on 
the  cathode  fall  field  or  voltage.  This  can  be  observed  in  Fig  1 1 1-7 
and  IV-13.  In  Xe/HCl  mixtures  a  slight  increase  in  electric  field  as 
well  as  voltage  drop  across  the  cathode  fall  was  predicted  by  both 
models.  This  increase  can  be  seen  in  Fig  III-ll  and  IV-22.  Data  from 
both  methods  agree  that  this  field  enhancement  is  due  to  the  loss  of 
electrons  to  attachment  near  the  cathode,  although  this  is  perhaps  most 
easily  seen  from  the  theory  described  in  section,  "New  Equilibrium 
Analysis  of  the  Cathode  Fall  Region  including  Negative  Ions"  in  Chapter 
III.  This  loss  of  electrons  to  attachment  close  to  the  cathode  is 
analogous  to  reducing  the  electron  secondary  emission  coefficient. 

Fig  V-l  through  V-3  compare  the  distribution  of  the  electric  field 
through  a  discharge  for  He,  Ar,  and  Xe.  For  each  gas,  the  slope  of  the 
electric  fields  are  approximately  the  same  as  a  result  of  the  same 
field  relationship  being  used  for  the  positive  ion  drift  velocity. 

Also,  the  results  of  the  SHEATH  model  are  consistently  higher  than 
those  of  the  GLOW  model  until  the  field  begins  to  merge  into  the 
negative  glow.  This  is  a  result  of  the  nonequilibrium  behavior  of  the 


FIG  V-l  COMPARISON  OF  ELECTRIC  FIELDS  FROM  GLOW  AND  SHEATH  CODES  IN  100%  He 
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FIG  V-2  COMPARISON  OF  ELECTRIC  FIELDS  FROM  GLOW  AND  SHEATH  CODES  IN  100%  Ar 


electrons  near  the  cathode.  The  best  agreement  is  obtained  for  Xe  and 
the  worst  for  He,  since  the  electrons  have  shorter  mean  free  paths  in 
Xe  than  in  He  and  thus  spatially  reach  equilibrium  faster. 

The  Townsend  ionization  coefficient  for  each  model  is  plotted  in 
Fig  V-4  along  with  the  electric  field  calculated  by  the  SHEATH  model. 
These  coefficients  were  calculated  for  the  given  electric  field  using 
the  equations  in  Table  V-l.  As  can  be  seen,  the  plots  for  these 
coefficients  are  not  the  same  shape  due  to  the  nonequilibrium  electron 
behavior  in  the  cathode  fall  region.  In  the  SHEATH  model,  the  elec¬ 
trons  require  several  mean  free  paths  before  they  reach  the  ionization 
energy  and  can  begin  to  multiply.  This  spatial  delay  in  ionization 
results  in  a  higher  cathode  fall  voltage  and  higher  field  at  the 
cathode.  The  ionization  coefficient  predicted  by  the  SHEATH  model 
actually  lags  the  field,  ie.  the  rate  is  low  where  the  field  is  highest 
and  then  overshoots  and  is  high  where  the  field  merges  into  the  nega¬ 
tive  glow.  The  GLOW  model  assumes  the  electrons  to  be  in  equilibrium 
with  the  field  and  thus  predicts  the  highest  ionization  at  the  highest 
fields  adjacent  to  the  cathode.  Thus  the  cathode  fall  voltage  required 
to  produce  sufficient  ionization  for  the  discharge  maintenance  condi¬ 
tion  (equation  II-l)  is  lower  for  the  equilibrium  model  than  for  the 
nonequilibrium  model. 

The  Townsend  attachment  coefficient  for  each  model  is  similarly 
plotted  in  Fig  V-5  along  with  the  electric  field  calculated  by  the 
SHEATH  model.  These  coefficients  were  calculated  for  the  given  elec¬ 
tric  field  using  the  equations  in  Table  V-l.  Again  the  plots  for  these 
coefficients  are  not  the  same  shape  due  to  the  nonequilibrium  electron 
behavior  in  the  cathode  fall.  The  electrons  close  to  the  cathode  have 


FIG  V-4  COMPARISON  OF  IONIZATION  COEFFICIENTS 
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a  low  energy  and  undergo  several  collisions  before  reaching  the  ioni¬ 
zation  energy  of  the  gas.  This  results  in  a  large  peak  in  the  attach¬ 
ment  coefficient  close  to  the  cathode  in  the  results  predicted  by  the 
SHEATH  model.  The  Townsend  attachment  coefficient  predicted  by  the 
SHEATH  model  actually  leads  the  field,  i.e.  the  rate  is  high  where  the 
field  is  highest  and  then  undershoots  and  rises  again  where  the  minimum 
field  merges  into  the  negative  glow.  Note  that  the  attachment  coeffi¬ 
cient  predicted  by  the  nonequilibrium  model  is  more  than  an  order  of 
magnitude  smaller  than  that  predicted  by  the  equilibrium  model  except 
immediately  in  front  of  the  cathode.  Thus  the  ionization  and  attach¬ 
ment  coefficients  predicted  by  the  equilibrium  model  do  not  represent 
the  actual  spatial  distribution  of  the  coefficients  predicted  by  the 
SHEATH  model . 

The  growth  of  the  electron  current  density  is  portrayed  for  each 
of  the  rare  gases  in  Fig  V-6  through  V-8.  The  GLOW  model  consistently 
predicts  a  faster  rising  electron  current.  This  difference  is  again  a 
result  of  the  nonequilibrium  behavior  of  the  electrons  close  to  the 
cathode.  The  electrons  first  have  to  be  accelerated  to  the  ionization 


potential  of  the  gas  before  they  can  begin  to  multiply.  This  causes 
the  electron  current  density  predicted  by  the  SHEATH  model  to  lag  the 
current  density  predicted  by  the  GLOW  model.  To  make  the  equilibrium 
model  agree  more  closely  with  the  nonequilibrium  model,  Tran  Ngoc,  et 
al.  (ref  86)  have  suggested  that  the  mean  electron  energy  and  the 
Townsend  ionization  coefficient  be  related  through  an  energy  balance 
equation.  This  yielded  good  agreement  with  their  nonequilibrium 
ionization  coefficient  calculated  by  a  Monte  Carlo  technique.  Simi¬ 
larly,  relating  the  Townsend  ionization  coefficient  to  the  mean 
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FIG  V-6  COMPARISON  OF  ELECTRON  CURRENT  DENSITIES  IN  100%  He  FROM  GLOW  AND  SHEATH  METHODS 


FIG  V-7  COMPARISON  OF  ELECTRON  CURRENT  DENSITIES  IN  100%  Ar  FROM  GLOW  AND  SHEATH  METHODS 
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COMPARISON  OF  ELECTRON  CURRENT  DENSITIES  IN  100%  Xe  FROM  GLOW  AND  SHEATH  METHODS 


electron  energy  should  also  predict  a  slower  rising  electron  current 
density  more  closely  approximating  the  nonequilibrium  results. 

Another  example  of  the  inability  of  the  equilibrium  model  to 
describe  the  nonequilibrium  behavior  in  the  cathode  fall  region  is 
shown  in  the  plots  of  the  charged  particle  densities  in  Fig  V-9  through 
V-ll.  There  is  good  agreement  between  the  methods  on  the  positive  ion 
density  near  the  cathode.  The  GLOW  results  agree  within  13%  for  He,  6% 
for  Ar,  and  5%  for  Xe  of  the  positive  ion  density  predicted  by  the 
SHEATH  model  at  the  cathode.  This  agreement  in  positive  ion  density 
results  in  the  electric  fields  being  almost  parallel  in  Fig  V-l  through 
V-3.  However,  not  as  good  agreement  is  obtained  as  the  cathode  fall 
transitions  into  the  negative  glow  region.  This  is  due  to  the  differ¬ 
ent  rates  of  growth  of  the  electron  density.  The  GLOW  method  consis¬ 
tently  predicts  an  electron  density  that  is  higher  than  the  SHEATH 
method  even  close  to  the  cathode.  At  the  cathode  where  the  electron 
current  densities  agree  for  both  methods,  the  electron  density  cal¬ 
culated  by  the  GLOW  method  is  larger  than  the  SHEATH  results  because 
the  electron  drift  velocity  calculated  by  the  GLOW  model  is  lower  than 
that  calculated  from  the  distribution  function  in  the  SHEATH  model. 

The  difference  between  the  electron  number  densities  calculated  by  the 
GLOW  and  SHEATH  models  does  not  affect  the  slope  of  the  field  until 
near  the  negative  glow  region  because  the  electron  density  is  3-4 
orders  of  magnitude  less  than  the  positive  ion  density.  Thus  the 
equilibrium  theory  is  not  capable  of  describing  the  particle  densities 
as  a  function  of  distance  through  the  cathode  fall  region.  Through 
adjustment  of  boundary  conditions  such  as  the  secondary  emission 
coefficient  and  the  electron  drift  velocity  at  the  cathode  better 
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FIG  V-9  COMPARISON  OF  POSITIVE  ION  AND  ELECTRON  DENSITIES  IN  100%  He  FROM  GLOW  AND  SHEATH  METHODS 
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FIG  V-10  COMPARISON  OF  POSITIVE  ION  AND  ELECTRON  DENSITIES  IN  100%  Ar  FROM  GLOW  AND  SHEATH  METHODS 
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FIG  V-ll  COMPARISON  OF  POSITIVE  ION  AND  ELECTRON  DENSITIES  IN  100%  Xe  FROM  GLOW  AND  SHEATH  METHODS 


agreement  with  experimental  results  could  probably  be  obtained,  but  it 
would  still  not  represent  the  actual  nonequilibrium  kinetics  involved 
there. 

The  macroscopic  parameters  such  as  discharge  current,  cathode  fall 
voltage  and  length,  and  electric  field  each  contain  one  or  more  aver¬ 
ages  or  integrals  over  the  electron  distribution  function.  It  was 
thought  that  these  integrals  might  tend  to  dilute  the  effects  of 
nonequilibrium  on  these  macroscopic  parameters.  Thus  good  agreement  on 
these  macroscopic  parameters  should  be  expected  between  the  GLOW  and 
SHEATH  models.  This  was  found  not  to  be  the  case.  Table  V-2  summa¬ 
rizes  these  macroscopic  parameters  for  the  GLOW  and  SHEATH  models.  The 
voltage  and  thus  the  electric  field  at  the  cathode  are  determined  by 
the  net  multiplication  of  the  electrons  as  they  cross  the  cathode  fall. 
For  an  electronegative  gas,  equation  1 1 -2  for  the  multiplication  of  the 
electrons  becomes 

/ 

M*  exp  '1(x))c/k  ( V-l ) 

o 

From  the  results  of  the  nonequilibrium  analysis  in  Chapter  IV,  the 
attachment  coefficient  is  always  much  smaller  than  the  ionization 
coefficient  except  immediately  in  front  of  the  cathode.  The  attachment 
rate  for  the  cases  examined  was  only  a  small  perturbation  on  M  in 
equation  V-l.  This  also  agrees  with  the  results  presented  for  Xe/HCl 
mixtures  where  the  ionization  coefficient  effectively  did  not  change  as 
a  function  of  gas  mixture. 

Returning  to  Fig  V-4,  the  area  under  the  ionization  curves  are 
different  for  the  GLOW  and  SHEATH  models.  The  GLOW  model  assumed  the 
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TABLE  V-2 

Summary  of  Macroscopic  Parameters  from  the  GLOW  and  SHEATH  Models 
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electrons  were  in  equilibrium  with  the  field  and  results  in  the  elec¬ 
trons  being  able  to  ionize  immediately  upon  leaving  the  cathode.  The 
SHEATH  model  allows  the  electrons  to  not  be  in  equilibrium  with  the 
field.  As  a  result,  several  electron  collisions  occur  before  the 
majority  of  the  electrons  have  attained  the  energy  necessary  for 
ionization  of  the  gas.  Thus  the  nonequilibrium  model  predicts  a  higher 
cathode  fall  voltage  than  the  equilibrium  model  in  order  to  maintain 
the  same  net  multiplication. 

Analyzing  this  through  the  mathematical  relationships,  the  slope 
of  the  electric  field  is  determined  by  the  mobility  (k+)  of  the  posi¬ 
tive  ions  through  Poisson's  equation.  This  is  because  the  positive  ion 
current  density  at  the  cathode  is  determined  by  the  boundary  conditions 
at  the  cathode  involving  the  secondary  emission  coefficient  and  the 
total  discharge  current  density.  The  length  of  the  cathode  fall  is 
determined  when  the  multiplication  of  the  electrons  is  sufficient  to 
maintain  the  discharge  and  equals  1  +  l/y(equation  1 1 - 1 ) .  The  slope  of 
the  electric  field  and  its  extrapolated  x-axis  intercept  (dc)  then 
determines  the  voltage  for  the  cathode  fall.  i:us  the  difference  in 
the  macroscopic  parameters  in  Table  V-2  is  a  result  of  the  nonequi¬ 
librium  character  of  the  electrons  which  alters  the  ionization  coef¬ 
ficients  from  those  predicted  by  the  equilibrium  (GLOW)  model.  Thus 
any  agreement  of  the  equilibrium  method  with  experiment  is  entirely 
fortuitous  as  von  Engel  suggested  (ref  89:224). 

The  scaling  factor  SF(y,k+)  listed  in  Table  V-2  should  be  identi¬ 
cal  for  both  the  SHEATH  and  GLOW  models  since  they  both  used  the  same 
secondary  emission  coefficient  and  same  field  relationship  for  the  pos¬ 
itive  ion  drift  velocity.  The  slight  differences  in  the  data  between 


the  models  are  primarily  due  to  extrapolation  errors  in  determining  dc 

and  to  round  off  errors  in  d  and  V  . 

c  c 

Although  the  inability  of  the  GLOW  model  to  describe  nonequilibri¬ 
um  effects  accounts  for  most  of  the  disagreement  in  the  preceding 
figures  and  in  Table  V-2,  another  factor  affecting  the  agreement  of  the 
two  methods  is  the  use  of  two  different  data  bases.  The  GLOW  model 
bases  its  electron  kinetics  on  curve  fits  of  empirical  data  for  the 
Townsend  ionization  and  attachment  coefficients.  These  coefficients 
are  normally  derived  from  drift  tube  experiments  where  a  swarm  of  elec¬ 
trons  drift  between  two  electrodes  and  are  in  equilibrium  with  the 
field.  The  SHEATH  model  on  the  other  hand  bases  its  electron  kinetics 
on  scattering  cross  sections.  Under  ideal  conditions  both  determina¬ 
tions  of  a,  either  from  drift  tube  experiments  or  cross  sections, 
should  be  consistent.  Some  attempt  was  made  in  Appendix  C  to  correlate 
the  sets  of  cross  sections  used,  with  measured  transport  properties 
using  an  equilibrium  Boltzmann  code.  It  was  impossible,  however,  to 
correlate  them  over  the  whole  range,  since  transport  data  was  lacking 
at  some  high  E/p's.  In  addition  the  equilibrium  Boltzmann  code  was 
limited  to  calculations  at  low  and  medium  fields  because  it  did  not 
include  the  newly  created  secondary  electrons  in  the  distribution 
function.  Thus  the  data  bases  for  E/p's  less  than  about  160v/cm-torr 
differed  less  than  20%  but  for  higher  E/p  values  the  difference  could 
be  more. 

In  summary,  both  the  SHEATH  and  GLOW  models  predicted  the  same 
trends  in  He/HCl,  Ar/HCl ,  and  Xe/HCl  mixtures  and  both  models  scale 
according  to  equation  I 11-71 .  However,  the  discussion  in  this  chapter 
has  shown  that  the  GLOW  model  should  not  be  used  to  describe  any  of  the 


cathode  fall  parameters.  In  the  cathode  fall  region,  the  electrons 
have  been  shown  not  to  be  in  equilibrium  with  the  field.  As  a  result, 
the  cathode  fall  voltage  and  length  are  larger  in  the  SHEATH  model  than 
the  GLOW  model  in  order  to  compensate  for  the  nonionizing  region 
(Primary  Dark  Space)  near  the  cathode  where  the  electrons  are  being 
accelerated  up  to  the  ionization  energy  of  the  gas. 

This  investigation  of  the  cathode  fall  region  sets  the  stage  for 
further  investigations  of  the  negative  glow.  Due  to  the  slightly 
bimodal  character  of  the  electron  distribution  function  at  the  cathode 
fall -negative  glow  boundary  found  in  the  SHEATH  results,  a  nonequilib¬ 
rium  technique  should  also  be  used  to  investigate  the  negative  glow. 

In  this  region,  a  small  number  of  electrons,  that  have  managed  to  tra¬ 
verse  the  cathode  fall  with  one  or  no  collisions,  lose  their  energy  in 
the  negative  glow  forming  new  electrons  and  positive  ions  through 
ionization.  In  essence,  the  electron  distribution  function  is  still 
changing  faster  in  coordinate  space  than  it  is  in  velocity  space. 


Anode  Fall  Region 

This  section  discusses  the  applicability  of  both  methods  to  the 
sheath  region  at  the  opposite  end  of  the  discharge  near  the  anode.  In 
the  GLOW  model  calculations  the  anode  was  placed  after  a  chosen  length 
of  positive  column.  However,  due  to  a  limitation  of  memory  space  the 
anode  for  the  SHEATH  cases  was  placed  in  the  negative  glow.  This  cor¬ 
responds  to  an  "obstructed  discharge"  (ref  6).  The  properties  of  the 
anode  are  very  different  when  it  is  placed  in  a  slightly  negatively 
charged  region  such  as  the  Faraday  dark  space  as  compared  to  being  in 
contact  with  a  neutral  plasma  such  as  the  positive  column.  These  prop- 


erties  are  discussed  further  below. 

The  boundary  conditions  for  the  anode  fall  are  fairly  simple  in 
the  more  usual  case  when  the  anode  fall  and  positive  column  are  adja¬ 
cent.  Sufficient  electrons  must  be  accelerated  across  the  anode  fall 
and  create  sufficient  ionization  such  that  the  number  of  positive  ions 
crossing  the  anode  fall-positive  column  boundary  is  equal  to  the  number 
crossing  the  positive  column-Faraday  dark  space  boundary. 

It  is  commonly  believed  that  each  electron  entering  the  anode  fall 
must  be  accelerated  to  the  ionization  potential  of  the  gas.  Each  elec¬ 
tron  has  to  create  an  ion  before  striking  the  anode  in  order  to  main¬ 
tain  charge  neutrality  in  the  positive  column.  It  was  found  experi¬ 
mentally  (data  is  summarized  in  ref  1  and  6)  that  the  anode  fall  volt¬ 
age  was  usually  about  the  same  magnitude  as  the  ionization  potential, 
except  in  electronegative  gases.  However,  it  is  also  known  that  the 
electrons  in  the  positive  column  have  an  equilibrium  distribution  of 
energy,  with  a  few  having  an  energy  greater  than  the  energy  equivalent 
to  the  ionization  potential.  Thus  the  potential  is  different  for  each 
electron  to  be  accelerated  to  the  ionization  potential  of  the  gas.  In 
addition,  charge  neutrality  and  uniform  density  in  the  positive  column 
are  maintained  by  requiring  that  an  ion  enter  the  Faraday  dark  space 
each  time  one  enters  the  positive  column  from  the  anode  region. 
Similarly  an  electron  is  required  to  enter  the  positive  column  from  the 
Faraday  dark  space  each  time  one  leaves  the  positive  column  and  enters 
the  anode  fall  region.  Using  these  hypotheses,  it  is  not  necessary 
under  ideal  conditions  for  the  anode  fall  voltage  to  always  be  equal  to 
or  slightly  larger  than  the  ionization  potential. 

The  anode  acts  likes  a  probe  collecting  the  random  positive  and 
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negative  currents  from  the  positive  column  and  so  its  area  is  also  very 
important.  If  the  anode  area  is  so  small  that  the  negative  random  cur¬ 
rent  is  less  than  that  required  by  the  external  circuit,  then  a  region 
of  space  charge  forms  to  pull  more  negative  charges  from  the  positive 
column  and  make  the  collecting  area  of  the  anode  larger.  This  requires 
a  potential  on  the  order  of  the  ionization  potential  or  less,  since  a 
potential  much  larger  than  the  ionization  potential  would  tend  to 
create  numerous  electron-ion  pairs  reducing  the  effectiveness  of  the 
sheath.  Anode  fall  voltages  smaller  than  the  ionization  potential  have 
been  observed  (ref  6)  using  large  concave  anodes.  This  type  of  anode 
reduces  the  losses  to  diffusion  and  enables  more  negative  charges  to  be 
collected  from  the  positive  column. 

The  boundary  conditions  for  the  anode  fall  are  more  complicated 
when  the  anode  is  in  the  Faraday  dark  space.  Near  the  negative  glow 
boundary  there  may  be  sufficient  fast  electrons  from  the  cathode  fall 
to  carry  most  of  the  discharge  current.  These  electrons  are  already  of 
sufficient  energy  that  little  or  no  acceleration  is  required  for 
ionization  to  produce  positive  ions.  In  this  case  no  anode  fall  exists 
and  the  field  remains  almost  constant  right  up  to  the  anode.  As  the 
electrode  separation  is  increased,  the  anode  fall  region  must  of  course 
eventually  transition  into  the  normal  type  previously  described. 

Francis  (ref  6)  claims  a  normal  anode  fall  develops  suddenly  at  a  crit¬ 
ical  distance  from  the  cathode  lying  near  where  the  Faraday  dark 
space-positive  column  boundary  would  normally  occur.  This  critical 
distance  should  be  approximately  equal  to  the  stopping  distance  in  a 
gas  for  an  electron  with  an  energy  approximately  equal  to  the  cathode 
fall  voltage. 


Other  factors  described  by  Francis  (ref  6)  which  affect  the  volt¬ 
age  across  the  a node  fall  are: 

a)  stepwise  ionization.  Stepwise  ionization  results  in  a  smaller 
voltage  since  the  electrons  do  not  have  to  be  accelerated  to  the  full 
ionization  potential  in  order  to  create  an  ion. 

b)  the  electron  current  density.  This  is  closely  related  to  the 
anode  geometry.  Increasing  the  current  density  beyond  what  the  anode 
can  normally  carry  is  tantamount  to  decreasing  the  anode  area  and  hence 
increases  the  anode  voltage.  At  high  currents,  regularly  spaced  spots 
form  on  the  anode  which  are  regions  of  high  conductivity. 

c)  electronegative  gases.  Adding  an  electronegative  gas  to  an 
electropositive  gas  may  or  may  not  cause  the  anode  voltage  to  increase 
for  the  same  reasons  the  cathode  fall  may  or  may  not  contract  as  dis¬ 
cussed  in  previous  chapters.  The  potential  drop  across  the  anode  would 
depend  upon  the  partial  pressure  averaged  ionization  potential  of  the 
gas  mixture  in  addition  to  the  loss  of  electrons  to  negative  ion  forma¬ 
tion.  Generally,  a  high  coefficient  of  electron  attachment  greatly 
reduces  the  ionization  efficiency.  For  example,  anode  fall  voltages  in 
the  halogens  have  been  observed  to  be  several  times  the  ionization 
potential . 

d)  artificially  produced  space  charge.  Specially  coated  anodes 
which  yield  positive  ions  upon  electron  impact  obviously  have  smaller 
anode  voltage.  If  the  number  of  positive  ions  is  more  than  that 
required  for  the  flux  entering  the  positive  column,  the  fall  can  be 
negative. 

e)  proximity  of  side  walls.  Changing  the  distance  of  the  side 
walls  can  change  the  rate  of  diffusion  and  radial  fields  which  can 


result  in  increasing  or  decreasing  the  anode  fall  voltage. 

The  anode  fall  in  the  GLOW  code  calculations  was  adjacent  to  a 
positive  column  region.  The  anode  fall  voltage  in  Fig  1 1 1 -3 ,  1 1 1-7 , 
and  III-ll  were  in  the  range  of  l-3v.  These  values  were  determined  to 
be  the  difference  in  voltage  at  the  anode  and  at  the  minimum  electric 
field.  This  is  much  smaller  than  the  normal  anode  fall  voltage  which 
is  on  the  order  of  the  ionization  potential  of  the  gas.  The  use  of 
Ward's  formula  for  the  drift  velocity  (ref  92)  could  lead  to  a  more 
correct  anode  fall  voltage.  Since  the  electric  field  near  the  cathode 
fall-positive  column  boundary  is  less,  when  Ward's  formula  is  used, 
this  would  also  result  in  a  smaller  electron  current  density  entering 
the  positive  column.  A  broader  anode  fall  region  would  then  be  re¬ 
quired  for  the  electron  current  density  to  reach  the  total  discharge 
current  limit.  So  even  though  the  electric  field  gradient  would  also 
be  less  in  the  anode  fall  region,  the  fact  that  it  would  be  over  a 
longer  distance  could  increase  the  anode  fall  voltage. 

The  anode  in  the  SHEATH  calculations  was  placed  in  the  region  of 
the  negative  glow.  Any  correlation  in  the  data  close  to  the  anode  is 
indeterminant  due  to  the  changing  location  of  the  anode  in  the  negative 
glow  region.  This  occurred  because  the  cathode  fall  width  changed 
while  the  electrode  separation  was  held  constant  for  each  rare  gas 
mixture. 

In  addition,  the  SHEATH  code  changes  from  constant  energy  bins  to 
constant  spatial  bins  very  near  or  at  the  minimum  in  the  electric 
field.  At  present,  this  change  of  variable  is  not  smooth  for  all 
parameters  calculated  by  the  code,  so  c.iy  predictions  by  the  SHEATH 
method  on  the  anode  fall  will  be  held  in  reservation.  Thus  in  their 
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present  configuration,  neither  code  properly  describes  a  normal  anode 
fall  region. 

In  order  for  the  SHEATH  model  to  be  able  to  correctly  describe  a 
normal  anode  fall,  several  new  modifications  would  have  to  be  made. 

The  model  could  be  changed  to  reflect  the  boundary  conditions  at  the 
anode  instead  of  the  cathode,  so  that  the  calculation  started  at  the 
anode  and  proceded  towards  the  positive  column. 

J„(J  )*0  (V-2a ) 

J,U)  ♦  UJ)-  T  (V-2b) 

This  also  requires  that  E  (d),  j_(d),  and  the  electron  distribution 
striking  the  anode  be  known.  Potential  problems  with  the  point  at 
which  the  code  changes  from  equal  energy  increments  to  spatial  incre¬ 
ments  could  also  occur  again,  since  the  energy  gained  in  a  mean  free 
path  is  much  less  in  the  anode  fall.  Similarly  there  would  be  poten¬ 
tial  problems  with  greater  than  or  less  than  statements  when  the  sign 
convention  was  changed  from  x  to  -x. 

If  a  new  computer  system  is  installed  which  has  2-3  times  larger 
accessible  memory  space,  it  would  also  be  feasible  to  expand  the  array 
so  that  discharges  with  larger  electrode  separations  could  be  modeled. 
This  is  perhaps  the  easiest  solution  to  modeling  the  anode  fall  using 
the  SHEATH  code. 


Conclusions 

In  the  present  study,  two  very  different  methods  of  analyzing  glow 
discharges  (GLOW  and  SHEATH)  were  developed  and  applied  to  three  dif¬ 
ferent  cases  represented  by  He/HCl ,  Ar/HCl ,  and  Xe/HCl  mixtures.  Even 
though  the  two  methods  did  use  different  input  data  bases  founded  on 
the  prior  work  of  different  authors  using  different  experimental  ap¬ 
proaches,  the  results  of  the  two  methods  had  the  same  trends  and  both 
models  scaled  according  to  equation  I I 1-71 .  The  following  conclusions 
are  based  on  the  results  presented  in  Chapters  III,  IV,  and  the  compar¬ 
ison  of  the  models  in  Chapter  V: 

a)  Electrons  were  found  not  to  be  in  equilibrium  with  the  elec¬ 
tric  field  throughout  the  cathode  fall  region.  Close  to  the  cathode, 
the  whole  distribution  is  accelerated  up  to  an  energy  equal  to  the 
ionization  potential.  This  predominantly  dark  region  is  called  the 
Primary  dark  space.  Next  the  electron  distribution  function  becomes 
bimodal  as  some  of  the  electrons  lose  their  energy  to  ionization  and 
excitation.  This  luminous  region  is  called  the  cathode  glow.  In  this 
region,  the  majority  of  the  electrons  have  a  mean  energy  greater  than 
the  ionization  energy.  As  more  and  more  of  the  electrons  lose  energy 
to  inelastic  processes,  the  cathode  glow  merges  into  a  relatively 
darker  region  called  the  cathode  dark  space.  The  energy  difference 
between  the  decreasing  number  of  fast  and  the  increasing  number  of  slow 
electrons  gets  larger  as  the  fast  electrons  continue  to  be  accelerated 
by  the  linearly  decreasing  electric  field.  As  the  electric  field  be- 
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comes  very  small  and  almost  constant  in  the  negative  glow,  the  remain¬ 
ing  fast  electrons  resemble  an  electron  beam.  They  deposit  their 
energy  in  the  gas  crea+'ng  many  highly  excited  species  which  cascade 
down  and  produce  the  bright  luminescence  characteristic  of  the  negative 
glow  region.  As  a  result  of  this  nonequilibrium,  the  GLOW  or  equilib¬ 
rium  model  was  found  in  Chapter  V  to  be  inadequate  in  describing  both 
the  microscopic  parameters  such  as  charged  particle  densities  and  ioni¬ 
zation  and  attachment  coefficients  as  well  as  macroscopic  parameters 
such  as  cathode  fall  voltage,  cathode  fall  length,  electron  current 
density  and  electric  field.  A  model  which  allows  the  electrons  to  not 
be  in  equilibrium  with  the  electric  field  such  as  the  SHEATH  model  or 
Monte  Carlo  technique  must  be  used  to  describe  the  parameters  in  the 
cathode  fall  region.  The  negative  glow  region  should  also  be  modeled 
by  a  technique  which  allows  nonequilibrium.  The  model  for  the  negative 
glow  however  could  be  simplified  by  separating  the  electron  distribu¬ 
tion  function  into  a  low  energy  component  and  a  high  energy  component. 

b)  Better  agreement  between  theory  and  experiment  can  be  ob¬ 
tained  for  pure  gases  by  a  judicious  choice  of  y  and  k+.  Equation 
111-71  can  be  used  to  related  y  and  k+  to  experimental  cathode  fall 
voltage,  currents,  and  cathode  fall  widths.  These  values  of  y  and  k+ 
should  then  be  used  as  inputs  in  the  SHEATH  model  since  the  SHEATH 
model  was  shown  in  Chapter  IV  to  obey  the  scaling  relationship  given  by 
equation  1 1 1-72 . 

c)  The  slope  of  the  electric  field  in  the  cathode  fall  is  deter¬ 
mined  by  the  positive  ion  density.  If  Poisson's  equation  is  written  in 
terms  of  the  ratios  of  the  charged  particle  current  density  and  drift 
velocity,  then  the  slope  of  the  field  close  to  the  cathode  is  deter- 


mined  by  the  mobility  (k+)  of  the  positive  ions.  The  positive  ion  cur¬ 
rent  density  is  constant  close  to  the  cathode  and  is  determined  by  the 
total  discharge  current  density  and  y.  The  length  of  the  cathode  fall 
is  determined  by  the  maintenance  condition  for  the  discharge  involving 
the  multiplication  of  electrons  through  the  cathode  fall  (equations 
I I -1  and  V-l).  The  slope  of  the  electric  field  and  its  extrapolated 
x-axis  intercept  (dc)  then  determines  the  voltage  drop  across  the  cath¬ 
ode  fall  region. 

d)  The  formation  of  small  amounts  (n  <  n  )  of  negative  ions  in 
the  cathode  fall  region  was  found  not  to  be  the  reason  that  the  cathode 
sheath  length  contracts  in  the  axial  direction.  Actually,  the  forma¬ 
tion  of  negative  ions  causes  the  cathode  fall  voltage  to  increase 
slightly  which  in  turn  results  in  a  slight  increase  in  the  cathode  fall 
length.  These  trends  could  be  seen  in  the  results  of  both  the  equilib¬ 
rium  and  nonequilibrium  analyses  in  Chapters  III  and  IV. 

e)  When  an  attaching  gas  which  has  a  low  ionization  potential  is 
added  to  He  or  Ne,  the  observed  contraction  of  the  cathode  fall  length 
was  found  to  be  due  to  an  increase  in  ionization  in  this  region.  This 
increase  in  ionization  is  a  result  of  the  lowering  of  the  ionization 
threshold  and  the  difference  in  cross  section  shape  for  the  two  gases 
in  the  mixture. 

f)  The  SHEATH  method  which  allowed  the  electrons  to  not  be  in 
equilibrium  with  the  field  predicted  a  Townsend  attachment  coefficient 
which  contained  a  peak  at  the  cathode,  reached  a  minimum  at  the  region 
where  the  electron  mean  energy  and  Townsend  ionization  peaked  and  then 
grew  again  as  the  field  as  well  as  the  electron  mean  energy  decreased. 
Thus  the  cathode  fall  in  an  electronegative  gas  mixture  was  found  to  be 
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characterized  by  a  thin  region  close  to  the  cathode  in  which  attachment 
dominated,  followed  by  a  region  in  which  ionization  of  HC1  dominated 
the  ionization  of  He  or  A r  but  not  Xe,  before  uniform  ionization  of  the 
mixture  began. 

Although  this  analysis  was  mainly  concerned  with  modeling  phenome¬ 
na  in  the  cathode  fall  region,  the  conclusions  d  and  e  above  could  also 
be  applied  to  the  anode  fall  region.  In  contrast  to  the  cathode  fall 
region,  the  negative  charge  density  is  much  larger  than  the  positive 
charge  density  in  the  anode  fall  region.  Thus  the  width  of  the  anode 
fall  would  depend  on  how  fast  the  electrons  are  accelerated  to  the 
ionization  potential.  Thus  the  relationship  of  the  electron  drift 
velocity  to  the  field  is  more  important  in  this  region. 

Theoretical  Considerations  for  Future  Study 

In  the  course  of  this  study  several  areas  can  be  identified  which 
warrant  further  theoretical  investigation: 

a)  Since  ionization  of  HC1  caused  the  cathode  fall  to  contract 
in  He,  further  study  of  the  effects  resulting  from  changing  the  ioni¬ 
zation  cross-section  would  be  of  value.  Both  changes  in  the  threshold 
and  initial  slope  of  the  cross  section  should  be  examined  to  see  which 
affects  changes  in  the  cathode  sheath  structure  the  most. 

b)  A  parameter  mapping  of  the  secondary  emission  coefficient  and 
positive  ion  mobility  (k+)  is  recommended.  This  has  not  been  inves¬ 
tigated  using  nonequilibrium  techniques. 

c)  Using  the  present  results  from  the  cathode  fall  as  inputs, 
the  nonequilibrium  method  could  be  used  to  investigate  the  negative 
glow-Faraday  dark  space  region  provided  the  field  does  not  go  negative. 
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d)  The  linear  analysis  presented  in  this  study  can  be  expanded 
to  describe  the  kinetics  in  different  electrode  geometries  represented 
by  thyratrons,  plasma  reactors  for  depositing  thin  films,  or  hollow 
c'thode  discharges. 

e)  The  results  of  this  study  can  also  be  applied  to  understand¬ 
ing  the  field  variations  and  charged  particle  kinetics  in  striations. 
Different  boundary  conditions  would  be  required  for  striations.  The 
charged  particle  densities  on  each  side  of  the  striation  would  replace 
the  boundary  conditions  on  the  current  densities  at  the  electrodes. 

f)  Hot  cathodes  could  be  modeled  using  either  the  GLOW  or  SHEATH 
techniques.  Because  the  accelerating  voltage  is  less  in  the  cathode 
fall,  there  should  be  a  smaller  degree  of  nonequilibrium.  Thus  the 
GLOW  model  should  give  much  better  agreement  in  this  case.  In  addi¬ 
tion,  the  boundary  condition  at  the  cathode  relating  the  positive  ion 
to  the  electron  current  density  would  have  to  be  replaced  by  a  single 
source  of  electrons. 

g)  Due  to  limited  computer  storage  space,  the  only  electron  col¬ 
lision  processes  calculated  as  functions  of  distance  were  the  ioniza¬ 
tion  and  attachment  coefficients.  Further  investigation  of  the  other 
excitation  coefficients  and  coefficients  for  vibrational  excitation 
would  be  useful  to  complete  the  understanding  of  the  cathode  fall  re¬ 
gion. 

h)  Operation  of  the  SHEATH  code  would  be  much  more  efficient  as 
well  as  convenient  if  a  numerical  technique  were  added  to  converge  to 
the  lowest  value  of  Eq  for  which  a  complete  solution  exists,  before  the 
field  is  iterated  for  consistency  with  the  electron  distribution  func¬ 
tion.  This  technique  could  be  similar  to  the  modified  interpolation 


technique  which  was  used  in  the  GLOW  code. 

i)  More  detailed  ion  kinetics  options  would  make  the  SHEATH  model 
more  universal.  Electron-ion,  and  ion-ion  recombination,  detachment, 
and  two  step  ionization  could  be  added  such  that  the  ion  or  metastable 
densities  could  be  updated  after  each  calculation  of  the  electron  dis¬ 
tribution  function  for  a  self-consistent  solution.  Recombination  and 
detachment  will  be  more  important  in  the  negative  glow  and  Faraday  dark 
space.  These  processes  can  not  be  incorporated  in  the  present  GLOW 
model  because  they  are  nonlinear  with  respect  to  j  . 

j)  Cathodes  can  be  coated  with  special  films  which  emit  negative 
ions.  Further  investigation  of  the  effect  of  a  cathode  as  a  source  of 
negative  ions  on  the  cathode  fall  voltage  and  width  would  be  interest¬ 
ing.  A  further  increase  in  voltage  and  perhaps  width  of  the  cathode 
fall  should  be  anticipated  according  to  the  theory  in  Chapter  III  if 
negative  ions  are  released  at  the  expense  of  electrons. 

Experimental  Considerations  for  Future  Study 

Several  experimental  topics  have  also  been  identified  in  the 
course  of  this  study  which  warrant  further  investigation.  At  present 
the  first  two  topics  are  being  pursued. 

a)  No  experimental  evidence  could  be  found  for  the  non-contrac¬ 
tion  of  the  cathode  fall  in  mixtures  of  Ar  or  Xe  with  HC1  or  any  other 
attaching  gas.  A  small  discharge  tube  is  presently  under  construction 
in  the  Plasma  Physics  Group  (AFWAL/POOC-3)  which  will  allow  confirma¬ 
tion  of  this  phenomenon.  A  mass  spectrometer  will  be  used  to  monitor 


any  change  in  the  neutral  particle  chemistry. 

b)  It  is  now  possible  to  consider  measuring  the  electric  field  as 


a  function  of  distance  through  the  cathode  fall  region.  Some  potential 
techniques  are  non-perturbing  Doppler-free  laser  polarization,  Doppler- 
free  spectroscopy ,  or  Stark  induced  mixing  of  different  parity  levels 
in  an  excited  molecular  electronic  state. 

c)  The  energy  distribution  of  the  electrons  could  be  measured  in 
an  obstructed  discharge  by  sampling  the  electrons  through  a  small  hole 
in  the  anode.  These  results  could  then  be  compared  with  theory. 

d)  Similarly  the  energy  distribution  of  the  positive  ions  could 
be  sampled  through  a  small  hole  in  the  cathode.  Both  major  and  minor 
ions  should  be  sampled  because  of  the  differences  in  charge  exchange 
cross  sections.  By  knowing  their  energy  and  charge  transfer  rates, 
their  origin  in  the  cathode  fall  or  negative  glow  could  be  determined 
and  checked  with  theory. 

e)  Laser  induced  fluorescence  could  be  used  to  search  for  fast 
neutral  molecules  which  would  be  an  indication  that  charge  transfer  has 
occurred.  This  technique  would  compliment  those  described  in  c  and  d 
above  and  could  be  accomplished  simultaneously. 

It  is  hoped  that  this  study  has  provided  a  better  understanding  of 
the  cathode  fall  and  the  influence  of  negative  ions  in  this  region. 

The  generic  results  of  this  study  can  be  applied  to  excimer  lasers, 
thyratrons,  or  plasma  reactors  for  the  deposition  of  thin  films.  As  a 
final  guideline,  it  is  just  as  important  to  understand  the  positive  ion 
kinetics  in  the  cathode  sheath  region  as  it  is  to  understand  the  elec¬ 
tron  kinetics,  because  it  is  the  ion  mobility  which  determines  the 
slope  of  the  electric  field. 
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j_x  A:  Review  of  Neg 
Ion  Formation 


In  order  to  provide  continuity  in  the  descriptions  of  the  models, 
the  mechanisms  and  important  features  of  negative  ion  formation  are  re¬ 


viewed.  Several  reactions  occur  which 
listed  and  described  below: 

X  +  e”  X"  +  h  v 
X  +  e  +  M  -*•  X  +  M 
XY  +  e"  -  X  +  Y' 

XY  +  e"  -+  X+  +  Y~  +  e" 


form  negative  ions  and  these  are 

radiative  capture 
3  body  electron  capture 
dissociative  attachment 
ion  pair  formation 


(X  =  an  atom,  e"  =  an  electron,  hv  -  a  photon,  M  =  an  electron  or  an 
atom,  XY  =  a  molecule  made  of  atoms  X  and  Y) 


Dissociative  Attachment 

The  most  common  process  of  forming  negative  ions  is  dissociative 
attachment.  In  this  case,  an  electron  collides  with  a  molecule  XY  and 
dissociates  the  molecule  into  a  neutral  atom  X  and  a  negative  ion  Y~. 

To  aid  the  visualization  of  transitions  between  states,  dissociative 
attachment  can  be  divided  into  3  cases  using  molecular  potential  energy 
curves.  These  three  cases  will  be  discussed  next.  The  discussion  will 
include  the  influence  of  vibrational  excitation  on  dissociative 
attachment. 

In  all  three  cases  represented  in  Fig  A-la,  b,  c,  the  heavy  solid 
line  refers  to  the  ground  state  of  the  neutral  molecule  XY.  The 
nuclear  separation  oscillates  then  between  R^  and  R2.  The  thinner 
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Fig  A-l 

Potential  Energy  Curve  Allowing  Dissociative  Attachment 


curve  refers  to  the  molecular  ion  state  XY~  which  at  infinite 
separation  (X  +  Y")  tends  to  a  limit  which  lies  below  the  potential 
energy  curve  of  the  neutral  molecule  by  an  amount  of  energy  equal  to 
the  electron  affinity  (E  )  of  atom  Y.  The  solid  region  represents  the 

a 

range  over  vhich  the  negative  ion  state  is  stable.  The  dashed  region 
represents  the  range  over  which  the  negative  ion  state  is  unstable.  In 
this  region  the  negative  ion  may  give  up  the  electron  and  return  to  a 
ground  state  neutral  molecule  XY,  via  autodetachment.  If  v  ( R ) /f)  is  the 
rate  (sec-*)  for  autodetachment,  then  the  energy  in  this  unstable 
region  is  uncertain  by  an  amount  T  (R)  and  the  dashed  line  actually 
should  be  a  band  with  a  width  V(R).  The  molecular  ion  state  curve 
crosses  the  ground  state  curve  at  Rc>  Applying  the  Franck-Condon 
principle  to  the  transition  from  the  ground  molecular  state  to  the 
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Table  A-l.  Threshold  for  Dissociati ve 
Attachment  and  Ion  Pair  Production 


Molecule 


(eV) 

(ref.  8) 

y.oF-r- 


(eV)  (eV) 
(ref.  51) _ 


(eV) 

(ref.  51) 
10.43-Hg 


Fi  mole  FTh 


HgBr 

3.36-Br 

.71 

0. 

10.43-Hg 

7.78 

HgCl 

3 . 61 -Cl 

1.04 

0. 

10.43-Hg 

7.86 

F2 

3.45-F 

1.60 

0. 

17.42-F 

15.69 

15.57 

HgF 

3.45-F 

1.8 

0. 

10.43-Hg 

8.78 

IBr 

3.36-Br 

1.82 

0. 

10.44-1 

9.85 

8.89 

3.06-1 

0. 

11.84-Br 

10.59 

h 

3.06-1 

1.54 

0. 

10.44-1 

9.31 

8.92 

IC1 

3.61-C1 

2.15 

0. 

10.44-1 

10.08 

8.98 

3.06-1 

0. 

13.01-C1 

11.8 

Br^ 

3.36-Br 

1.971 

0. 

11.84-Br 

10.52 

10.45 

BrCl 

3.36-Br 

2.23 

0. 

11.84-C1 

11.1 

10.71 

3.61-C1 

0 

13.01-Br 

11.63 

FO 

3.45-F 

2.23 

0. 

13.61-0 

12.79 

12.39 

1.46-0 

.77 

17.42-F 

18.18 

ci2 

3.62-C1 

2.48 

0. 

13.01-C1 

11.5 

11.8 

BrO 

3.36-Br 

2.40 

0. 

13.61-0 

12.65 

1.47-0 

.S3 

11.84-Br 

12.  /  7 

IF 

3.45-F 

2.52 

0. 

10.44-1 

9.51 

3.06-1 

0. 

17.42-F 

16.58 

BrF 

3.45-F 

2.548 

0. 

11.84-Br 

11.78 

10.94 

3.36-Br 

0. 

17.42-F 

16.60 

CIO 

3.61 -Cl 

2.75 

0. 

13.61-0 

11.0 

12.76 

1.47-0 

1.29 

13.01-C1 

14.30 

HgH 

.  75-H 

.37 

0. 

10.43-Hg 

10.05 

Table  A-l.  (cont.) 


Ea  Do  EDA  Ei  Ei  mole  ETh 

(eV)2  (eV)2  (eV)  (eV)  (eV)  (eV) 

(ref.  8)  ref.  51)  (ref.  51) _ . _ 


Molecule 

(ref.  8) 

(ref.  51) 

(ref.  51) 

HI 

3.06-1 
. 754-H 

3.05 

0. 

2.79 

13.60-H 

10.44-1 

10.38 

13.59 

13.23 

10 

1.47-0 

3.06-1 

1.8 

.34 

0. 

10.44-1 

13.61-0 

10.78 

12.05 

HBr 

3.36-Br 
.  75-H 

3.76 

.40 

3.0 

13.60-H 

11.84-Br 

11.67 

13.99 

14.84 

HC1 

3.61-C1 

.75-H 

4.43 

.82 

3.68 

13.60-H 

13.1-C1 

12.75 

14.42 

16.69 

SH 

2.08-S 

.75-H 

3.55 

1.47 

2.80 

13 . 61 -H 
13.60-S 

10.43 

15.09 

16.39 

HF 

3.45-F 

.75-H 

5.94 

2.49 

5.18 

13.60-H 

17. 42- F 

16.06 

16.08 

22.6 

OH 

1.47-0 

.75-H 

4.45 

2.99 

3.70 

13.60- H 

13.61- 0 

12.91 

16.58 

17.31 

°2 

1.47-0 

5.12 

3.65 

13.61-0 

12.07 

17.26 

h2 

.75-H 

4.48 

3.73 

13.60-H 

15.43 

17.33 

HD 

.75-H 

4.51 

3.75 

13.60-D 

15.45 

13.35 

SO 

1.47-0 

2.08-S 

5.36 

3.89 

3.28 

10 . 36- S 
13.61-0 

10.29 

14.25 

16.90 

NO 

1.47-0 

1.5 

5.03 

14.54-N 

9.26 

19.57 

CN 

1.27-C 

7.76 

6.49 

14.54-N 

14.17 

21.03 

CO 

1.47-0 

1.27-C 

11.09 

9.63 

9.82 

11.26-C 

13.61-0 

14.01 

20.89 

23.44 

'DA 


Ei  mole 
ETh 
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Molecular  Ionization  Potential 
Threshold  Ion  Pair  Formation 
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negative  ion  state  leads  to  a  final  state  represented  by  some  point 
lying  between  a  and  b  on  the  potential  energy  curve.  The  threshold  for 


dissociative  attachment,  E^,  is  the  difference  between  the  molecular 
dissociation  energy  and  the  electron  affinity,  E  .  of  the  negative 
atomic  ion.  Thresholds  have  been  tabulated  for  several  molecules  and 
they  are 

EDA  =  °o-Ea  (fcl> 

compared  to  thresholds  for  ion  pair  formation  shown  in  Table  A-l. 

The  competition  between  autodetachment  and  dissociation  in 
molecular  negative  ions  leads  to  a  strong  dependence  of  the  cross 
section  on  the  mass  of  the  nuclei  and  the  rotational  and  vibrational 
states  of  the  neutral  molecule.  In  case  I  represented  in  Fig  la,  first 
suppose  autodetachment  did  not  occur.  Then  all  transitions  from  the 
ground  state  to  the  ion  state  would  result  in  dissociation  into  an  atom 
X  and  ion  Y"  with  a  relative  kinetic  energy  between  Ej  and  E 2.  If 
autodetachment  is  allowed  to  occur  at  some  internuclear  separation  R  at 
a  rate  V  ( R ) /f» ,  then  after  the  transition  to  the  negative  ion  state, 
the  atomic  nuclei  begin  to  separate  with  a  relative  velocity  V(R),  and 
the  fraction  of  negative  ion  molecules  remaining  (ref  8)  or  "surviving" 


(A-2) 


If  the  separation  between  the  atomic  nuclei  reaches  Rc  before  auto¬ 
detachment  takes  place,  then  the  neutral  molecule  dissociates  into 
X  +  Y“.  The  cross  section  for  this  process  can  then  be  approximated  by 


(A- 


t 


rJR)=cr€  dR 


where  qc  is  the  cross  section  for  the  initial  capture  (ref  8).  Because 

heavy  atomic  nuclei  typically  oscillate  more  slowly  than  lighter  nuclei, 

the  heavier  nuclei  spend  more  time  near  the  turning  points  and  have 

steeper  potential  wells.  This  increases  the  probability  that  the  inter- 

nuclear  separation  R  of  the  molecular  ion  will  exceed  Rc  before  it  auto- 

detaches.  Typical  peak  cross  section  values  from  the  ground  vibrational 
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state  for  dissociative  attachment  range  from  10  cm  for  light  mole- 
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cules  such  as  H2,  HD,  and  D2  to  10  or  10  cm  for  heavy  molecules 
such  as  HI,  and  I2«  (See  Table  A-2). 

In  case  II  represented  in  Fig  1  (b),  only  electrons  with  energies 
between  and  can  dissociate  the  molecule  and  produce  atomic  nega¬ 
tive  ions  through  the  dissociative  attachment  process.  The  resulting 
negative  ion  Y~  and  atom  X  thus  will  have  relative  final  energy  ranging 
from  0  to  Ej-Eg.  Electrons  with  energies  between  E2  and  E^  can  be  cap¬ 
tured,  but  they  will  form  excited  molecular  negative  ions  XY~.  These 
molecular  ions  will  eventually  autodetach  unless  the  excess  energy  can 
be  absorbed  through  a  collision.  Thus  the  total  attachment  cross  sec¬ 
tion  for  this  case  will  be  independent  of  pressure  for  pressures  greater 
than  some  critical  pressure,  but  will  be  proportional  to  the  pressure 
for  pressures  less  than  the  critical  pressure.  The  critical  pressure 
occurs  when  the  average  time  for  the  vibrational ly  excited  molecule  to 
transfer  its  excess  energy  to  a  gas  molecule  is  comparable  to  the  time 
for  autodetachment. 


Table  A-2.  (cont.) 

Process 

Diatomic 

Threshold 

Energy  At 

Molecules 

(eV)  Peak  (cm2) 

Peak  (eV) 

Reference 

Tri atomics 
e"  +  N20  0'+ 


e  +  OCS  S  +  . . . 


e  +  CS2  S  +  . 


e"  +  C02  0*  +  ... 


e~  +  H20  Neg  Ions 

e  +  H20  H  +  . . . 

e  +  C02  C  +  . . . 

Polyatomics 

e"  +  SFg  Neg  Ions 

e"  +  NF^  F”  +  . . . 

e"  +  NH3  Neg  Ions 

e  +  NFg  NF  2  +  . 

e  +  NF^  F  2  +  . . 

e  +  CF,  CF  ~  +  . 
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In  case  III  represented  in  Fig  1  (c),  negative  ion  molecules  with 
vibrational  states  v<v'  are  stable,  but  those  with  v>v'  are  unstable 
and  can  result  in  autodetachment.  These  states  may  be  stabilized 
through  collisions  similar  to  the  process  described  for  case  II. 
Similarly,  the  cross  section  for  attachment  will  exhibit  resonance 
maxima  which  are  associated  with  the  vibrational  levels  with  v>v' 
provided  the  lifetime  of  the  excited  vibrational  state  is  long  compared 
to  its  period.  This  situation  occurs  in  02  and  NO,  resulting  in  the 
production  of  0~2  and  NO"  after  collisions  with  other  gas  molecules. 
Table  A-2  tabulates  experimental  data  for  thresholds  and  cross  section 
peaks  for  several  attaching  molecules. 

Increased  vibrational  excitation  enhances  the  dissociative 
attachment  cross  section  for  all  three  cases  since  electron  capture  can 
generally  occur  at  larger  and  larger  values  of  R.  Thus  the  molecular 
negative  ion  is  more  likely  to  "survive",  the  closer  R  approaches  Rc 
through  vibrational  excitation.  Rotational  excitation  can  also  enhance 
the  cross  section  for  dissociative  attachment  due  to  the  centrifugal 
stretching  occurring  in  the  excited  states,  but  in  a  less  dramatic 
manner  than  vibrational  excitation.  Bardsley  and  Wadehra  (ref  17,  91) 
have  shown  that  the  magnitude  of  the  dissociative  attachment  cross 
section  from  the  v=  5  level  in  H2  and  D2  increases  three  to  four 
orders  of  magnitude.  Allan  and  Wong  (ref  12)  have  measured  similar 
Increases  in  their  experiments  with  HC1.  The  change  in  the  "survival 
factor"  is  thought  (ref  91)  to  be  more  important  than  the  change  in  the 
Franck  Condon  factor  in  explaining  this  dramatic  effect.  If  the 
lifetime  of  the  vibrational ly  excited  molecule  is  long  compared  to  the 
period  of  vibration,  the  attachment  probability  will  contain  resonance 


peaks.  This  is  a  result  of  the  initial  electron  capture  cross  section 


o  exhibiting  resonance  peaks  due  to  the  ratios  of  populations  in  the 
various  vibrational  levels.  Experimentally  this  variation  with  respect 
to  rotational  and  vibrational  population  appears  as  a  temperature 
dependence. 

The  remaining  processes,  including  radiative  capture,  three  body 
capture,  and  ion  pair  formation  are  generally  less  important  than 
dissociative  attachment  and  will  be  discussed  only  briefly. 

Radiative  Capture: 

Radiative  capture  is  the  simplest  process  for  creating  negative 
ions.  Electrons  are  captured  directly  by  neutral  atoms  with  the  excess 
energy  being  dissipated  in  radiation.  The  amount  of  energy  released  by 
this  process  is  equal  to  E  +  E.  where  E  is  the  kinetic  energy  of  the 

O 

Incident  electron  and  E_  is  the  electron  affinity  of  the  atom.  This 

H 

capture  process  produces  a  continuous  emission  spectrum  extending 
Indefinitely  from  the  long  wavelength  limit  at 


-22  2 

A  typical  magnitude  for  this  process  is  on  the  order  of  10  cm  with 
a  large  peak  near  zero  impact  energy  (ref  8).  This  is  two  orders  of 
magnitude  smaller  than  the  smallest  cross  section  for  dissociative 
attachment.  The  relatively  large  low  energy  peak  occurs  because  the 
electron  spends  a  longer  time  within  the  range  of  the  atomic  field  at 
lower  incident  energies. 


Three  Body  Electron  Capture: 

In  three  body  electron  capture  the  energy  released  by  the  process 

is  converted,  not  into  a  photon,  but  into  kinetic  energy  of  the  third 

body.  The  magnitude  of  the  attachment  cross  section  for  this  process 

can  be  determined  by  arguments  based  on  two  body  collisions.  The  third 

body  will  be  effective  in  absorbing  energy  only  if  it  makes  an 

effective  collision  with  the  capturing  atom,  when  the  electron  is  also 

-8 

within  an  atomic  radius  (r  «  10"  cm)  of  the  capturing  atom.  If  X  is 
the  mean  free  path  for  effective  collisions  between  the  atom  and  a 
single  third  body,  then  the  probability  for  an  effective  collision  is 
to  a  good  approximation  rQ/X  .  This  results  in  an  effective  three  body 
attachment  cross  section  that  can  be  approximated  by  the  following 
expression  in  terms  of  an  effective  cross  section  per  unit 
concentration  of  the  third  body. 


°r  “  (A-5) 

Thus  it  remains  to  determine  X  ,  which  varies  considerably  depending  on 
the  nature  of  the  third  body.  Two  cases  will  be  considered: 

1)  when  electrons  are  third  bodies,  and 

2)  when  atoms  and  molecules  are  third  bodies. 

When  electrons  are  the  third  body,  the  effective  cross  section  for 
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energy  transfer  is  on  the  order  10  cm  ,  and  the  mean  free  path  for  a 
unitary  third  body  density  simply  becomes 

X  =  1016cm  (A-6) 

and  the  third  body  attachment  cross  section  is  on  the  order  of 

o  =  10" 

a 
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In  comparison  to  dissociative  attachment,  an  electron  number  density 

20  3 

greater  than  or  equal  to  10  cm  would  be  required  to  produce  a 
comparable  number  of  negative  ions  per  sec.  For  most  purposes  this  is 
an  exceptionally  large  electron  density. 

When  atoms  or  molecules  are  the  third  body,  two  cases  exist: 
resonant  and  non-resonant  energy  transfer.  In  the  resonant  case,  the 
energy  transferred  is  absorbed  totally  as  internal  energy  in  the  third 
body.  The  cross  section  for  this  process  is  known  to  be  quite  large  on 
the  order  10”*4cm3  or  greater(ref  8).  This  yields  a  mean  free  path  and 
effective  attachment  cross  section  of 

A  =  1014cm  a  a  =  10‘38cm5  (A-8) 

O 
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Thus  third  body  densities  of  10  cm  or  larger  would  be  necessary  to 
make  three  body  attachment  comparable  to  that  of  dissociative 
attachment.  Thus  molecules  possessing  several  internal  degrees  of 
freedom  are  the  most  effective  as  third  bodies.  However,  energy 
transfer  into  vibrational  motion  often  is  not  as  efficient  as  transfer 
into  rotational  motion. 

Ion  Pair  Formation: 

The  last  process  to  be  discussed  which  forms  negative  ions  is  ion 
pair  formation  or  polar  excitation.  During  this  type  of  collision,  an 
electron  excites  a  molecule  (XY)  to  an  unstable  electronic  state  which 
dissociates  spontaneously. 

XY  +  e"  -*■  X+  +  Y"  +  e“  (A-9) 

This  process  differs  from  the  previous  processes  in  that  the  electron  is 
not  captured  during  the  collision,  but  corresponds  to  an  inelastic 
collision  in  which  the  electron  acts  as  a  source  of  energy  for 


electronic  excitation.  Thus  polar  excitation  occurs  at  much  higher 
electron  energies  than  the  capture  processes.  Minimum  threshold 
energies,  usually  range  from  10  to  25  eV. 

A  typical  potential  energy  curve  for  this  process  is  depicted  in 
Fig  A-2. 

I  * 

§ 

o 

R,  R, 
distance 

Fig  A-2.  Potential  Energy  Curve  For  Ion  Pair  Formation 

Again  the  heavy  solid  line  refers  to  the  ground  state  of  the  neutral 
molecule.  The  thinner  line  refers  to  an  unstable  electronic  state  which 
dissociates  spontaneously  into  ions.  The  energy  at  infinite  nuclear 
separation  tends  to  a  value  of  E.  -  E  .  where  E.  is  the  ionization 

id  1 

energy  of  X  and  Ea  is  the  electron  affinity  of  Y.  The  final  relative 

d 

kinetic  energy  (T)  of  the  ions  is  given  by  the  following  equation: 

T  ■  ETh  *  Ea  '  Ei  -  Do  (A-10) 

where  Do  is  the  dissociation  energy  of  the  ground  state.  Thresholds  for 

ion  pair  formation  for  several  diatomics  are  previously  listed  in  Table 

A-l. 

The  shape  of  the  cross  section  for  polar  dissociation  differs  from 
the  previous  capture  process.  Actually  it  resembles  the  shape  of  other 
electronic  excitation  processes  and  thus  is  typified  by  a  sharp  rise 


near  threshold  to  a  maximum  on  the  order  of  10  -10"A/cnr  at  2  to  4 
times  the  threshold  energy,  after  which  it  declines.  The  decline  being 
inversely  proportional  to  the  energy  or  more  slowly  if  optical 
transitions  are  allowed  between  the  ground  state  and  the  excited  state 
(ref  68:497). 
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Appendix  B.  Parameters  for  GLOW 
Code  Calculations 

This  appendix  describes  the  parameters  used  to  characterize  the 
ionization,  attachment,  diffusion,  and  ion-ion  recombination 
coefficients  as  well  as  the  ion  drift  velocity.  Most  of  t1  data  for 
the  rare  gases  was  taken  from  the  literature.  Effort  was  2n  to 
obtain  the  best  parameters  for  HC1  from  drift  tube  measures  ■r  of 
Davies  (ref  33:31-55). 

The  parameters  for  the  rare  gases  are  summarized  in  Table  B-l.  A, 
B,  and  s  refer  to  the  parameters  in  equation  1 11-38  and  were  taken  from 
Ward  (ref  95:2790).  The  transverse  diffusion  coefficient  was  obtained 
from  data  tabularized  by  Dutton  (ref  39:6491  (He),  651  (Ar),  654  (Xe)). 
The  mean  free  path  X  was  estimated  using  the  formula: 


where  o  is  the  momentum  transfer  cross-section.  Momentum  transfer 
cross-sections  have  been  tabulated  by  Kieffer  (ref  54:2(He),  13  (Ar),  17 


Table  B-l 


Parameters 

Used  to  Describe  Rare  Gas 

Processes 

in  the  GLOW  Code 

A 

B  s 

De» 

X 

r  k+ 

1 

Volt* 

cm2torr 

cm 

cm! torr* 

cm  torr 

cirrtorr2 

sec 

volt*sec 

(Xe)).  The  value  of  the  momentum  transfer  cross  section  at  leV  was  used. 
This  energy  is  representative  of  electrons  diffusing  radially  in  the 
positive  column  region  of  a  discharge.  The  value  of  r  was  also  taken 
from  the  data  given  by  Ward  (ref  92:2791)  for  v+  and  vg. 

The  ion-ion  recombination  rates  were  estimated  using  the  approxi¬ 
mate  scaling  formula  developed  by  Hickman  (ref  49:4875).  The  formula  is 
given  by 

l  -  _ t _  ( B-2 ) 

"  mi  (£./ I.)* 

where  T  is  the  gas  temperature  in  degress  Kelvin,  m  is  the  reduced  mass 
of  the  diatomic  molecule  in  a.u.,  and  E.A.  is  the  electron  affinity  of 
the  negative  ion  in  eV.  The  values  of  these  parameters  and  the  result¬ 
ing  ion-ion  recombination  rate  for  a  300°K  gas  is  given  in  Table  B-2. 


The  ionization  and  attachment  coefficients  for  HC1  were  curve  fit 
from  the  experimental  data  of  Davies  (33:33).  Using  equation  1 1 1-38 , 
both  values  of  s  =  i  and  s  =  1  were  tried.  The  value  s  =  1  gave  the 
best  fit,  which  also  yielded  values  for  A  =  14.34  and  B  =  243.94.  In 
order  to  model  the  attachment  coefficient  over  the  complete  range,  both 


equations  I I 1-39  and  I I 1-40  had  to  be  used.  A  fifth  order  (n=5)  poly¬ 
nomial  was  used  for  E/p  <  90v/cm-torr  with  the  following  values  for 
coefficients 

a.  =  -.748878814, 
a?  =  .1274582651, 
ai  =  -.0052731662, 
ai  =  9.48227228E-5, 
a.  =  -7.86191615E-7,  and 
a*  =  2.46552683E-9. 

For  E/p  >.  90v/cm-torr,  the  best  fit  was  given  by  C  =  .254193623029, 
D  =  9.14964650482E-3,  and  t=l.  This  concludes  the  data  inputs  used  for 
the  GLOW  code  except  for  the  discharge  parameters,  such  as  current, 
pressure,  and  electrode  separation  which  are  discussed  in  the  Chapter 
III  with  the  results  of  the  calculations. 
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This  appendix  describes  the  parameters  used  in  the  SHEATH  code. 

The  initial  values  of  code  parameters  will  be  summarized  first  followed 
by  a  description  of  the  electron  impact  cross-sections  used  in  calcula¬ 
ting  the  electron  distribution  function. 

The  discharge  parameters  are  summarized  in  Table  C-l.  The  anode 
was  placed  just  outside  the  cathode  fall  region  where  the  electric  field 
leveled  off  to  a  constant  value.  This  was  done  to  keep  the  energy  step 
size  as  small  as  possible,  less  than  1.6eV,  in  order  to  maximize  the 
energy  resolution  in  calculating  the  distribution  function.  With  this 
step  size  limitation,  lengthening  the  electrode  distance  was  not  feasi¬ 
ble  since  most  of  the  storage  space  on  a  Cyber  74  was  already  being 
used.  For  100%  He,  the  electrode  distance  was  2.2  cm  and  for  the  re¬ 
maining  He/HCl  mixtures  the  electrode  distance  was  2.1  cm  due  to  the 
contraction  that  accurred.  J  Is  the  discharge  current  density  and  y  is 
the  secondary  emission  coefficient.  The  distribution  of  electrons 

Table  C-l 

Discharge  Parameters  Used  in  the  SHEATH  Code 

D  J  Y 

Mixtures _ (cm) _ (  amp/cm2) _ 


v\ 


leaving  the  cathode  was  the  same  as  that  used  by  Long  (ref  60)  and  is 
displayed  in  Fig  C-l.  The  x  axis  is  plotted  in  terms  of  energy  bins 
which  were  typically  1.  -  1.6eV  in  width  depending  on  the  gas  mixture. 
This  distribution  is  higher  in  energy  than  the  one  used  by  Boeuf  et  al. 
(ref  14)  which  was  a  flat  isotropic  energy  distribution  over  the  range 
of  0  to  5eV.  Boeuf  and  Marode  (ref  13)  later  used  a  narrower  dis¬ 
tribution  which  was  isotropic  and  flat  between  4  and  5eV. 

The  gas  parameters  are  summarized  in  Table  C-2.  The  gas  density 
was  constant  throughtout  the  discharge  and  was  determined  from  the  gas 
pressure,  pQ  and  temperature,  Tq.  The  positive  ion  mobility,  k+  was 
obtained  from  Ward  (ref  92:2791)  and  was  the  same  as  that  used  in  the 
equilibrium  analysis.  determined  the  point  at  which  the  code 

transforms  from  energy  space  to  spatial  coordinate  as  described  earlier. 
It  was  varied  so  the  coordinate  transition  was  as  close  as  possible  to 
the  minimum  of  the  electric  field  within  the  bounds  of  equation  IV-11. 


Table  C-2 

Gas  Parameters  Used  in  the  SHEATH  Code 


Mixtures 

Po 

(torr) 

To 

(°K) 

k+ 

cm  -torr 

volt  sec 

He/HCl 

1 

300 

4.1E+4 

Ar/HCl 

1 

300 

8.25E+3 

Xe/HCl 

1 

300 

4E+3 

The  cross 

sections  for  He,  Ar, 

Xe,  and  HC1 

which  were  used  in  the 

SHEATH  code  calculations  are  discussed  next.  The  most  important  cross 
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sections  In  the  calculations  are  the  ionization  cross  sections  of  the 
rare  gases.  This  cross  section  determines  how  quickly  the  electrons 
multiply,  and  thus  directly  effects  ng  and  dE/dx  in  Poisson's  equation. 
For  the  rare  gases,  the  ionization  cross  section  was  obtained  from 
measurements  by  Rapp  et  al.  (ref  75)  over  the  complete  range.  The 
ionization  cross  section  for  HC1  had  to  be  extended  above  lOOeV.  Since 
HC1  was  used  in  small  concentrations  and  the  mean  electron  energy  in 
the  HC1  mixtures  was  always  less  than  87  eV,  any  error  in  the  approxi¬ 
mation  should  be  negligible.  More  specifically  the  peak  electron  mean 
energy  ranged  in  He  mixtures  from  59eV  to  99eV  in  pure  He,  for  Ar  mix¬ 
tures  it  ranged  from  20.5eV  to  21. eV,  and  for  Xe  mixtures  it  ranged 
from  25.4eV  to  25.9 eV. 

The  six  cross  sections  used  to  model  the  electron  kinetics  in  He 
are  plotted  in  Fig  C-2.  The  threshold  and  references  for  each  range  are 
listed  in  Table  C-3  below.  The  electronic  cross  sections  were  extended 
using  the  formulas  listed  under  the  reference  column. 

The  cross  section  for  Ar  are  plotted  in  Fig  C-3  and  the  references 
and  thresholds  are  given  in  Table  C-4.  The  three  electronic  cross  sec¬ 
tions  were  derived  using  a  standard  technique  of  backing  out  cross  sec¬ 
tion  from  transport  data  using  a  code  which  solves  the  equilibrium 
Boltzmann  equation.  Ratios  between  these  cross  sections  were  obtained 
from  data  of  Eggarter  (ref  41)),  Ganas,  eta!  (ref  44),  and  Chutjan 
(ref  28).  Keeping  the  ratios  constant,  the  sum  of  the  three  cross  sec¬ 
tions  was  then  varied  and  the  predicted  electron  drift  velocity  and 
Townsend  ionization  rate  were  compared  with  transport  data  listed  in 
Dutton  (ref  39)  and  Specht,  et  al.  (ref  83).  Up  to  20eV,  the  agreement 
with  transport  data  is  within  9%  and  the  total  electronic  cross  section 
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Table  C-3 


Electron  Impac*  Cross  Sections  In  He 


Process 


Momentum 

Transfer 

Electronic 

23S  +  2*S 


Ionization 


Threshold 


0. 

19.805 

20.949 

21.203 

23.071 

24.586 


Range 


0.-6. 

6.-500. 


0.-23.54 

24.-100 

100.-500. 

0.-199.6 

200.-500. 

0.-196.1 

200.-500. 

0.-207.8 

210.-500. 

0.-500. 


References 


Kieffer  (ref  54) 
Hayashi  (ref  48) 


Kieffer  (ref  54) 
Constant  . 
o=  .624/fe* 

Kieffer  (ref.54) 
o-  .02755A4 

Kieffer  (ref  54) 
a~  1.179/e:  * 

Kieffer  (ref  54) 
o  =  .3605£* 

Rapp,  et  al.  (ref  75) 


Process 


Table  C-4 

Electron  Impact  Cross  Sections  in  Ar 


Threshold 


Range 


Reference 


agrees  well  with  Jacob  and  Mangano  (ref  53)  up  to  17eV  where  their 


cross 

section  stopped.  The  shape  of  the  total  electronic  cross  section  for 
higher  cross  energies  was  obtained  from  de  Heer,  et  al.  (ref  35). 

The  cross  sections  for  Xe  are  plotted  in  Fig  C-4  and  the  refer¬ 
ences  and  thresholds  are  given  in  Table  C-5.  As  for  Ar,  the  electronic 
and  a  section  of  momentum  cross  section  for  Xe  were  derived  by  backing 
them  out  of  transport  data. 


Table  C-5 

Electron  Impact  Cross  Sections  in  Xe 


Process 

Threshold 

(eV) 

Range 

(ev) 

Reference 

Momentum 

0. 

o 

• 

1 

<T\ 

• 

CJ1 

Kieffer  (ref  54) 

Transfer 

6.5-20. 

Derived 

20.-500. 

Hayashi  (ref  48) 

Excitation 

8.32 

o 

• 

1 

CJ1 

o 

o 

• 

Derived 

Ionization 

12.13 

o 

• 

1 

CJ1 

o 

o 

• 

Rapp,  et  al.  (ref  75) 

The  momentum  transfer  cross  section  was  lowered  in  the  region  from 
about  6.5  to  20eV  in  order  to  get  the  drift  velocity  to  agree  within 
20%  of  that  measured  by  Huang  (ref  50)  and  the  data  listed  in  Dutton 
(ref  39).  Huang  found  that  the  electron  mobility  for  fields  above  4Td 
remained  constant,  given  by  yN  *  5.5E21  molecules/(v-sec-cm)*.  The 
shape  of  the  electronic  cross  section  below  lleV  was  taken  from 
Dowell's  unnormalized  data  (ref  37).  The  ionization  rate  predicted 
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from  the  Boltzmann  code  agreed  within  20%  with  the  data  in  Dutton  (ref 
39)  up  to  60Td  using  this  electronic  cross  section  up  to  20eV.  The 
shape  at  higher  energies  was  obtained  from  Ganas,  et  al.  (ref  44)  and 
de  Heer  et  al.  (ref  35)  and  normalized  to  the  value  derived  for  20eV. 

The  cross  sections  for  HC1  are  plotted  in  Fig  C-5  and  Fig  C-6.  The 
thresholds  and  references  for  each  range  are  listed  in  Table  C-6.  The 
electronic  cross  section  for  the  Band  C  state  was  assumed  to  have  a  peak 
at  lOOeV  which  is  slightly  less  than  the  energy  of  the  peak  of  the  ioni¬ 
zation  cross  section.  The  cross  section  was  assumed  to  decrease  accord¬ 
ing  to  the  formula  a  -  18./ e^.  The  shape  of  the  ionization  cross 
section  from  Compton,  et  al.  (ref  79)  was  used  over  the  range  100-280eV 
and  normalized  to  Davies'  data.  At  higher  energies  the  shape  of  the 
ionization  cross  section  for  Og  (ref  75)  was  used  which  also  peaks  at 
the  same  energy.  It  was  normalized  to  the  HC1  ionizaton  at  280eV.  This 
concludes  the  discussion  of  the  data  inputs  used  for  the  SHEATH  code. 

Table  C-6 

Electron  Impact  Cross  Sections  in  HC1 


Process 

Threshold  Range 

References 

(eV)  (eV) 

Momentum  CT  O.-lOO.  Davies  (ref  33:45) 


Transfer  100.-500.  Logrithrimic  extrapolation 


Vibration 

v  =  1 

.36 

0.-7. 

Davies  (ref  33:46) 

v  =  2 

.70 

0.-7. 2 

Davies  (ref  33:47) 

Attachment 

Cl" 

.67 

0.-2. 7 

Davies  (ref  33:47) 

H“ 

5.6 

0.-12. 

Davies  (ref  33:48) 

Excitation 

A 

5.5 

0.-11. 

Davies  (ref  33:47) 

B  +  C 

9.3 

0.-100. 

100.-500. 

Davies  (ref  33:48) 
Derived 

Ionization 


12.74 


0.-100. 

100.-500. 


Davies  (ref  33:48) 
Derived 


FIG  C-5  HIGH  ENERGY  ELECTRON  CROSS  SECTIONS  IN  HC1 
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